SOME ASPECTS OF 
NEUTRINO MIXING AND OSCILLATIONS 



o 

O ■ THESIS SUBMITTED TO 



(N 



X 



THE UNIVERSITY OF CALCUTTA 
FOR THE DEGREE OF 



<. 

. DOCTOR OF PHILOSOPHY (SCIENCE) 

(N 

(N ■ By 
> 

^ SANJIB KUMAR AGARWALLA 

(N 

■ DEPARTMENT OF PHYSICS, 

oo ■ 

O : UNIVERSITY OF CALCUTTA 

o 



SEPTEMBER 2008 



ii 



- To My Parents - 



iii 



iv 



Acknowledgements 



First and foremost, I would like to express my deepest respect and most sincere gratitude to my 
supervisor, Professor Amitava Raychaudhuri, for introducing me into the fascinating world of 
neutrinos and for being an excellent teacher and a perfect guide. I convey my regards to him not 
only for his continuous support and patience throughout this study, which was a great comfort, 
but also for the new dimensions he brought into my life. Without his inspiring discussions, 
invaluable guidance and exact encouragement, this accomplishment would have never been 
fulfilled. For almost everything I know about neutrino phenomenology, I owe to him. I consider 
myself very fortunate to have him as my supervisor who gave me immense freedom to pursue 
my interests and always encouraged me to take part in various summer schools, workshops, 
conferences in India and abroad from the very starting point of my research career. This was 
quite helpful in building my confidence to discuss physics issues with other scientists which most 
of the times culminated into fruitful collaborations. 1 am greatly indebted to his assistance and 
understanding in matters of non-academic concern which have helped me endure some difficult 
times during my study period. His suggestions will always help me in moving the things ahead 
in my future research career. 

There are no words to describe such a great honour that I have the chance to convey my sincere 
gratefulness to my senior collaborator Dr. Sandhya Choubey with whom I have done most of 
my work. She has been my constant source of inspiration for gaining grounds in research work 
and I strongly feel indebted to her as nothing would have been possible without her support. 
I would like to thank her for her close daily guidance, her patience in assisting me during the 
course of this work, her kindness, and her opportune encouragements. Her commitment to 
research and dedication to physics, and even her unique style of language and her sense of 
humour have all greatly influenced me. 

I am very much thankful to Dr. Srubabati Goswami, my senior collaborator, with whom a part 
of the research work contained in this thesis was done and whose extraordinary energy and 
enthusiasm for physics has set a great example for me. 

I would like to thank Dr. Abhijit Samanta for my first successful collaboration and for helping 
me regarding the use of computers. He also shared his expertise on various issues related to 
the ICAL@1N0 detector simulation with me. 

Mention must be made of my senior research colleague Dr. Rathin Adhikari, who extended a 
helping hand both in academic and non-academic issues whenever it was necessary. 

I would like to thank Dr. Subhendu Rakshit for successful collaboration and providing a friendly 
environment for discussing physics. 

I express my gratitude to Dr. Walter Winter and Dr. Patrick Huber for providing me useful 
informations regarding the simulations related to long baseline neutrino experiments at various 
stages of this work. 



V 



I would like to express my sincere gratitude and appreciation to Professor Mina Ketan Parida 
who introduced me to the topic of renormalization group evolution analysis during his stay 
at Harish-Chandra Research Institute and provided me an opportunity to collaborate with 
two big names in the field of neutrino physics - Professor R. N. Mohapatra and Professor G. 
Rajasekaran. 

The constructive suggestions and useful comments from Professor Manfred Lindner, Professor 
A. Smirnov, Professor G. Raffclt, Professor Andre de Gouvea, Professor Mats Lindroos, Profes- 
sor Mauro Mezzetto, Professor J. J. Gomez Cadcnas, Dr. Andrea Donini, Dr. Pilar Hernandez, 
Dr. Osamu Yasuda and Dr. Francesco Terranova helped me a lot at different stages of this 
work. 

I gratefully acknowledge generous helps offered by Professor Naba Kumar Mondal who is the 
spokesperson of the INO collaboration. At this stage I would like to extend my thanks to 
Professor Kamales Kar, Professor Sudeb Bhattacharya, Professor G. Rajasekaran, Professor 
H. S. Mani, Professor M. V. N. Murthy, Professor Raj Gandhi, Dr. Amol Dighe, Dr. D. 
Indumathi, Dr. Brajcsh Choudhary, Dr. Gobinda Majumder, Dr. Debasish Majumdar, Dr. 
Ambar Ghosal, Dr. Satyajit Saha, Dr. Subhasish Chattopadhyay, Dr. Abhijit Bandhopadhyay 
for having useful discussions on INO. I also acknowledge the moral support from my friends 
Saikat Biswas, Tapasi Ghosh, Satyajit Jena who are working for the INO collaboration. 

I benefitted from discussing physics with Professor Amitava Datta, Professor Ashoke Sen, 
Professor Rajesh Gopakumar, Professor Dileep Jatkar, Professor Biswarup Mukhopadhyaya, 
Professor Soumitra Sengupta, Professor Anjan S. Joshipura, Professor Sreerup Raychaudhuri, 
Professor Gautam Bhattacharyya, Dr. Subhendra Mohanty, Dr. C. S. Aulakh, Dr. Anindya 
Datta, Dr. Asesh Krishna Datta and Dr. Raghavan Rangarajan at various stages. 

I am much obliged to Professor Tapan Kumar Das, Professor Subinay Dasgupta, Dr. Anirban 
Kundu, Dr. Parongama Sen and Dr. Gautam Gangopadhyay from Calcutta University who 
have helped me in various ways during this work. 

I wish to thank Dr. Arunansu Sil who helped me at the initial stage of my research work to 
grab the basics of neutrino oscillations. I would also like to thank my friend Shashank Shalgar 
for sharing his insights in neutrino physics and computer related issues. I admire the moral 
support that I received from Swarup Majee. 

My heartiest thanks to Soumitra Nandi, Biplob Bhattacharjcc, Kamalika Basu Hazra, Anasuya 
Kundu, Pratap K. Das, Anjan Kumar Chandra from Calcutta University for providing me the 
moral support during the initial stage of this work. I also would like to thank Professor S. 
D. Adhikari, Dr. R. Thangadurai, Anamitra Mukherjee, Arijit Saha, Pomita Ghosal, Atashi 

Chatterjee, Kalpataru Pradhan, Sudhir Kumar Gupta, Subhaditya Bhattacharya, Priyotosh 
Bandyopadhyay, R. Srikanth H., Shailesh Lai, Mahender Singh for sharing light moments with 
me and for their constant encouragement during my stay at Harish-Chandra Research Institute. 

I gratefully acknowledge the cluster facilities of Harish-Chandra Research Institute for com- 
putational work. I received help from Professor Jasjeet Singh Bagla regarding the use of this 



vi 



cluster facility. 1 would like to thank Mr. Sanjai Verma, Mr. Shahid Ali Farooqui and Mr. 
Chandan Kanaojia for providing help regarding computers at various stages during this work. 



I would like to acknowledge support from the project (SP/S2/K-10/2001) of the Department of 
Science and Technology, India for providing the financial support during the initial part of this 
work. I also acknowledge support from the Xlth Plan Neutrino Project of the Harish- Chandra 
Research Institute. 

Last but not least with all my heart, my deepest thanks are devoted to my beloved family : my 
father, my mother and my brother, Rakesh, for their ever patience and constant support, their 
daily encouragements and boundless love. All these things I could do, just for them. 



Department of Physics, Sanjib Kumar Agarwalla 

University of Calcutta. September 2008. 



vii 



viii 



ABSTRACT OF THE THESIS 



Neutrino physics is a very intense field of research having implications in different branches of 
physics, such as high energy physics, quantum field theory, cosmology, astrophysics, nuclear 
physics and geophysics. Spectacular results on neutrino oscillations in the last several years 
have triggered a lot of interest in neutrinos, from experimental as well as theoretical point of 
view, and many future neutrino experiments are in preparation or under discussion to sharpen 
our understanding about these tiny particles. This thesis addresses several aspects of these 
issues. 

We have studied the physics reach of an experiment where neutrinos produced in a beta-beam 
facihty at CERN are observed in a large magnetized iron calorimeter (ICAL) at the India- 
based Neutrino Observatory (INO). The idea of beta-beam is based on the production of a 
pure, intense, collimated beam of electron neutrinos or their antiparticles via the beta decay 
of accelerated radioactive ions circulating in a storage ring. Interestingly, the CERN-INO 
distance of 7152 km happens to be tantalizingly close to the so-called "magic" baseline where 
the sensitivity to the neutrino mass ordering (sign of Amg^ = mg — ml) and more importantly, 
^13, goes up significantly, while the sensitivity to the unknown CP phase is absent. This 
permits such an experiment involving the golden channel to make precise measurements of 
the mixing angle ^13 and neutrino mass hierarchy avoiding the issues of intrinsic degeneracies 
and correlations which plague other baselines. 

We propound the possibility of using large matter effects in the survival channel. Pee, at long 
baselines for determination of the neutrino mass ordering and the mixing angle ^13. Matter 
effects in the transition probabilites Pg^ Per act in consonance to give an almost two-fold 
effect in the survival channel. In addition, the problem of spurious solutions due to the leptonic 
CP phase and the atmospheric mixing angle 6^23 does not crop up. Thus a beta-beam enables 
one to exclusively study the Pee survival probability with the help of proposed megaton class 
water detectors like UNO, HyperKamiokande, MEMPHYS. 

We have also explored the possibility of detecting new physics signals in dedicated experiments 
using a near and a far detector and a beta-beam source. We focus on the possible impact 
of flavor-changing and flavor-diagonal neutral current interactions that might crop up at the 
production point, in the oscillation stage, or at the detection point when one will deal with 
these upcoming facilities. For an example, the R-parity violating Supersymmetric model allows 
these kind of interactions and long baseline neutrino oscillation experiments may well emerge 
as test beds for this kind of models. 
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Neutrino of Love 



I go undetected 
In all my interactions 

I cannot be seen 
From any point of view 
You won't know if I'm here 
Except when I'm gone 
I'm the neutrino of love 
And I'm coming over you 

You cannot keep me in a cage 
No matter how thick the walls 

I will escape 
You cannot hold me in a box 
Cannot bind me with a lock 
Cannot keep me anyway 
I'm not afraid of the dark 
I'm the neutrino 
Neutrino of love 
I'm the neutrino 
Neutrino of love 

I go undetected 
In all my interactions 

I cannot be seen 
Prom any point of view 
You won't know if I'm here 
Except when I'm gone 
I'm the neutrino of love 
And I'm coming over you 

I'm the neutrino baby 
Neutrino of love 

Cannot inhibit my infiltration 
Neutrino 

Cannot prevent my penetration 
Neutrino 

I am the neutrino 



Dylan Casey: guiteir and vocals, 2001 
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Chapter 1 

A Preamble to Neutrino Physics 



The musical description of neutrinos by Dylan Casey in his song "Neutrino of Love" is really 
fantastic. Yes indeed, neutrinos are elusive, mysterious, yet abundant. Despite that (or be- 
cause of that!), even after fifty years of it's discovery, it still poses many mysteries and creates 
challenges to the physicists who want to detect it. Like electrons, they are elementary particles. 
F. Reines would narrate neutrino as, it is "...the most tiny quantity of reality ever imagined by 
a human being" . 

Neutrino physics is a very intense and exciting field of research having wide range of implications 
in high energy physics, quantum field theory, cosmology, astrophysics, nuclear physics, and 
geophysics. Marvellous results on neutrino oscillations in the last several years have triggered 
a lot of enthusiasm and interest in neutrinos, from experimental as well as theoretical point of 
view. One of the most important facts is that neutrino physics is a data driven field - for several 
years now, new data are pouring at an outstanding rate. Our understanding of neutrinos has 
improved dramatically in the past ten years and there is no doubt that neutrino oscillation 
is an exclusive example of experimental evidence for physics beyond the Standard Model of 
particle physics. This success sets a fantastic example of a road-map in which both theoretical 
understanding and experimental achievements have walked hand in hand to provide us with the 
first evidence of physics beyond the Standard Model. These developments culminated in the 
Nobel prize for physics in the year 2002, which was awarded to two pioneers in neutrino physics. 
Masatoshi Koshiba was awarded the prize for the detection of neutrinos from a supernova and 
Ray Davis Jr. for his detection of solar neutrinos. 

Neutrino physics is now poised to move into the precision regime. Active attempts are under 
way to commence the era of precision neutrino measurement science which will surely widen 
the horizon of our knowledge about neutrinos. A number of high-precision neutrino oscillation 
experiments have been contrived to sharpen our understanding about these tiny particles. This 
is the right time to ask how different planned/proposed next generation experiments in the 
coming decades would perform to explicate the nature of neutrinos and our thrust for new 
physics. This thesis is an effort to have a look on several aspects of these issues. 
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1.1 Neutrino In a Nutshell 



Neutrinos are electrically neutral particles of spin | with a very tiny mass, almost 500 000 times 
smaller than the mass of the electron, which itself is 2000 times smaller than the proton mass. 
There are at least three species (or flavours) of very light neutrinos, z/g, and Ur, which are 
left handed, and their antiparticles, u^, i^^i and which are right handed. After the photon, 
the neutrino is the most abundant particle in the Universe : each cubic meter of the Universe 
contains about 30 million neutrinos, which are remnant from the Big Bang, similar to the well 
known cosmic microwave background. It also arrives "unscathed" from the farthest reaches of 
the Universe, carrying information about its source. The interactions of neutrinos are mediated 
by heavy and bosons and therefore at low energies they talk feebly with ordinary matter 
and pass through the Earth very much like hght through a crystal. If a matter target as big 
as Earth is placed in front of 100 billion neutrinos, only one of them is likely to interact with 
it. The mean free path of a 1 MeV neutrino in lead is about 1 light year! Therefore neutrino 
detection requires very large detectors and/or very intense neutrino beams. 



1.2 Neutrino Odyssey 



Let us have a look at the incredible journey of discovery into one of Nature's most elusive 
particles. In a letter to colleagues on 4th December, 1930, Wolfgang Pauli [1] postulated the 
existence of neutrinos to guarantee the energy conservation in radioactive beta-decay. After 
the discovery of the neutron by James Chadwick two years later, it was first speculated that 
the particle predicted by Pauli could be the neutron. However, soon it was realized that 
Pauli's particle had to be much hghter than the neutron. In 1933, Enrico Fermi introduced 
the name neutrino, where he used the Italian syllable "-ino" to indicate "small neutron". 
More than two decades after Pauli's letter proposing the neutrino, in 1956, Clyde Cowan and 
Frederick Reines [2] observed the antineutrinos (the antimatter partners of neutrinos) emitted 
by a nuclear reactor. This neutrino is later determined to be the partner of the electron. In 1969, 
neutrinos produced by the Sun's burning were detected by Ray Davis with a detector based 
on Chlorine in an underground laboratory in the Homestake mine in USA. This experiment 
reported that less than half the expected neutrinos were detected. This originated the long- 
standing "solar neutrino problem". The scope that the missing electron neutrinos may have 
transformed into another type (undetectable to this experiment) was soon suggested, but the 
lack of our knowledge of the solar model on which the expected neutrino rates were based was 
initially considered a more likely explanation. 

In 1987, neutrinos from a supernova in the Large Magellanic Cloud were also detected. Only 
19 events were observed [3-5] and they established the standard picture of core-collapse super- 
novae. In recent years, several experiments could confirm the existence of neutrino oscillations. 
In 1998, the Super-Kamiokande experiment [6] reported the evidence for oscillations of atmo- 
spheric neutrinos. It was a crucial juncture for neutrino physics. Neutrino oscillation demands 
that neutrinos do have a mass and the observation of large mixing angle was completely beyond 
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the range of our imagination because in analogy to the quark mixing, it was the common behef 
that if neutrinos mixed at all then the mixing should be small. 

The year 2002 was a spectacular year for neutrino physics. The neutral current (NC) data of 
the SNO [7] solar neutrino experiment provided an independent determination of the total flux 
of active neutrinos from the Sun. The combined SNO and other solar neutrino data finally 
could establish an explanation of the longstanding solar neutrino problem in terms of neutrino 
oscillations. The KamLAND [8, 9] reactor neutrino experiment confirmed the oscillation hy- 
pothesis observing disappearance of z/g and constrained the mixing parameters to the so-called 
LMA-MSW [10] solution. The K2K [11] experiment is the first long baseline experiment which 
uses novel man-made accelerator beams. This experiment provided an independent confirma- 
tion of the explanation for the observed atmospheric neutrino anomaly by supplying the data 
which is consistent with it's oscillation interpretation. 

All the results from these dedicated experiments forced us to recall the pioneering work by 
Gribov and Pontecorvo [12, 13]. In 1968, they showed that flavour conversions can arise if 
neutrinos arc massive and mixed. Neutrino oscillation is the only phenomenon which can 
describe the outcome of all these experiments and it can easily explain the disappearance of 
both atmospheric i/^^'s and solar UgS. There are two recent very good reviews [14, 15] on this 
topic in general which can give us more insight. 



1.3 What are the Main Sources of Neutrinos? 



Neutrinos arc the most common matter particles in the universe. Neutrinos are produced via 
weak interactions (like beta-decays in atomic nuclei). In number, they exceed the constituents of 
ordinary matter (electrons, protons, neutrons) by a factor of ten billion. We have used neutrinos 
from many different origins with different energy ranges to study neutrino oscillations and the 
properties of neutrino sources. Mainly we can classify the neutrino sources into two categories : 



1. Natural sources of neutrinos. 

2. Artificial man-made sources of neutrinos. 



1.3.1 Natural Sources of Neutrinos 



The Sun 



One of the strongest neutrino sources is our Sun. The Sun shines not only in hght but also 
in electron neutrinos produced in the thermonuclear reactions which generate the solar energy. 
These reactions take place via two main chains, the pp chain and the CNO cycle. There are 
five reactions which produce in the pp chain and three in the CNO cycle. Both chains result 
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in the overall fusion of protons into ^He : 

4p He + 2e+ + 2i/e + 7, (1-1) 

where the energy released in the reaction, Q = Anip — nuHe — 2me — 26 MeV, is mostly 
radiated through the photons and the neutrinos carry only a small fraction of it, {E2uJ = 0.59 
MeV. Hence, the observation of solar neutrinos provides direct evidence for the nuclear process 
in the center of the sun. Moreover, current solar neutrino data allow a quantitative test of 
the Standard Solar Model (SSM) [16]. By the measurement of the solar neutrino flux the 
temperature in the center of the sun can be determined with the impressive accuracy of 1%. 

The Sun emits about 2 x 10^^ electron neutrinos per second, leading to the neutrino flux at the 
surface of the earth of ~ 6 x 10^° cm~^s~^ in the energy range E < 0.42 MeV and ~ 5 x 10^ 
cm^^s"^ in the energy range 0.8 MeV ^ < 15 MeV. The detection mechanism of solar 
neutrinos is very sophisticated and involves mainly the radiochemical processes (Homestake, 
Gallex, Sage, GNO detectors) [17, 18] and the water Cerenkov techniques (Kamiokande, Su- 
perKamiokande, SNO detectors) [6,7,19,20]. 



The Ecirth's Atmosphere 



Earth's atmosphere is another crucial source of electron and muon neutrinos and their antipar- 
ticles which are created in the hadronic showers induced by primary cosmic rays. Atmospheric 
neutrinos were first detected in the 1960 's by the Kolar Gold Field experiment in India [21] and 
the underground experiments in South Africa [22] . The following chain of reactions depicts the 
main mechanism of production of the atmospheric neutrinos : 

p{a,...) + Atr Tr^{K^) + X 

7r±(K±) ^ /x± + i/^(P^) (1.2) 

Atmospheric neutrinos can be directly detected in large mass underground detectors predomi- 
nantly by means of their charged current (CC) interactions : 

i/e(Pe)+^^e-(e+)+X, 

v^{V^)+A^li-{li+)+X. (1.3) 

Atmospheric neutrinos cover a wide range of energy starting from few MeV to hundreds of 
GeV. The typical flux of atmospheric neutrinos at the earth's surface is ~ 10~^ cm~^s~^. 
These neutrinos are observed in underground experiments with bigger and better detectors 
using different detection techniques and leading to different type of events depending on their 
energy. In the last ten years, the high precision and large statistics data from the Super- 
Kamiokande [6] experiment (using water Cerenkov detectors) has played an important role 
to solve the atmospheric neutrino puzzle through the concept of neutrino oscillation. It has 
received important confirmation from the iron calorimeter detectors Soudan2 [23] and MACRO 
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[24]. In June 1998, in the Neutrino98 conference, Super-Kamiokande collaboration presented 
evidence of oscillations based on the angular distribution for their contained event data 
sample. 

Recently, in India, to observe atmospheric neutrinos, the proposal for a large magnetized iron 
calorimeter detector (ICAL) with charge identification capability is being evaluated by the 
INO [25] collaboration. We will discuss about this detector in detail later. 



The Earth's Crust 



The Earth contains a certain amount of natural radioactivity, and the decay of these radioactive 
elements is an important and perhaps main source of gcothermal heat. These same decays also 
generate particles known as geoneutrinos. "Geoneutrinos" are electron antineutrinos produced 
by beta-decays of the unstable, radioactive nuclei in the decay chains of ^^^U, ^^^Th, "^^K. The 
amount of neutrinos coming from this natural radioactivity is huge : about 6 millions per second 
and per cm^ with the energy E ^ 1 MeV. 

A careful observation of the arrival directions of neutrinos generated in the decay of natural 
radioactive elements in the Earth's interior can give us a three-dimensional view of the Earth's 
composition and shell structure. This will provide a new and detailed understanding of the 
origin of the Earth's gcothermal heat, and will finally answer the question of how much heat 
comes from radioactive decays, and how much is "primordial" heat leftover from the birth of the 
Earth. The mapping of the Earth's interior might also help give answers to such questions as 
"What powers the magnetic field of the Earth?" and "What dominates the geo dynamo?" . To 
actually take a neutrino picture of the Earth is quite challenging technically, but not impossible. 

Recently the KamLAND experiment, which was primarily designed to measure anti-neutrinos 
from nuclear reactors, reported 9 events [26] due to geoneutrinos. This marks the first detec- 
tion of neutrinos from the Earth's interior, and already demonstrates that radioactivity is an 
important heat source for the Earth. 



The Supernovae within our galetxy 



The life of a star ends often with a huge explosion called Supernova, which can be even brighter 
than a whole galaxy. However, only a tiny amount of the total released energy is emitted as 
light. About 99% of the energy is released in the form of neutrinos having energies in the range 
10-30 MeV. Roughly, they emit ~ 6 x 10^^ neutrinos and antineutrinos of all fiavours over the 
time interval of about ten seconds. The neutrino luminosity of a gravitational collapse-driven 
supernova is typically 100 times its optical luminosity. The neutrino signal emerges from the 
core of a star promptly after core collapse, whereas the photon signal may take hours or days 
to emerge from the stellar envelope. The neutrino signal can therefore give information about 
the very early stages of core collapse, which is inaccessible to other kinds of astronomy. 
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In 1987, there was a Supernova in the Large Magellanic Cloud within our galaxy. Indeed, 
about 19 neutrinos have been observed, which confirmed our basic understanding of the Super- 
nova explosion mechanism. Neutrinos from gravitational collapse can be detected in various 
ways. For water Cerenkov detectors, such as Super-Kamiokande, the most important detection 
reaction is the absorption of electron antineutrinos on protons : 

Ue+P^n + e^. (1.4) 

The positron from this reaction, which retains most of the energy of the incoming neutrino, is 
detected from it's Cerenkov light. 



The Big-Bang 



The most extraordinary explosion of all, the Big Bang, created more neutrinos than any other 
source which has existed since. The "standard" model of the Big-Bang predicts, like for the 
photons, a cosmic background of neutrinos. These relic neutrinos still exist, have a number 
density of about 110 cm~^ for each neutrino species and a black-body spectrum with the average 
energy of about 0.0005 eV. The energy of these neutrinos is too small so that no experiment, 
even very huge, has been able to detect them. 



The Ultra-High Energetic Cosmic Neutrino Sources 



Astrophysical neutrinos can also be produced with remarkably high energies. Neutrinos born 
in the center of an active galaxy can arrive on Earth with more energy than we will ever be able 
to create with a terrestrial accelerator. Gamma-ray bursts (GRB) and active galactic nuclei 
(AGN) jets have been suggested as sources of high-energy, > 10^^ eV, neutrinos, with fluxes 
that may be detectable in a kilometer-squared effective area telescope [27]. 

1.3.2 Artificial Man-made Sources of Neutrinos 



The Nucle£tr Reactors (Power Plants) 



Electron antineutrinos with E ~ MeV arc produced copiously in the process of generating 
electrical power in nuclear power plants using controlled fission technique. A 3 GW plant 
releases about 7.7 x 10^*^ Pg per second and creates a flux of ~ 6 x 10^^ cm~'^s~^ at 100 m. Due 
to the low energy, e's are the only charged leptons which can be produced in the neutrino CC 
interaction. If the z/g oscillated to another flavour, its CC interaction would not be observed. 
Therefore only disappearance experiments can be performed with reactors. 

The KamLAND experiment [8], a 1000 ton liquid scintillation detector, is currently in operation 
in the Kamioka mine in Japan. This underground site is located at an average distance of 150- 
210 km from several Japanese nuclear power stations. This experiment has played a crucial 
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role to establish the fact that the solar neutrino puzzle can be explained by the so-called 
LMA-MSW [10] solution. Gosgen [28], Krasnoyarsk [29], Bugey [30], CHOOZ [31] and Palo 
Verde [32] are the examples of reactor experiments which are performed at relatively short or 
intermediate baselines. It is worth-while to mention here that none of these experiments find a 
positive evidence of flavour mixing. 



The Pcirticle Accelerators 



Now neutrino physics has entered into precision age from discovery era. Only in a well tuned, 
fully optimized environment will it be possible to perform precision measurements of neutrino 
oscillation parameters. A high intensity neutrino source with known spectrum is most desirable 
for precision measurements, the consensus direction for the future. Man-made accelerator based 
neutrino beams of the energy ranging typically between 30 MeV to 30 GeV are the novel, 
intense sources of neutrinos, an important tool for studying neutrino properties. Beam shape 
parameters play a very key role for the measurement of oscillation length, while the absolute 
normalization is crucial for the determination of the mixing angle. One of the most important 
features of the neutrino beam is that one can control the flux of the produced neutrinos and 
can tune the main parameters that govern the systematic uncertainties on the neutrino fluxes. 
Currently, there are three widely different schemes for producing neutrino beams and they 
mainly differ from each other on the issue of what particle is decaying (pion decay, muon decay 
and radioactive ion decay) to give rise to the neutrinos. For a detailed discussion of future 
beams and their comparison see [33, 34] and references therein. 



1. Conventional Neutrino Beams 



Conventional neutrino beams are produced by shooting a target with as many protons as 
can be provided, and then focusing the produced mesons (mostly pions and some kaons) 
into an evacuated decay pipe where they are allowed to decay. The mesons will decay 
primarily to muons and muon neutrinos. This decay chain is quite similar in nature to 
the decay process through which the atmospheric neutrinos are produced. But one should 
keep in mind that the energy spectrum of the conventional neutrino beam is quite different 
compared to that of atmospheric one. Prom meson decay (two-body decay) kinematics it 
follows that the neutrino energy is given by 

2 2 

where 7 is the Lorentz boost of the parent meson, Et^{^k) it's energy and 9 the angle of the 
neutrino with respect to the meson flight direction. The polarity of the focusing device 
has to be reversed to produce a beam of v^. In both the cases, there is always some 
contamination of v^. or due to the three-body decays of the kaons and daughter muons. 

We can classify the conventional neutrino beams into three categories : the Wide Band 
Beams (WBB), the Narrow Band Beams (NBB) and the Off- Axis Beams (OAB). 
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The main feature of WBB is that they have wide energy spectrum with high neutrino flux. 
They are perfectly suited to make discoveries. But they have some limitations. In a WBB, 
the irreducible fraction of z/g originating from the meson decays is quite significant and 
these produce electrons inside the detector which play the role of intrinsic beam related 
background when we try to extract interesting results using the Vf. appearance 

channel. Another problem is that WBB comes with a tail of high energy neutrinos and 
often these neutrinos produce 7r° inside the detector via NC process. Now, the early 
showering of gamma's from the 7r° decay can be misidentified for a z/g CC interaction. If 
the signal corresponds to a small part of the energy spectrum, it could be overwhelmed 
by the beam induced background coming from the region which is outside the signal. 

The NBB are quite opposite in nature. This facility can provide us almost monochro- 
matic energy spectra by judiciously choosing a small momentum bite of the parent vr and 
K. However, the neutrino yield is quite suppressed which causes problem for oscillation 
searches. 



The OAB [35] is a classic example of a neutrino beam with high flux and a narrow energy 
spectrum. This technique requires designing a beam-line which can produce and focus 
a wide range of mesons in a given direction (as in the WBB case), but then placing the 
detectors at an angle with respect to that direction. Since the pion decay is a two-body 
process (see Eq. II. 5p . we can obtain neutrinos of a given energy at a given angle between 
the pion direction and the detector location. But the most important feature to be noted 
that at this given angle, the energy of the neutrino produced in the decay of the pion 
becomes practically independent of the boosted pion energy. Now, if we place a detector 
at this particular angle with respect to the decay pipe then it will see a neutrino beam with 
a very narrow energy spread compared to the on-axis beam. Furthermore, this off-axis 
technology helps us to reduce the background coming from the intrinsic i^e contamination 
in the beam and a smaller fraction of high energy tails reduces the background from NC 
events. As a result, it improves the signal-to-background ratio a lot. 

But it is important to mention that, independent of the adopted solution, all conventional 
neutrino beams have some common problems. The first major drawback is that the hadron 
production in the proton-target interaction has large uncertainties due to lack of data and 
theoretical difficulties in describing hadronic processes. It creates obstacles in predicting 
the neutrino flux and spectrum with good accuracy. Secondly, in addition to the dominant 
flavour in the beam (typically z/^) there is a contamination (at the few percent level) from 
other flavours (z/^, v^. and z/g) resulting into a "multiflavour" neutrino beam. 

A recently terminated conventional beam experiment is the K2K experiment [11] where 
a neutrino beam was directed towards the Super-Kamiokande detector from the KEK 
accelerator. This experiment has already confirmed the disappearance of as predicted 
by atmospheric neutrino data. The MINOS [36] experiment in US, and the CERN to 
Gran Sasso (CNGS) experiment OPERA [37] are the two conventional beam experiments 
which are collecting data now. Another CNGS experiment ICARUS [38] will start taking 
data soon. 



2. SuperBeams 
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Figure 1.1: The schematic layout of the neutrino factory set-up. 
http://www.cap.bnLgov/mumu/project/ISS/ is the source of this figure. 



The web address 



The technology of conventional beam experiments with some technical upgrades is known 
as superbeam. These experiments are 'super' in the sense that they will use proton beams 
of unprecedented strength around 1-4 MW and detectors with large fiducial mass. All 
superbeams use a near detector for a better control of the systematics. The most advanced 
superbeam proposals are the J-PARC to Super-Kamiokande experiment (T2K) [39] in 
Japan, and the NuMI off-axis experiment (NOi/A) [40], using a neutrino beam produced 
at Fermilab in US. 



3. 



Neutrino Factory 



The term "neutrino factory" [41] has been associated to describe neutrino beams created 
by the decays of high energy muons which are circulated in a storage ring with long 
straight sections. The decay of muons in these straight sections produces an intense, well 
known and pure beam of z/^ and z/g. If /i"*" are stored, /i"*" e^VeV^ decays generate a beam 
consisting of equal numbers of Ve and v^. The overall layout is shown in Fig. 11.11 The 
muons are obtained via pion decay. To make a muon beam which can be accelerated, first 
the muons have to be cooled in phase space with the help of ionization cooling technique 
which is being studied by the MICE experiment [42]. There are other ways also being 
examined to tackle this problem. Several design studies have been performed in this 
direction in Europe, the United States and Japan [33,43]. Typical neutrino factories are 
being considered with muon energies ranging from 20 GeV to 50 GeV with ~ 10^^ useful 
muon decays per year inside the storage ring. 

In a neutrino factory experiment with stored the appearance channel, Ve give rise 
to /x"*" in the detector via CC deep inelastic scattering. These anti-muons are called wrong 
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Table 1.1: Beta decay parameters: lifetime r, electron total end-point energy E'o, /-value or 
Fermi integral and decay fraction for various ions [74]. 




■ Baseline Beta Beam 
facility comprises these 
sections 

— Proton Driver 

oSPL(»4SeV) 

— ISOL Target 

spallation neutrons or 
direct protons 

— Ion Source 

pulsed ECR 



- Acceleration 

.linoc, RCS, PS, SPS 
■ Decoy Ring 

c7000 m; 2500 m straight 



Ring 



Vg BclH lii:;im 



Baseline concept assumes 
am PS, SPS 


ECR 


SPS ^ 




PS 




Use of Tevotron also 


RCS 




being considered 







Figure 1.2: Left panel depicts the beta-decay process which is the source of pure beta- 
beam. The proposed schematic layout of the beta-beam set-up at CERN (http://beta- 
beam.web.cern.ch/beta-beam/task/index.asp) has been displayed in the right panel. 



sign muon events, since they have the opposite charge relative to the muons produced by 
the in the beam itself. These wrong sign muons have to be cleanly separated inside the 
detector from the muons created by the surviving i/^. This charge identification can be 
achieved with a magnetized iron calorimeter detector like the proposed ICAL@INO [25]. 

Muon decay is well known and it will provide the neutrino flux and spectrum with minimal 
systematic uncertainties, compared to conventional neutrino beams. Radiative effects on 
the muon decay is negligible. Ultimately, the flux from a neutrino factory is expected 
to be known with a precision of the order of 10~^. Another vital point to be noted is 
that a neutrino factory beam has a sharp cut-off at the energy of the stored muons which 
reduces the background from NC events. In a neutrino factory, the neutrinos with higher 
energy opens up the possibility to study the oscillation channels like z/^ and v^. — > 
because one can produce tau leptons inside the detector [44]. 



4. Beta-beam 
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Zucchelli [45] put forward the novel idea of a beta-beam [46-67], which is based on 
the concept of creating a pure, well understood, intense, coUimated beam of v^, or v^, 
through the beta-decay (see left panel of Fig. 11.21) of completely ionized radioactive ions. 
Firstly, radioactive nuclides are created by impinging a target by accelerated protons. 
These unstable nuchdes are collected, fully ionized, bunched, accelerated and then stored 
in a decay ring (see for e.g. [68,69]). The decay of these highly boosted ions in the 
straight sections of the decay ring [70] produces the so-called beta-beam. Feasibility of 
this proposal and its physics potential is being studied in depth [71], and will take full 
advantage of the existing accelerator complex at CERN and/or FNAL (Fermi National 
Accelerator Laboratory). The proposed schematic layout of the beta-beam set-up at 
CERN0 is given in the right panel of Fig. II. 2[ The main future challenge lies in building 
an intense proton driver and the hippodrome-shaped decay ring which are essential for 
this programme. 

It has been proposed to produce v^. beams through the decay of highly accelerated ^^Ne 
ions (i^Ne + e+ + z/e) and from "^He (^He ^Li + e" + z/e) [45,70]. More recently, 

(^B — > |Be + e+ + z/e) and ^Li (^Li — > |Be + e" + z7e) [72,73] with much larger end-point 
energy have been suggested as alternate sources since these ions can yield higher energy z/g 
and z/g respectively, with lower values of the Lorentz boost 7 [49,50,52-55,61,66]. Details 
of the four beta-beam candidate ions can be found in Table [LTl ®B and ®Li decay to the 
broad 2+ first excited state of ^Be which has an energy width of ~ 1.5 MeV. It may be 
possible to store radioactive ions producing beams with both polarities in the same ring. 
This will enable running the experiment in the z/g and Vf. modes simultaneously. 

In the low 7 design of beta-beams, the standard luminosity taken for the ^'^Ne and ^He are 
1.1 X 10^^ (z^e) and 2.9 x 10^® (z/g) useful decays per year, respectively [75]. We will discuss 
about this newly proposed facility later in great detail. 



1.4 Neutrino : "NU" Horizons 



We live in an exciting time when the light of new discoveries is breaking apart our long-held 
picture of the Standard Model. This revolution began in part with the widely confirmed 
assertion that neutrinos have mass, and it will continue to be waged by upcoming neutrino 
experiments. Spectacular results from a series of experiments over the last four decades [6, 
8,11,17-19,31,36,76-79] have firmly established the phenomenon of neutrino oscillation and 
paved the way for the "golden" age of neutrino physics. Since neutrino oscillations can occur 
only if there is a mass difference between at least two neutrinos, an observation of this effect 
proves that at least one non-zero neutrino mass exists. 

Neutrinos are strictly massless in the Standard Model of particle physics and the finite neutrino 
masses required by the experimental data provide the first hint for physics beyond the Standard 
Model, and make an extension of the theory necessary. No doubt that this has put the Standard 



detailed R&D on this issue is being pursued by Mats Lindroos and his coUaborators. 
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Model in a paradoxical situation. Moreover, the fact that neutrino masses are so tiny (very 
much smaller than that of any other known fermion) should find an explanation in the new 
theory. 

Recent discoveries on neutrinos might provide unique information on a more complete theory 
of elementary particles. The sensitivity of neutrino experiments to very tiny mass scales might 
provide the scope to learn something about physics at very high energy scales (i.e., at very 
small distances), which will never be accessible in particle accelerator experiments. Therefore, 
information from neutrinos is complementary to the one from accelerator experiments, and it 
may provide a key to a so-called Grand Unified Theory, in which the electromagnetic, the weak 
and the strong interactions are unified to one fundamental force. 

Another puzzle of modern physics is the origin of matter. In the so-called Leptogenesis mecha- 
nism the origin of matter in the very first moments after the Big Bang is related to neutrinos. 
In that theory the small asymmetry between matter and anti-matter is generated by processes 
involving neutrinos in the early stage of the Universe. In this way a theory of neutrino may 
even provide the reason for our existence. Neutrinos have played a key role in shaping the 
Universe as we see today. We have just started our journey in the mysterious world of neutri- 
nos, a tiny creature of Nature. A long journey is waiting for us ahead and many experimental 
approaches are required to get the full view. In the near future, the Large Hadron Collider 
(LHC) will start it's quest for Higgs and it is expected that the LHC will explore the mechanism 
of electroweak symmetry breaking and provide clues of new heavy degrees of freedom. This 
will certainly boost up the future road map of the neutrino physics programme and it is for 
sure that neutrino physics, a bit player on the physics stage in yesteryears, has now donned a 
central role and will play a crucial part in the high energy physics programme. 



1.5 Layout of the Doctoral Work 



We organize the description of the doctoral work in the following way. The first part of chapter 
2 deals with the basic introduction to the qiiantum mechanics of neutrino oscillation in vacuum 
under both two and three flavour frameworks. Then we discuss the importance of matter 
effects in neutrino oscillations. In the second half of chapter 2 we take a look at our present 
understanding of neutrino parameters and we identify the major unknowns in the neutrino 
sector. Finally we close chapter 2 by giving a brief note on the future neutrino road-map based 
on long baseline experiments. In chapter 3 we underscore in detail the physics advantage of 
an experimental set-up where neutrinos produced in a beta-beam facility at CERN would be 
observed in the proposed large magnetized iron calorimeter detector (ICAL) at the India-based 
Neutrino Observatory (INO). The CERN-INO beta-beam set-up offers an excellent avenue to 
use the "Golden" channel {i/e — > J^n) oscillation probability for a simultaneous determination of 
the neutrino mass ordering and ^13. The merit of the earth matter effects in the z/g z/g survival 
probability at long baselines in order to cleanly determine the third leptonic mixing angle ^13 
and the sign of the atmospheric neutrino mass squared difference, Aml^, using a beta-beam as 
a source has been discussed in chapter 4. In chapter 5 we study the possibility of detecting 
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new physics signals in a dedicated neutrino beta-beam experiment with the source at CERN 
and the detector at the proposed INO. These new physics signals arise in the R-parity violating 
supersymmetric models (RPVSM) [80] due to the flavor-changing neutral current (FCNC) and 
flavor-diagonal neutral current (FDNC) interactions of neutrinos with matter in the oscillation 
stage. In chapter 6 we show that a detector placed near a beta-beam storage ring can probe 
lepton number violating interactions that might crop up at the production and the detection 
point, as predicted by supersymmetric theories with R-parity non- conservation. Finally we end 
with the conclusion and an outlook in chapter 7. 
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Chapter 2 

Neutrino Oscillations Revisited 



This chapter is organized as follows. In the first section we discuss the role of neutrinos in the 
Standard Model. An introduction to the basic formalism of neutrino oscillations in vacuum has 
been given in the subsequent section. We consider both two and three flavour frameworks. In 
the following section we concentrate on the importance of matter effects in neutrino oscillations. 
Then we turn our focus on the present global understanding of the neutrino mass-mixing 
parameters. This is to set the stage for the next section where we take a drive to the unknown 
territories of the neutrino sector. Finally we take a look at the future neutrino road-map. 



2.1 Neutrinos in the Standard Model 



The Standard Model is a quantum field theory and it unifies strong, weak and electromagnetic 
interactions. The Standard Model is based on the gauge group SU{3)c x SU{2)l x U{1)y, 
broken spontaneously by the Higgs mechanism to SU{?))c x U{l)em- The transformation groups 
SU{?>)c, SU{2)l and U{1)y correspond to quantum numbers called colour, weak isospin and 
weak hypercharge {Y) respectively. Experimental observations suggest that the matter fields 
are fermionic, consisting of leptons and quarks. In the Standard Model, these leptons and 
quarks are chiral fermions : 

i^R,L = \il±l5)^, (2.1) 

where ip = ipi + '^R- The left-handed components transform as doublets under SU{2)l, whereas 
the right-handed components are SU{2)l singlets. Within the framework of the Standard 
Model there are no right-handed neutrinos or left-handed antineutrinos. The particle content 
of the Standard Model is depicted in Table 12.11 It can be readily seen from Table 12.11 that 
neutrinos do not take part in strong and electromagnetic interactions and they only undergo 
weak interactions, that is, they are singlets of SU{3)c x U{l)em- Precision data of the Z-decay 
width at the e~^e~ collider at LEP indicate that there are three neutrinos which take part in 
weak interactions. These are called the active neutrinos. Often a scenario with four neutrinos 
is discussed. So this fourth neutrino, if it exists, should be sterile to weak interactions. This 



15 



Partirlps 


TVJntation 










Leptons 




1 
1 


2 
1 


-1/2 
-1 




-1 
-1 


Quarks 


'^R 


3 

3 
3 


2 

1 
1 


1/6 

2/3 
-1/3 


2/3 
-1/3 

2/3 
-1/3 



Table 2.1: Matter contents of the Standard Model with their corresponding gauge quantum 
numbers. / = 1, 2, 3 is the generation index. The electromagnetic charge listed in the last 
column is defined SlS Qem = + Y ■ 

sterile neutrino does not have any Standard Model interactions and passes undetected through 
the experimental set up. 

In weak interactions, parity is maximally violated and only the left-handed states do couple 
with the gauge bosons while the right-handed states do not, unlike in the case of electro- 
magnetism where both the states couple to the photon with the same strength. The unbroken 
electromagnetic gauge symmetry ensures the masslessness of photon. But no such fundamental 
principle predicts the neutrino to be massless. In the mass-term, the left- and right-handed 
chiral states are coupled to each other in the following way : 

m^i) = m{4)L^pji + ijRipL), (2.2) 

which can be easily obtained with the help of Eq. 12. II and using the properties of 75. Therefore 
the existence of both left-handed and right-handed chiral components is mandatory to have a 
nonzero mass of a fermion. But within the Standard Model there is no right-handed neutrino 
and as a consequence the neutrino does not have any mass. So non-zero neutrino mass, as 
required by the neutrino oscillation data, indicates the existence of physics beyond the Standard 
Model. 

As we have already mentioned, the Standard Model has three active neutrinos accompanying 
the charged lepton mass eigenstates, e, /x and r. Therefore there is a possibility of weak CC 
interactions between the neutrinos and their corresponding charged leptons which is given by : 

- = 7f E ^LiY^liK + h-c. (2.3) 
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osc 



Baseline 

Figure 2.1: The two flavour z/g to oscillation probability (Eq. 12.141) . The mixing angle 6 
governs the oscillation amplitude and the frequency depends on Am^2 = ""^i ~ ""^i^ where mi 
and 7712 are the masses of the mass eigenstates ui and i>2 respectively. 



The Standard Model neutrinos also have NC interactions : 

- ^Nc = Yl ^LiY^LiZ^. (2.4) 

2 cos Ow J 

Thus, within the Standard Model, all the neutrino interactions are described by Eqs. 12.31 and 



2.2 Oscillations in Vacuum 



Neutrino oscillation is a simple quantum mechanical phenomenon in which neutrino changes 
flavour as it propagates. In 1957, Pontecorvo [81] gave the concept of neutrino oscillation based 
on a two-level quantum system. This phenomenon arises if neutrinos have masses and there is 
mixing. 

In the course of time evolution of a stationary state |\I'/c) (which is one of the eigenstates of the 
hamiltonian of a system with the energy eigenvalue Ek), it will remain in the same state and 
will be associated with a phase in the form {h = c = 1) : 

\^dt)) = e-''''^'\^km. (2.5) 

Now if we deal with an arbitrary non- stationary state of the system which is a superposition 
of say two eigenstates (with eigenenergies Ei and E2) of the hamiltonian of the system, then it 
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will not be in the same state after time evolution. The chances that the system will remain in 
it's initial state will be an oscillatory function of time with frequency {E2 — Ei). 

In order to implement this mechanism in case of neutrinos, we consider the fact that neutrinos 
(fe, J^/^, t^t) are produced or detected via weak interactions and therefore they are referred to 
as weak-eigenstate neutrinos (denoted as z/q,) that means they are the weak doublet-partners 
of e~,fi~,T~ respectively. In such a scenario, if we assume that neutrinos are massive, then in 
general, it is not mandatory that the mass-matrix of neutrinos written in this weak (flavour) 
basis will have to be diagonal. So, it follows that the mass eigenstate neutrinos Uiji = 1,2,3 
(the basis in which the neutrino mass matrix is diagonal) are not identical to the weak or flavour 
basi^ and we have for n number of light neutrino species 

n 

Wa) = J2u:M), (2.6) 

1=1 

where U is the unitary leptonic mixing matrix known as the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) matrix [82]. This matrix is analogous to the CKM matrix in the quark sector. Because 
of Eq. 12.61 the probability of finding a neutrino created in a given flavour state to be in the 
same state (or any other flavour state) oscillates with time. Suppose a given source is producing 
a neutrino flux of given flavour |z/q) at t = x = 0, the neutrino state at a later time t is then 

n n 

Wait)) = Y.u:Mit)) = Y.^:^e-''''>^m, (2.7) 

i=l 1=1 

and propagates as an admixture of the mass eigenstates. The neutrino oscillation probability, 
i.e. the transformation probability of a flavour eigenstate neutrino {I'a) into another one lujj) is 
then 

n n 

Pap = KMMm' = (2.8) 

i=i j=i 

To obtain the oscillation probability for antineutrinos one has to replace U hj U*. For ultra- 
relativistic neutrinos with small mass one can assume pi c::^ pj = p c::^ E and we have 

E. = ^/^^+^,^P+^, (2.9) 
where E^ and rrii are, respectively, the energy and the mass of the neutrino mass eigenstate Ui. 



Now using the orthogonality relation = Sij and with the help of Eq. 12.91 we obtain the 

transition probability that an initial of energy E gets converted to a z/^ after traveling a 
distance L{= t) in vacuum as 



Pap = 5o.f3 - 4^Re[f/:,f/;^.f/;3if/ai]sin2x,, 

i<j 
n 

+ 2^Im[f/:,f/;^.f/;3if/„,]sin2X,,-, (2.10) 



^T\ie charged Icpton mass-matrix is diagonal in this basis. 
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where 

(mf -m3)L Amf- LIE , , 

Xi. = — ^ = 1.27 ^ — / . 2.11 

Am^j = mf — is known as the mass sphtting and neutrino oscillations are only sensitive 
to this mass squared difference but not to the absolute neutrino mass scale. The transition 
probability (depicted by Eq. I2.10p has an oscillatory behaviour with oscillation lengths 

n.. 47rE E(MeV) , E (GeV) 

L°f = -—y- - 2.48 m ^ \ ^ ; = 2.48 km ^ \ ^ ' 2.12 
Am\ Amf.{eV^) Am^. (eV^) ^ ^ 

and the amplitudes are proportional to the elements in the mixing matrix. Since neutrino 
oscillations can occur only if there is a mass difference between at least two neutrinos, an 
observation of this effect proves that at least one non-zero neutrino mass exists. 

In general, to construct a unitary n x n matrix we need n{n — l)/2 angles and n{n + l)/2 
phases. If neutrinos are Dirac type in nature then 2n — 1 phases can be absorbed by a proper 
re-phasing of the left-handed fields, leaving {n — l){n — 2)/2 physical phases. Therefore CP 
violation is only possible in the case of n > 3 generations. In the Majorana case there is no 
freedom to re-phase the neutrino fields; only n phases can be removed using the charged lepton 
fields, leaving n{n — l)/2 physical phases. Out of these (n — 1)(?t, — 2)/2 are the usual Dirac 
phases, while n — 1 are specific to the Majorana case and are called Majorana phases. The 
latter do not lead to any observable effects in neutrino oscillations and are not discussed further 
in this thesis. 



2.2.1 Two-Flavour case 

In a simple case with only two families of neutrinos, the mixing matrix depends on a single 
parameter 6 (known as the mixing angle), 

U=( ^"^M (2.13) 

\^ — sm cos (J J 

and there is a single mass-squared difference Am^ . Using this form of U in Eq. 12.101 one 
obtains 

P,^^,, = sm^2e sin2(1.27Am2--), (2.14) 

E 

Pu^^u^ = 1 - sin^ 29 sin'(1.27Am24), (2.15) 

E 

where Am^ is in eV^, L is in m (km) and E in MeV (GeV). In the above equations, 9 and Am^ 
(these are fundamental constants like the electron mass or the Cabibbo-angle) determine the 
oscillation amplitude and frequency respectively (see figure [2711) . The various neutrino sources 
like the Sun, the atmosphere, reactors and accelerators (discussed in the previous chapter) 
provide neutrinos with average energy varying in a wide range and the distances between the 
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sources and the detectors are also quite different which offer the possibihty to probe various 
oscillation frequencies. In order to penetrate a given value of Am^, the experiment has to be 
performed with E/L k, Am? which is equivalent to L ~ L°'^'^. 

In an oscillation experiment, if we search for a new flavour then we call it an appearance 
experiment. On the other hand if we look for a reduction in the neutrino flux from the source 
then it is known as a disappearance experiment. In both the cases, the measurement of Am^ 
and 9 are correlated with each other. Therefore an error in measuring one of the parameters 
causes an additional uncertainty on the other one. To locate the position of the peak clearly, 
we need an experiment with good enough energy resolution [AE -C LA-m?) otherwise the 
information on Am^ gets blurred and the experimental signal becomes energy independent 
(~ |sin^26'). The proper energy calibration of the detector (it provides the absolute energy 
scale) is also a very crucial issue while measuring the mass splitting. We can probe very small 
values of 9 with an appearance experiment because the measurement is done relative to zero. On 
the contrary, a disappearance experiment measures relative to unity. Therefore the impact of 
certain systematical errors is quite different in both the cases. In an appearance experiment, the 
signal is proportional to 9 and therefore an accurate understanding of the level of background 
is very vital while measuring small values of 9. But in case of a disappearance measurement 
the issue of total normalization is crucial because a large normalization error makes it tough to 
detect deviations from unity. 



2.2.2 Three-Flavour case 



Let us now consider the case of three neutrino flavours. In a three neutrino framework, the 3x3 
unitary mixing matrix U depends on three mixing angles 6'i2, ^13 and ^23 and one CP- violating 
phase 5c p (ignoring Majorana phases). The mixing matrix U can be parameterized as 



t/ = ^23^^13^12, (2.16) 

where 

/ C12 S12 0\ / ci3 si3e-^^^^\ /I 

V12 = -S12 C12 , 1^13 = 1 , \^23 = C23 S23 

V 01/ \~sne'^^^ ci3 / \0 -S23 C23 

(2.17) 

Cij = cos 9ij and Sij = sin 9ij. The neutrino mixing matrix takes the form 

C12C13 S12C13 sise-'^'^p 

[/ = ( -S12C23 - Ci2S23Sl3e*''^'P C12C23 - Si2S23Sl3e*'^'='-f' S23C13 |- (2.18) 

S12S23 - Cl2C23Si3e"^c;-P -C12S23 - Si2C23Si3e*'^C-P C23C13 



The probabilities of oscillations between various flavour states are given by Eq. 12.101 For 
three neutrino species, one can identify Am^ = Am^i > as the neutrino mass squared 
difference responsible for the solar neutrino oscillations and the dominant atmospheric neutrino 
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oscillations are caused by the |Am^| = |Am|^| = lArriggl ^ Aml^. 6*12 = 9q and 6*23 = 6'a 
are the solar and atmospheric neutrino mixing angles, respectively. The angle ^13 is the so- 
called "CHOOZ mixing angle" which connects the solar sector with the atmospheric one and 
determines the impact of the three flavour effects. Three flavour effects leave their imprints in 
an experiment which is sensitive to both the mass splittings. It is true that for the precision 
measurements of the mass squared differences or the two large mixing angles, these effects are 
not that vital but they can play a key role in the determination of the small mixing angle and 
in many cases, they decide the fate of an experiment. Like in the quark sector, in the case 
of three generations there can be CP-violation in the neutral lepton sector provided that ^13 
and the CP phase 6cp are non-zero. If CP is not conserved, the oscillation probabilities for 
neutrinos are different from those for antineutrinos. The CP-odd asymmetries are defined as 

AP„^ = P(z/„ ^ up; L) - P(z7„ ^ z7^; L) . (2.19) 
CPT invariance ensures that APap = — AP^^. With the help of Eq. 12.181 we have 



APe^ = AP^, = AP.e = 4JcP X 

where ^ 

Jcp = -cos6'i3 sin26'i3 sin 26^23 sin26'i2 sin ^c-p (2-21) 
8 

which is known as Jarlskog CP-odd invariant [83] . One can easily draw the following conclusions 
by observing Eq. 12.201 

1. The CP-odd asymmetry vanishes if 6c p is zero or 180° and it comes with full strength if 
6c P equals to 90° or 270°. 

2. It vanishes if any of the mixing angles 612, O13 or 6'23 is zero or 90°. This indicates that 
we need at least three generations to observe CP-violation which is suppressed by the 
smallest of the mixing angles, ^13. 

3. The mass squared differences satisfy the relation Am^i + Am^g + Am^g = 0. Therefore, if 
any one of the Amfj is zero then the CP-odd asymmetry vanishes immediately. It ensures 
the fact that the three flavour effects are essential for probing the leptonic CP-violation. 

There is no doubt that it would be of great importance if we can establish the CP- 
nonconservation in the leptonic sector and future long baseline neutrino oscillation experiments 
will provide a crucial hint along this direction. 
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Figure 2.2: Feynman diagrams showing neutrino scattering inside the matter. Left panel depicts 
CC interactions whereas right panel describes NC processes. 



Medium 


Matter Density 


Vcc 


Solar core 


~ 100 g/cm^ 


~ 10-1^ eV 


Earth core 


~ 10 g/cm^ 


~ 10-13 gv 


Supernova 


~ 10^^ g/cm^ 


~ eV 



Table 2.2: The impact of matter potential in various mediums. 



2.3 Oscillations in Matter 



Neutrino propagation through matter can modify oscillations significantly. When neutrinos 
travel through matter {e.g. in the Sun, Earth or a supernova), the weak interaction couples 
the neutrinos to matter and besides few hard scattering events there is also coherent forward 
elastic scattering of neutrinos with matter particles they encounter along the way. This motion 
of neutrinos can be compared with the visible light traveling through glass. The crucial point 
to be noted here is that the coherent forward elastic scattering amplitudes are not the same 
for all neutrino flavours. The ordinary matter consists of electrons, protons and neutrons but 
it does not contain any muons or tau-leptons. One can readily see from the right panel of Fig. 
12.21 that neutrinos of all three flavours (z/g, and Vt) interact with the electrons, protons and 
neutrons of matter through flavour independent NC interaction mediated by Z*^ bosons and 
these contributions are same for neutrinos of all three flavours and therefore these interactions do 
not have any impact on neutrino oscillation probabilities. Interestingly, the electron neutrinos 
have an additional contribution due to their CC interactions (see left panel of Fig. 12. 2p with 
the ambient electrons of the medium which are mediated by the exchange. This extra 
matter potential comes in the form 



A = ±2V2GFNeE , (2.22) 

where Gj? is the Fermi coupling constant, A'^g is the electron number density inside the Earth 
and E is the neutrino energy. The + sign refers to neutrinos while the — to antineutrinos. The 
electron density (A^g) is connected to the matter density (p) in the following way 

Vcc = V2G,N, ^ 7.6Fe— -^eV , (2.23) 
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where = j^^^ is the relative number density. Np, are the proton and neutron densities in 
Earth matter respectively. In an electrically neutral, isoscalar medium, we have N^. = Np = 
and Ye comes out to be 0.5. The strength of Vcc inside the Sun, Earth and a supernova is given 
in Table [2^ If we compare Vcc with /S.m? /2E then we can evaluate the importance of matter 
effects on neutrino oscillations. If we consider a neutrino of 5 GeV passing through the core of 
the Earth then Vcc is comparable with ISw? jlE (= 2.5 x 10^^^ eV if Am^ = 2.5 x 10~^ eV^). 

Staying in the two-flavour formalism, the time evolution of the flavour eigenstates in matter is 
given by the following Schrodinger equation 



where f/ is the mixing matrix defined by Eq. 12.131 In the constant matter density case the 
problem boils down to a stationary one and a trivial diagonalization of the Hamiltonian is 
required to obtain the solution. In the presence of matter, the vacuum oscillation parameters 
are mapped to the new parameter^ in the following way 



(Am^)'" = v^(Am2 cos 2^ - Af + (Am2 sin 2^)^, 
sin 2^"^ = sin2e AmV(Am2)'^. 

The so-called MSW-resonance [10] condition is met at 

Am^ cos 26* = A. 



(2.25) 



(2.26) 



At MSW-resonance, sin 2^™ = 1 (from Eq. 12.251 and I2.26P which immediately signifies that 
independent of the value of the vacuum mixing angle ^, the mixing in matter is maximal i.e. 
6"^ = 7r/4. This resonance occurs for neutrinos (antineutrinos) if Am^ is positive (negetive). 
So the sign of Am^ determines the oscillation probability. From Eq. I2.26[ the resonance energy 
can be expressed as 
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11.16GeV 
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'cos 20" 




2.8 g/cm^ 


2.5 X 10-W^ 




0.95 




P 



(2.27) 



In course of neutrino propagation through the upper Earth mantle with p = 2.8 g/cm^, the 
resonance occurs at roughly 11.2 GeV provided that [Am^l = 2.5 x 10~^eV^ and cos 26* = 0.95. 
The matter effect arises due to matter and not anti-matter and this fact is responsible for the 
observed asymmetry between neutrino and antineutrino oscillation probabilities even in the 
two neutrino case. 

In the presence of three neutrinos the time evolution of flavour eigenstates in matter can be 
written as 
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(2.28) 



^The new parameters in matter carry a superscript m. 
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Parameter 


Best fit 


3(T (1 d.o.f) 


A 2 n n— 5 A^l 

Am2i [lU eV J 
|Am|i| [10-3 eV^] 
sin^ 9i2 
sin^ 023 
sin^ 6*13 


'~7 n 

7.0 
2.4 
0.32 
0.50 
0.007 


7.1-8.3 
2.0-2.8 
0.26-0.40 
0.34-0.67 
< 0.050 



Table 2.3: Best-fit values and 3o"(l d.o.f) constraints on the oscillation parameters under three- 
fiavour scheme from global data including solar, atmospheric, reactor (KamLAND and CHOOZ) 
and accelerator (K2K and MINOS) experiments. This table has been taken from [79]. 



where f/ is the mixing matrix defined by Eq. 12.181 Even in the constant matter density 
case, it is very tough to obtain the analytical solution of the above equation. Under certain 
approximations, it is possible to obtain some expressions which are quite useful in explaining 
the underlying physics. In the next chapter, we will discuss a few such cases in detail. In 
the case of three fiavours, besides the genuine CP asymmetry caused by the CP phase ^cp-, 
we also have fake CP asymmetry induced by matter which causes obstacles in extracting the 
information on (J^p- Now let us see what is our current understanding of these various neutrino 
mass-mixing parameters? 



2.4 Present Status 



Neutrino physics has entered the precision era, with the thrust now shifting to detailed under- 
standing of the structure of the neutrino mass matrix, accurate reconstruction of which would 
shed light on the underlying new physics that gives rise to neutrino mass and mixing. The 
full mass matrix is given in terms of nine parameters, the three neutrino masses, the three 
mixing angles and the three CP violating phases. Neutrino oscillation experiments are sensi- 
tive to only two mass squared differences, all the three mixing angles and the so-called Dirac 
CP phase. The remaining parameters, comprising of the absolute neutrino mass scale and the 
two so-called Majorana phases, have to be determined elsewhere. We already have very good 
knowledge on the two mass squared differences and two of the three mixing angles. Results 
from solar neutrino experiments [17, 19] which have been collecting data for more than four 
decades have now culminated in choosing the Large Mixing Angle (LMA) solution. The latest 
addition to this huge repertoire of experimental data is the result from the on-going Borexino 
experiment [18], and this result is consistent with the LMA solution. This conclusion from 
solar neutrino experiments has been corroborated independently by the KamLAND reactor 
antineutrino experiment [8,9], and a combined analysis of the solar and KamLAND data gives 
as best-fit Am^^ = 7.6 x 10-^ eV^ and sin^ = 0.32 [9,78,79]. The other mass squared 
difference Amg^ and mixing angle ^23 are now pretty well determined by the zenith angle de- 
pendent atmospheric i/^ data in SuperKamiokande [6,19,84] and the long baseline experiments 
K2K [11,77] and MINOS [36]. The combined data from the atmospheric and long baseline 
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experiments have pinned down |Am|]^| = 2.4 x 10~^ eV^ and sin^2^23 = 1- Our knowledge on 
the third mixing angle 6*13 is restricted to an upper bound of sin^6'i3 < 0.033 from the global 
analysis of all solar, atmospheric, long baseline and reactor data, including the CHOOZ [31,76] 
results in particular. We do not have any information on 5c p and it is fully unconstrained 
yet. All the best-fit values and 3(j(l d.o.f) constraints on the oscillation parameters under a 
three-flavour scheme from global data including solar, atmospheric, reactor (KamLAND and 
CHOOZ) and accelerator (K2K and MINOS) experiments have been summarized in Table [2^31 
Recent cosmological observations (WMAP 5- year and other data) predict that the sum of three 
neutrino masses would be less than 0.19 eV [85]. In the next section, we will focus on the major 
unknowns in the neutrino sector that we would like to resolve in the next ten to fifteen years 
down the line. 



2.5 Missing Links 



Despite the spectacular achievements in the last ten years or so, a lot of information is still 
required to complete our understanding of the neutrino sector. In the following, we have to 
tried to list them up : 

• Are the neutrinos Dirac or Majorana particles? 

• How many neutrino species are there? Do sterile neutrinos exist? Are three- flavour 
oscillations enough? 

• What is the mass scale of the neutrinos? Why are neutrino masses so small? 

• Does the neutrino have a non-zero magnetic moment? 

• Why is the pattern of the neutrino mixing so different from that of the quarks? Is there 
any connection between quarks and leptons? 

• What are the precise values of ^Tn\^ and sin^ ^12? 

• What is the precise value of |Am|]^|? What is the sign of Am|]^ or the character of the 
neutrino mass hierarchy? 

• Is the ordering of the neutrino mass states hierarchical or quasi-degenerate? 

• Is sin^ 26'23 exactly maximal (= 1)? If sin^ 2^23 7^ 1, what is its octant? 

• How tiny is sin^ 2^13? Is it zero? If it is non-zero then what is the precise value of ^13? 

• Does the behaviour of neutrinos violate CP? Can we probe the existence of the three CP 
odd phases? 

• Is there any link between the low energy CP violation in the neutral lepton sector and 
the observed matter anti-matter asymmetry (baryon asymmetry) in the universe ? 
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Figure 2.3: The sign of Am|^ = 
normal or inverted hierarchical. 



is not known. The neutrino mass spectrum can be 



• Can we probe the sub dominant effects due to possible new physics in the next generation 
high precision neutrino oscillation experiments? 

• What is the importance of neutrino physics in our understanding of dark matter and dark 
energy? 

• Is there any connection between neutrino mass and leptogenesis and galaxy-cluster for- 
mation? 

• What is the role of neutrinos in connecting the predictions at the grand-unification scale 
with low energy phenomena in the framework of see-saw mechanism and supersymmetric 
extensions of the Standard Model? 



Out of these several unanswered questions (listed above), there are several issues which can 
be answered in the future dedicated long baseline neutrino oscillation experiments. Let us 
discuss a few of them. While the solar neutrino data have confirmed that Am^i > at a 
CL. of more than 6a, we still do not know what is the sign of Am'^i, i.e., the neutrino mass 
hierarchj|§. Knowing the ordering of the neutrino masses is of prime importance, because it 
dictates the structure of the neutrino mass matrix, and hence could give vital clues towards the 
underlying theory of neutrino masses and mixing. Knowing the sgn{Am1i) could have other 

^Though we call this the neutrino mass hierarchy, what we mean is basically the ordering of the neutrino mass 
states (see Fig. 12. 3p . Therefore, our discussions are valid for both hierarchical as well as quasi-degenerate mass 
spectra. We define Ato|^ = mf — rrij and refer to sgn(Arn^^) as the neutrino mass hierarchy - sgn(Am|j^) > 
is called "normal hierarchy" (NH) while sgn{Am'^i) < is called "inverted hierarchy" (IH). 
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far-reaching phenomenological consequences. For instance, if it turns out the Am|^ < and 
yet neutrino-less double beta-decay is not observed even in the very far future experiments, 
that would be a strong hint that the neutrinos are not Majorana particles (see for e.g. [86] and 
references therein). A vital task to be undertaken in the immediate future is the determination 
of the hitherto unknown mixing angle ^13. Discovery of a non-zero value for ^13 would open up 
the possibility of observing CP-violation in the lepton sector. Non-zero ^13 also brings in the 
possibility of large Earth matter effects [10, 87] for GeV energy accelerator neutrinos traveling 
over long distances. Effect of matter on neutrino oscillations depends on the sgn{Am\^ and 
is opposite for neutrinos and antineutrinos. For a given sgn(Am3]^) it enhances the oscillation 
probability in one of the channels and suppresses it in the other. Therefore, comparing the 
neutrino signal against the antineutrino signal in very long baseline experiments gives us a 
powerful tool to determine sgii{Am\{) . A way to judge the merit of an experiment is by 
observing it's limiting reach in the small mixing angle ^13 where it is still sensitive to the mass 
ordering and the leptonic CP-violation. 



2.6 Future Neutrino Road-Map 



Determining the unknown mixing angle ^13, the CP phase Sep., and the neutrino mass ordering 
have emerged as the next frontiers in this field. All these three quantities can be probed 
by experimentally measuring the so-called "golden" channel [88] oscillation probability Pg^i 
(z/g v^) or its T-conjugate channel P^e- A series of ambitious projects are under discussion 
which plan to use this oscillation channel. The on-going and near future experiments include the 
MINOS experiment in the US [36], and the CNGS experiments ICARUS [38] and OPERA [37] 
in Europe. Next experiments in line will be T2K in Japan [39] and NOz/A in US [40]. All 
these experiments will use muon neutrino beams from conventional accelerator sources in order 
to observe P^e- Collectively and in combination with short-baseline reactor experiments, such 
as Double Chooz [89], these experiments are expected to improve the bound on 61^ to about 
sin^ 2^13 < 0.01 (90% CL) [90]. The mass hierarchy and CP violation, though in principle 
accessible using the combined data from the T2K and NOi/A experiments, can be determined 
only for valueqjof sin^ 26'i3(true) close to the current bound and for some fraction of the possible 
values of the CP phase (5c'p(true). The sensitivity of these experiments is mainly restricted 
by statistics, while for larger luminosity set-ups, the intrinsic z/g background poses a natural 
limitation for experiments sensitive to z/^ oscillations into z/g. Therefore, if Nature has not been 
very kind we will need larger sophisticated experiments to complete our understanding of the 
neutrinos, possibly using an alternate smart technology with powerful, well measured neutrino 
beams from upgraded accelerator facilities and improved detector technologies. The beta-beam 
and neutrino factory can play a crucial role in this direction in future. A summary of the 
potential of selected neutrino factory and beta-beam set-ups have been compiled by the physics 
working group of the International Scoping Study for a future neutrino Factory, superbeam 



■^We distinguish between the "true" values of the oscillation parameters, which are the values chosen by 
Nature, and their fitted values. Throughout this work we denote the true value of a parameter by putting 
"(true)" after the symbol for the parameter. 
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and beta-beam, in their report [71]. In the next chapter, we will expound in detail the physics 
reach of an experimental set-up in which the proposed large magnetized iron detector at the 
India-based Neutrino Observatory (INO) [25] would serve as the far detector for a CERN-based 
monoflavour beta-beam [45]. 
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Chapter 3 



Neutrino mass ordering and ^13 with a 
magical Beta-beam experiment at INO 



3.1 Introduction 



The most promising avenue for the determination of ^13, the CP phase 5cp and sgn(Am3^) is 
the Ve ^ oscillation channel (or its T-conjugate P^e)i often referred to in the literature 
as the "golden channel" [88]. The P^^ channel can be studied in experiments which use an 
initial Ue (or Ve) beam and a detector which can efficiently see muon^j. Also in order to probe 
small values of sin^2^i3, one needs the following : 



Low backgrounds 



As far as the low background requirement is concerned, it is an advantage to use a pure 
flavor neutrino beam without any intrinsic beam contamination. A beta-beam [45] is 
well suited for this purpose. An alternative approach is the so-called neutrino factory 
[41]. The detector backgrounds coming mainly from neutral current interactions and mis- 
identification of particles, can be reduced by imposing intelligent cuts. The atmospheric 
neutrino backgrounds, which can be important for beta-beams at lower energies, can be 
suppressed using timing and directional information. 



Large statistics 



Statistics can be increased by a higher beam power and the size and efficiency of the 
detector. 



3. Reduced systematical uncertainties 



^This is complementary to the standard accelerator beams where the initial (anti)neutrino flux consists of 
(or t'p) and where the relevant channel is the probability F^e. Produced from decay of accelerated pions, 
these conventional (anti)neutrino beams suffer from an additional hurdle of an intrinsic j/g (i/g) contamination, 
which poses a serious problem of backgrounds. A beta-beam is comprised of pure Ve (or v,,). 
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Beam-related systematic imcertainties can be reduced to a large extent by working with 
a two detector set-up, one very close to the beam line and another serving as the far 
detector. The systematic uncertainties coming from the lack of knowledge of the neutrino- 
nucleus interaction cross-sections are another important source of error. These can be 
controlled to some degree by the near-far two detector set-up, but they cannot be canceled 
completely [91]. 

For a beta-beam, a variety of plausible set-ups have been proposed in the literature [46-67]. 
The proposal which poses minimal challenge for the beta-beam design, is the commonly called 
CERN-MEMPHYS project [47,58,59]. It proposes to use the EURISOL ion source to produce 
the radioactive source ions ^^Ne and ^He, and demands a Lorentz boost factor 7 ~ 100 for 
them, which can be produced using the existing accelerator facilities at CERN. The far detector 
MEMPHYS, a megaton water detector with fiducial mass of 440 kton, will have to be built 
in the Frejus tunnel, at a distance of 130 km from CERN. Another possible beta-beam set-up 
using water detector but higher boost factors and an intermediate baseline option was put forth 
in [62,63] (see also [66]). In these papers authors have used a high 7 ^^Ne and ^He beta-beam 
option at CERN and 440 kton fiducial volume water detector at GranSasso or Canfranc, which 
corresponds to L = 730 and 650 km respectively. Excellent sensitivity to ^13 and CP violation 
is expected [62,63] from this proposal. Set-ups with a neutrino beam from CERN to GranSasso 
or CanFranc [60], from CERN to Boulby mine [67] and from Fermilab [64] (L ~ 300 km) have 
also been proposed, and their sensitivity reach has been explored. Set-up with two sets of 
source ions with different boost factor for each set but with the same baseline was proposed 
in [66]. In [61] the authors consider the complementary situation where they take only one set 
of source ions, ^B and ^Li, with 7 = 350 and two different baselines, L = 2000 km and 7000 km. 
Very high gamma beta-beam options have been studied in [48,63,65]. The physics potential 
of the low energy beta-beam option was probed in [92]. A comparison of the physics reach 
among different beta-beam experimental proposals can be found in [46,56]. Here we would like 
to discuss the the possibility of measuring to a very high degree of accuracy the mixing angle 
^13 and sgn(Am|]^) aka, the neutrino mass ordering, in an experimental set-up [48-50,52,53] 
where a pure and intense multi-GeV z/g (using ^B ion) and/or (using ^Li ion) beta-beam 
with the Lorentz boost 7 between 250 and 650 is shot from CERN to the India-based Neutrino 
Observatory (INO) [25] where a large magnetized iron calorimeter (ICAL) detector is expected 
to come-up soon. 

Even if both neutrinos and antineutrinos are used, a serious complication with all long baseline 
experiments involving the golden Pg/i channel arises from discrete degeneracies which manifest 
in three forms: 



1. 
2. 
3. 



the {0is,6cp) intrinsic degeneracy [93], 



the (sgn(Am|]^), 6cp) degeneracy [94], 



the (6*23,^2 - 6*23) degeneracy [95]. 
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This leads to an eight-fold degeneracy [96] . with several spurious or "clone" solutions in addition 
to the true one and severely deteriorates the sensitivity of any experiment. It has been shown 
[96,97] that the problem of clone solutions due the first two types of degeneracies can be evaded 
by choosing the baseline of the experiment equal to the characteristic refraction length due to 
the matter inside Earth [10,96-99]. This special value goes by the colloquial name "magic 
baseline" [97]. As we will discuss in detail later, at this baseline the sensitivity to the mass 
hierarchy and, more importantly, 6*13, goes up significantly [52,97,100], while the sensitivity to 
5cp is absent. 

Interestingly, the CERN-INO distance of 7152 km happens to be tantalizingly close to the magic 
baseline (~ 7500 km). This large baseline also enhances the matter effect and requires traversal 
through denser regions of the Earth. Thus, for neutrinos (antineutrinos) with energies in the 
range 3-8 GeV sizable matter effects are induced if the mass hierarchy is normal (inverted). A 
unique aspect of this set-up is the possibility of observing near-resonant matter effects in the 
i^fi channel. In fact, to our knowledge, what we propose here is the only experimental 
situation where near-resonant matter effects can be effectively used in a long baseline experiment 
to study the neutrino mass matrix. We show in this work that the presence of this near- maximal 
Earth matter effect not only maximizes the sensitivity to the neutrino mass hierarchy, it also 
gives the experiment an edge in the determination of the mixing angle ^13. The increase in the 
probability Pg^ due to near-resonant matter effects, compensates for the fall in the beta-beam 
flux due to the very long baseline, so that one can achieve sensitivity to 6*13 and mass hierarchy 
which is comparable, even better, than most other proposed experimental set-ups. Therefore 
a beta-beam experiment with its source at CERN and the detector at INO could emerge as a 
powerful tool for a simultaneous determination of the neutrino mass hierarchy and ^13. 

In [48] such an experimental set-up was considered for the first time and the physics potential 
explored. It was demonstrated that both the neutrino mass hierarchy as well as ^^13 may be 
probed through such a set-up. The ions considered for the beta-beam in that work were the 
most commonly used ^^Ne (for z/e) and ^He (for 1/^)- One requires very high values of the Lorentz 
boost for these ions (7 ~ 10^) because the energy E of the beam has to be in the few GeV 
range for achieving near-resonant matter effects and to enable detection in the ICAL detector, 
which is expected to have a threshold of about 1 GeV. Such high values of 7, although possible 
in principle, might turn out to be very difficult to realize. Subsequently, two other ions, (for 
Ug) and ^Li (for 1/^), have been projected as viable options for the beta-beam source [72,73]. 
The main advantage that these ions offer is their substantially higher end-point energy, £"0. 
This allows one to access E ~ few GeV very easily with medium values of 7, that could be 
possible to achieve with either the existing CERN technology, or with the projected upgrades. 

In this chapter we will mainly present our results with and ^Li as the candidate ions for 
the beta-beam source at CERN and ICAL@INO as the far detector. We perform a detailed 
analysis of the projected data in a future CERN-INO beta-beam experiment and show our 
results with 7 ranging between 250 and 650. We begin in section 3.2 with a brief description of 
the beta-beam produced at the source. In section 3.3 we discuss the oscillation probability and 
highlight the importance of the magic baseline and the near-resonant matter effects in the Pg/i 
channel. Section 3.4 contains the results for the event rates expected in the ICAL detector and 
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a discussion addressing the issue of background rejection capability of ICAL. In section 3.5 we 
present the details of the statistical analysis used. In section 3.6 we focus on the sensitivity to 
the neutrino mass hierarchy. The sensitivity of the results on the mixing angle ^13 are presented 
in section 3.7. We end in section 3.8 with our conclusions. 



3.2 The Beta-beam Fluxes 



Very pure and intense z/g and/or Vf, beams can be produced by the decay of highly accelerated 
radioactive beta unstable ions, circulating in a storage ring [70]. This is what is called a "beta- 
beam" (see section 1.3.2) and was first proposed by Piero Zucchelli [45]. The selection of the 
beta unstable parent ion is determined by a variety of factors essential for efficiently producing, 
bunching, accelerating and storing these ions in the storage ring. A beta-beam comes with the 
following novel features : 



1. This neutrino beam would be very suitable for precision experiments because 
it is pure and mono-flavor and hence beam related backgrounds are almost 
absent. 

2. The neutrino spectrum depends only on the beta decay total end-point energy 

and the Lorentz boost of the radioactive ions 7. The spectral shape can 
therefore be very well determined which mcikes it almost free of systematic 
uncertainties. 

3. The flux normalization is determined by the number of useful ion decays in 
the straight section of the beam. 

4. The source is very intense. 

5. The beam divergence is controlled by the Lorentz boost 7. Hence by increasing 
7, we can produce a higher beam collimation and increase the beam intensity 
along the forward direction cx 7^. 

6. The z^e and z/g beta-beam fluxes could be produced simultaneously in distinct 
bunches from the storage ring. Therefore, in principle, it should be straight- 
forward to distinguish the neutrino from the antineutrino events at any de- 
tector using the nano-second time resolution that most detectors possess and 
the presence of magnetic fleld is not mandatory. 

7. It can be built using the existing CERN facilities. 
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3.2.1 The Beta-decay Spectrum 



In the rest frame of the radioactive ion the beta-decay spectrum is given by 

.f"^* =- ^(^0 - E*)E*^{Eo - E*y - ml (3.1a) 

dQ*dE* Anmlfr^ u ; v v u y e> \ / 

where nie is the electron mass, Eq the electron total end-point energy, E* and r are the neutrino 
energy and lifetime of the decaying ion respectively in the latter's rest frame and / is defined 
in the fashion 



five) = 



60?/, 



Vl-ye'(2 - 9yl - 8yt) + 15y,^Log 



1- 



Ve 



(3.1b) 



where ye = rrie/EQ. 

In the lab frame, the flux of the unoscillated beta-beam at the detector is given by 

_ 1 iV^ 1 



AttL"^ m-g / 7(1 — /5 cos ( 



-{Eo - E*)E*WiEo-E*)^-ml 



(3.2) 



where L is the distance between the source and detector, A^^ are the number of useful decays 
in the storage ring per unit time, 6 is the angle between the neutrino flight direction and the 
direction in which the ions are boostecH and 7 is the Lorentz boost such that E* = 7-E(l — 
(3 cos 6), E being the neutrino energy in the lab frame. The maximum energy of neutrinos 
produced by a beta-beam with a Lorentz factor 7 is given by 



E„ 



{Eo - rrie 



For the on- axis 



7(1 -/5 cos ^)' 

0) neutrino beam, the above expression looks like the following 



(3.3a) 



E„ 



7(Eo-me)(l + /3) 



on — axis 



(3.3b) 



3.2.2 Candidate Ions for the Beta-beam 

From Eq. 03.21) it is seen that the total flux and energy of the beta-beam at the far detector 
depends mainly on the end point energy Eq of the beta-decay ion and the Lorentz boost factor 
7. The flux increases as both Eq and 7 increase. Larger 7 results in better collimation of the 
impinging flux, thereby increasing the statistics. The spectrum also shifts to larger energies as 
£^0 and 7 increase. The neutrino cross section in the detector increases with energy, so for the 
same total flux a harder spectrum further enhances the statistics. Since the flux falls as 1/L^ as 

^We work with the on-axis flux for which 9 = 0. 



33 




I 

u 

to 
o 




1.5 2 2.5 3 3.5 

Energy (GeV) 



1.5 2 2.5 3 3.5 4 
Energy (GeV) 



Figure 3.1: The unoscillated beta-beam flux spectrum arriving at ICAL@INO. The upper panels 
are for ®B (left panel) and ^Li (right panel), while lower panels are for ^^Ne (left panel) and 
^He (right panel). Due to the finite energy width (~ 1.5 MeV) of the 2+ first excited state of 
^Be, the maximum neutrino energy available in the decay of ®Li and may be increased by 5 
-6%. 
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Figure 3.2: Left panel shows the total unoscillated flux in yr~^m~'^ expected at INO, as a 
function of the Lorentz factor 7. The solid (dashed) line corresponds to (^Li) and we have 
assumed 1.1 x 10^^ (2.9 x 10^^) useful ions decays per year. Right panel shows the energy at 
which the unoscillated flux peaks, as a function of 7. 



the source-detector distance increases, high values of either £'0 or 7 or both are needed to have 
sizable number of events at the far detector. Eq is an intrinsic property of the decaying ion 
while 7 is restricted by the nature of the accelerators and the machine design. The other desired 
properties which the ion should have include high production yield, large decay fraction, and 
a lifetime that is long enough to allow the ions to be adequately accelerated. It is also easier 
to store larger number of lower-Z isotopes in the storage ring [101]. The characteristic features 
of the four different ions which have been discussed in the literature as possible candidates for 
the beta-beam have already been shown in Table [TTT] (section 1.3.2). While ^^Ne and are 
emitters (producing a Ue beam), ^He and ^Li are /3~ emitters (producing a z/g beam). They 
all have comparable lifetimes, conducive to the requirements necessary for the beta-beam, very 
high (or maximal) decay fraction and very low A/Z ratio. The point to be noted that the 
end-point energies, Eq, for ^B and ^Li are much larger than those for ^^Ne and ^He. 

In this work we are interested in the physics potential of an experimental set-up with a beta- 
beam at CERN and a large magnetized iron detector in India. For such a large baseline, one 
needs a very high value of 7, both to cross the detector energy threshold as well as to get 
reasonable statistics in the detector. In particular, with ^He (or ^^Ne) as the source ion, one 
needs 7 ^ 1000. Such high values of 7 can only be achieved by using the LHC itself. As noted 
above, Eq for the alternative options for the radioactive ion source (^B and ^Li) is higher by a 
factor of more than 3. This means that it should be possible to produce high intensity, high 
energy beams with *B and ^Li for a much lower boost factoiE]. In the low 7 design, the standard 
numbers taken for the ^^Ne and ^He are 1.1 x 10^^ (z/g) and 2.9 x 10^^ (ue) useful decays per 
year, respectively [75]. Earher, only these ions were considered because it was believed that 

^The loss in collimation is not significant as we show later. 
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could not be produced with the standard ISOLDE techniques. Since most exercises focused 
on observation of CP-violation, it was necessary to have both i/g and i>e beams with similar 
spectra, so '^Li (though considered in [101]) was also generally ignored. Interest in both these 
ions have been rekindled recently [66,72,73], as it appears that having intense and ^Li fluxes 
should be feasible using the ionization cooling technique [72,73]. In what follows, we will vary 
the value of 7 to test the physics potential of the CERN beta-beam INO-ICAL set-up and, 
following the current practice, assume that it is possible to get 2.9 x 10^^ useful decays per 
yeaiQ for ^Li and 1.1 x 10^^ for ^B for all values of 7. Note, however, that new ideas suggest 
luminosities higher even by a factor of ten or so, depending on the isotopes used, by using a 
recirculating ring to improve the performance of the ion source [72,73]. Studies have shown 
that it is possible to accelerate ^He to 7 ^ 250 with the existing facilities at CERN, while 
7 = 250 — 600 should be accessible with the "Super-SPS" , an upgraded version of the SPS with 
super-conducting magnets [60,62,63]. The Tevatron at FNAL could in principle also be used 
to produce a beta-beam with 7 < 600. 

The upper panels of Fig. 13.11 depict the unoscillated '^B (left-hand panel) and ^Li (right-hand 
panel) beta-beam flux expected at ICAL@INO, for four different benchmark values of 7. The 
lower panels of Fig. 13.11 show the corresponding spectra for the ^^Ne (left-hand panel) and ^He 
(right-hand panel) beta-beam. Note that even though apparently it might seem from the figures 
that the (j^e) fiux is larger for ^^Ne (^He), in reality, for a given 7, the total fiux is given by 
the area under the respective curves. One can easily check that for a fixed 7 this is same for 
both the ions as expected, since we have assumed equal number of decays for both ^^Ne and 
^B for z/e (^He and ^Li for Also note that even though the total Ve (t'e) fiux remains the 
same for both the ions for a given 7, we would expect that the number of events produced in 
the detector is much enhanced for the ^B (^Li) beam since the energy of the beam is larger and 
the CC cross sections increase with the neutrino energy. In the left panel of Fig. 13.21 we show 
the energy integrated total number of unoscillated neutrinos in units of yr~^m~'^ arriving at 
ICAL@INO, as a function of the Lorentz boost 7. The solid (dashed) line corresponds to '^B 
(^Li). The figure shows that the energy integrated fiux arriving at the detector increases almost 
quadratically with 7. Note that with the same accelerator, the Lorentz boost acquired by '^B 
is 1.67 times larger than that by '^Li, determined by the charge to mass ratios of the ions. The 
right panel of Fig. 13.21 depicts the energy at which the unoscillated fiux peaks, as a function of 
the Lorentz boost 7. It turns out that this peak energy is roughly half the maximum energy of 
the beam, which is given as -Emax — 2(-Eo ^ f^e)l (see Eq. I3.3ap . 



3.3 Neutrino Propagation and the "Golden Channel" 



In this section we will very briefiy review some issues related to the neutrino oscillation prob- 
ability. 



"^Wherever not explicitly mentioned, these reference numbers of useful ion decays for v,, and are chosen. 
^Of course, the on- axis flux increases with 7 because of better coUimation of the beam. 
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3.3.1 The "Golden Channel" (z/g ^ i^^d 

The expression for Pg/x in matter [10, 87], up to second order terms in the small quantities 613 
and a = Amj^/Aml^, is given bj@ [88, 102] 

9 ^ sin2[(l-i)A] 
Pe^ - sin" 26^3 sin" 623- ^ ^ 



± a sin 2^13 ^ sin 5 sin(A) 



+ a sin 2^13 ^ cos 5 cos(A) 



1 - A)2 

sin(iA) sin[(l - A) A] 



A (i-A) 

sin(iA) sin[(l - i)A] 



where 



+ a' cos' sin' 29,2^^^^-^^, (3.4) 

A-^. (3.5) 

^ = COS 6*13 sin 26*12 sin 26*23, (3-6) 



A 



A = (3.7) 
Ami, 

Above, A is the matter potential (see Eq. I2.22p . The Eq. 13.41 has been derived under the 
constant matter density approximation. The first term of Eq. ( 13. 4p can be used to extract 
information about the value of 6*13. This is also the term which has the largest Earth effect and 
this effect of matter can be used to determine the sign of Am\,. In Eq. (13.41) the second term 
has the CP violating part. This term is positive for neutrinos and negative for antineutrinos. 
The third term, though 5cp dependent, is CP conserving. The last term is independent of both 
9i3 and 5c p and depends mainly on the solar parameters Am\, and 6*12. In very long baseline 
experiments, this term has sizable matter effects and if the true value of 6*13 turns out to be 
zero (or nearly zero), this would be the only surviving term in Pg^ which could still be used to 
study matter enhanced oscillations [103]. While we will use this formula to discuss our results 
in some cases, our simulation is based on the exact probabilities. 



3.3.2 Eight-fold Degeneracy 

Despite its advantage in determining the most interesting oscillation parameters, this channel 
is, however, rife with the problem of "degeneracies" - the (^13 - 5cp) intrinsic degeneracy [93], 



^This particular low order expansion of the transition probability is valid only in the range of L and E 
where the resonance condition A = 1 is never reached. For the E and L range that we consider in this work, 
this condition is satisfied and Eq. (|3.4p fails, as we discuss later in this section. Nonetheless, we present this 
expression to illustrate the effect of the magic baseline. 
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Table 3.1: Possibilities for long baselines from CERN. In the third column, the distances (in 
km) from CERN to various existing/proposed underground neutrino observatories spread over 
the entire world are given. The respective dip angles are also mentioned in the last column. 
The distance of INO from CERN is 7152 km with a dip angle of 34.0°. 



the {sgn{^m\^) - 6cp) degeneracy [94] and the octant of 623 degeneracy [95] - leading to an 
overall eight-fold degeneracy in the parameter values [96], of which, obviously, only one is true. 
These degeneracies always deteriorate the performance of an experiment. 

3.3.3 Remedy with "Magic" Baseline 

To tackle these problems, there are various schemes suggested in the literature which can be 
broadly categorized in the following fashion : 

1. combine data from several experiments observing the golden channel, but with different 
baselines L and neutrino energies E [93, 104, 105], 

2. combine data from accelerator experiments observing different oscillation channels [51, 
106,107], 

3. combine the golden channel data with those from atmospheric neutrino experiments [47, 
108], 

4. combine the golden channel data with those from reactor antineutrino experiments [109] , 

5. kill the spurious clone solutions at the "magic baseline" [97]. 
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Figure 3.3: The upper panel shows the constant line average density of the Earth for various 
baselines calculated using the Preliminary Reference Earth Model (PREM) [110] and the lower 
panel depicts the corresponding resonance energy (where 6"^ = '/r/4) at those baselines. The 
blue dashed vertical line is drawn at the CERN-INO baseline. 



Now let us see how the concept of "magic" baseline works. A particularly interesting scenario 
arises when the condition 

sin(iA) = (3.8) 

is satisfied. In such an event, the last three terms in Eq. (13. 4p drop out and the channel 
enables a clean determination of ^13 and s(y'n(Am|]^). 

Since A A = AL/4E by definition, the first non-trivial solution of Eq. 13.81 comes in the form 
pL = ^J^-K IGpYe- This gives 

r ^ 3in^i - 32725 . (3.9) 
[gm/cm'^J [km] 

From the upper panel of Fig. 13.31 one can immediately observe that the above condition is 
satisfied for the "magic baseline" [97] 

L^,gi, ~ 7690 km, (3.10) 

where the constant line average density of the Earth0 estimated from the PREM profile [110] 
comes out to be 4.25 gm/cm^. An enlightening discussion on the physical meaning of the magic 
baseline is given in [98]. 

Table [3TT] depicts the distances (in km) and the corresponding dip angles (in degree) from CERN 
to various existing/proposed underground neutrino observatories spread over the entire world. 



^We use the full PREM profile in our simulation without any approximations. 
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Benchmark Values 


la estimated error 


lAm" (true)l = 2.5 x lO^^ eV" 


(T(|Amii|) = 1.5% 


sin" 2^23 (true) 1 = 1.0 


a (sin" 2^23) = 1% 


Am^i(true) = 8.0 x 10"^ eV" 


a(Am^i) = 2% 


sin"^i2(true) = 0.31 


(T(sin" Ou) = 6% 


p(true) = 1 (PREM) 


a(p) = 5% 



Table 3.2: Chosen benchmark values of oscillation parameters and their la estimated errors. 
The last row gives the corresponding values for the Earth matter density. 



The CERN-INO distance corresponds to L = 7152 km, which is tantalizingly close to the magic 
baseline. We therefore expect that the beta-beam experiment we consider here should give an 
essentially degeneracy-free measurement of both 613 and sgn{Am'^i). 



3.3.4 Near-Resonant Matter Effects 

The large CERN-INO baseline, of course, results in very significant Earth matter effects in the 
Pefj, channel. In fact one can readily see from the upper panel of Fig. I3.3l that for the baseline of 
7152 km, the line average Earth matter density calculated using the PREM profile is pav = 4.17 
gm/cm^, for which the resonance occurs at (see the lower panel of Fig. 13.31) 

^ |Am^|cos2g,3 (3^^^) 

2V2GFNe 

= 7.45 GeV , (3.12) 

for |Am|i| = 2.5 x 10~^ eV" and sin" 26113 = 0.1. One can check (cf. Fig. [XT]) that this is 
roughly in the ballpark where we expect good enough beta-beam flux using and ^Li ions 
with the 7 in the range 350-650. The resonance energy is oc cos 2^13 and therefore it depends 
on ^13 mildly. At the CERN-INO baseline, the resonance energy can vary within a very small 
range of 7.1 to 7.85 GeV depending on the present allowed span of 6*13. 



3.3.5 The Benchmark Oscillation Parameters 

Note that the low order expansion of the probability given by Eq. fl3.4p is valid only for 
values of E and Earth matter density p (and hence L) where flavour oscillations are far from 
resonance, i.e., A <^ 1. In the limit A ~ 1, one can check that even though the analytic 
expression for given by Eq. (13. 4p remains finite, the resultant probability obtained is 
incorrect [111]. We reiterate that Eq. fl3.4p was presented only in order to elucidate the 
importance of the magic baseline. For all the numerical results presented in this thesis, we 
calculate the exact three generation oscillation probability using the full realistic PREM [110] 
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Figure 3.4: The transition probability Pe/j, as a function of E for four values of the baseline L. 
The band reflects the effect of the unknown 6cp- The dark (red) shaded band is for the NH 
while the light (cyan) shaded band is for the IH. We have taken sin^ 26'i3 = 0.1 and for all other 
oscillation parameters we assume the benchmark values given in Table 13.21 



profile for the Earth matter density. Unless stated otherwise, we have generated our simulated 
data for the benchmark values in the first column of Table 13.21 These values have been chosen 
in conformity with the status of the oscillation parameters in the light of the current neutrino 
data [78]. The values of sin^ 2^13 (true), 5c'p(true) and mass hierarchy which are allowed to 
vary in our study, will be mentioned wherever applicable. When we fit this simulated data, we 
marginalize over all the oscillation parameters, the Earth matter density, as well as the neutrino 
mass ordering, as applicable. We expect to have a better knowledge of all the parameters 
mentioned in Table 13.21 when the Beta-beam facility comes up. In particular, we assume that 
the la error on them will be reduced to the values shown in the second column of Table [3^ [90. 
112,113]. Therefore, we impose "priors" on these quantities, with the corresponding la error. In 
the fit, we allow for a 5% uncertainty in the PREM profile and take it into account by inserting 
a prior and marginalizing over the density normalization. We analyze the full spectral data 
from both the neutrino C^B) and antineutrino (^Li) run expected in the CERN-INO beta-beam 
set-up. 
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Figure 3.5: The dark (red) shaded band is the same as in Fig. 13.41 The hght (green) shaded 
band shows the corresponding Pe^ for sin^26'i3 = 0.05. Values of all the other oscillation 
parameters are same as in Fig. 13.41 and the hierarchy is assumed to be normal. 



3.3.6 Phenomenology with P^^ 

The exact neutrino transition probability using the full PREM density profile is given in Fig. 
13.41 as a function of the neutrino energy, for four different baselines (four panels) . We allow 
6c p to take on all possible values between 0-360 degrees and the resultant probability is shown 
as a band, with the thickness of the band reflecting the effect of 6cp on Pe^j- The figure is 
drawn assuming the benchmark values of the oscillation parameters given in Table 13.21 and 
sin^ 26'i3 = 0.1. We show the probability for both the NH (dark band) as well as the IH (light 
band). As discussed in detail above, for L = 7500 km, which is close to the magic baseline, the 
effect of the CP phase is seen to be almost negligible, while for all other cases the impact of 6cp 
on Pefj, is seen to be appreciable. In fact, for L = 1000 km, the probability corresponding to the 
NH and IH become almost indistinguishable due to the uncertainty arising from the unknown 
value of (5cp. As the baseline is increased. Earth matter density increases, enhancing the impact 
of matter effects. The probability for NH is hugely enhanced for the neutrinos, while for the IH 
matter effects do not bring any significant change. This difference in the predicted probability, 
evident in the panels corresponding to L = 4000, 7500 and 10000 km of Fig. 13.41 can be used 
to determine the neutrino mass ordering. 
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L (km) 

Figure 3.6: The values of sin^ 2^13 needed for the maximal matter effect (from Eq. 13.181 with p 
= 0) at different baselines using the constant line average density of the Earth. The blue dashed 
horizontal line shows the present upper bound on sin^ 2^13 which predicts that the maximal 
probability (= 0.5) can only be achieved for L > 7880 km (see the black dot-dashed vertical 
hne). 



In Fig. 13.51 we display the dependence of the neutrino probability Pg/x on the mixing angle ^13 
for four different baselines. The dark bands, as in Fig. 13.41 are for NH and sin^2^i3 = 0.1 with 
full variation of Sep, while the light bands are for NH and sin^2^i3 = 0.05. The impact of 
matter effect in increasing the ^13 sensitivity of a given experimental set-up is evident from the 
figure. The ^13 sensitivity for L = 1000 km can be seen to be much weaker than for the other 
cases, since matter effects are smaller. However, the most striking feature seen in Fig. 13.51 is 
the effect of the magic baseline in enhancing the sensitivity of the experiment to 61^. The figure 
clearly shows that the difference in the predicted Pg^ for the two values of sin^ 26'i3 is largest 
for L = 7500 km (since effect of 6cp is the least) and thus an experiment at this baseline is 
most suitable for probing ^13. Figs. 13 . 41 and [ 3.51 therefore reinforce our choice of the near-magic 
baseline as one of the best options for determining the neutrino mass hierarchy and ^13, since 
both these parameters are directly related to large matter effects and the uncertainty of Sep 
could prove to be a hindrance in their measurement at non-magic baselines. 



3.3.7 One Mass Scale Dominance 

Now let us discuss the "Golden" channel oscillation probability in the approximation where 
the solar mass squared difference Amg^ can be neglected. This is known as one mass scale 
dominance (OMSD) approximation. Under this condition with the constant matter density 
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approximation, the probability Pg^ can be expressed as 

= sin^ ^23 sin^ 2^^ sin^ [l .27 {Aml^)'^ L / E] , (3.13) 

where sin^ 26*]^ and [/\m\^)^ are given by Eq. 12.251 The probabihty Pe^ would be largest (= 
0.5, assuming 6'23 = 7r/4) if the conditions 

sin2 2^™ = l (3.14) 

and 



sin^ 



[l.27(Am2j"^L/E] = 1 (3.15) 

are satisfied simultaneously. The first condition is achieved at the resonance energy. The second 
condition gives the energy where the (Amg^)™ driven oscillatory term is maximal, 

_ 1.27(Am^,)-L 

- (2p + l)vr/2 (3-16) 
Maximum matter effect is obtained when [114-116] 

which gives (using Eqs. I2.25l fc [3.1ip . 

/ T^max (2p+ l)7r X 5.18 X 10^ ^ ,3 

(pL) = — km gm/cm"*. (3.18) 

tan 2^13 

Although both Eres and E^^^^ depend on the value of Amg^, the distance at which we get 
the maximum matter effect is independent of Amg^^. However, it is controlled by very 
sensitively. In Fig. 13. 6[ we have plotted the values of sin^ 2^13 required for the maximal matter 
effect at different baselines using the constant line average density of the Earth as predicted 
by the PREM profile. For p = 0, we can have maximal matter effect only if L > 7880 km 
because of the upper bound of 0.19 on sin^ 2^13 [31,78]. At L = 7152 km, for maximal matter 
effect sin^ 2^13 ~ 0.23, which is already ruled out by the present constraint. Therefore at the 
CERN-INO baseline, maximal matter effect cannot be achieved and the oscillation probability 
peaks roughly at 6 GeV (see the Figs. 13.41 & 13.51) instead of 7.45 GeV 



3.4 Event Rates in ICAL@INO 



3.4.1 The ICAL Detector at INO 

The proposed large magnetized iron calorimeter at the India-based Neutrino Observatory [25] is 
planned to have a total mass of 50 kton at startup, which might be later upgraded to 100 kton. 
The left panel of Fig. 13.71 shows the location of the INO facility which is expected to come up 
at PUSHEP (lat. North 11.5°, long. East 76.6°), situated close to Bangalore in southern India. 



44 




Figure 3.7: Left panel shows the location of the India-based Neutrino Observatory. The 
schematic view of the 50 kton iron calorimeter detector consisting of three modules each having 
140 layers of iron plates is given in the right panel. These figures are obtained from the official 
website of INO (http://www.imsc.res.in/~ino/). 



The distances and the corresponding dip angles from INO to various existing accelerator based 
laboratories are given in Table 13.31 The JAERI (Japan Atomic Energy Research Institute) to 
INO distance is 6477 km with a dip angle of 30.5°. The distances from INO to CERN and 
RAL (Rutherford Appleton Laboratory) are 7152, and 7653 km, respectively. These baselines 
are quite close to the magic baseline. In this sense, the location of INO is ideal for futuristic 
long baseline experiments. The BNL (Brookhaven National Laboratory) and the FNAL are 
the two major accelerator based facilities in US and they constitute the baselines of 11081, and 
11306 km, respectively from INO with the dip angles > 60°. The first phase of this detector 
(starting around 2012) will focus on atmospheric neutrino measurements. The ICAL detector 
will have a modular structure with a total lateral size of 48m x 16m, divided into three modules 
of 16m X 16m each (see the right panel of Fig. 13. 7p . Each of these modules will have 140 
horizontal layers of ~ 6 cm thick iron plates, separated from each other by a gap of ~ 2.5 
cm to hold the active detector material, giving a total height of 12 m for the full detector. 
The active detector elements will be 2 cm thick glass resistive plate chambers (RFC) and will 
be filled with a suitable gas mixture, which will be recycled with approximately one volume 
change per day. An internal magnetic field of ~ 1.3 Tesla would be applied over the entire 
detector. The detector will be surrounded by an external layer of scintillator or proportional 
gas counters which will act both as veto to identify external muon backgrounds as well as to 
identify partially contained events. 

We are interested in measuring the golden channel probability Pe^t- Since we have a z/g (z/g) 
flux in the beam, we need a detector which is sensitive to muons (antimuons). The detector 
should have a suitable energy threshold, depending on the energy spectrum of the beta-beam. 
In addition, it should have a good energy resolution and low backgrounds. According to the 
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X-Zc LcC LUi ±jUCc1LHJ11 




Til ci" n r»o 






JAERI 


6477 km 


30.5° 




CERN 


7152 km 


34.0° 


INO 


RAL 


7653 km 


36.8° 




BNL 


11081 km 


60.3° 




FNAL 


11306 km 


62.4° 



Table 3.3: The distances from INO to various accelerator based laboratories spread over the 
entire world. The corresponding dip angles are also noted. 



Total Mass 


50 kton 


Energy threshold 


1 GeV 


Detection Efficiency (e) 


80% 


Charge Identification Efficiency (//d) 


95% 


Energy Resolution (a) (GeV) 


0.15E(GeV) 


Bin Size 


1 GeV 


Background Rejection 


0.0001 


Signal error 


2.5% 


Background error 


5% 



Table 3.4: Detector characteristics of ICAL@INO used in the simulations [117]. The bin size 
is kept fixed, while the number of bins is varied according to the maximum energy. 



ongoing detector simulation study performed by the INO collaboration, the detector energy 
threshold for /i^ is expected to be around ~ 1 GeV and charge identification efficienc}{§ will 
be about 95%. The detection efficiency of ICAL after cuts is expected to be about 80%. In 
what follows, we will present our numerical results assuming an energy threshold of 1 GeV, 
detector charge identification efficiency as 95% and detection efficiency as 80% for /i^ (cf. Table 
13. 4p . However, we will also show the impact of changing the threshold. We have seen that our 
results remain unaffected if the energy threshold is raised to 2 GeV for the entire range of 
assumed Lorentz boost factor 7 = 250 — 650. For 7 > 350 the threshold can be even 3 GeV, 
while for 7 > 500 one can work with an energy threshold of 4 GeV, without changing the final 
results. This can be seen from Fig. 13. Ij for 7 = 350(500), the majority of neutrinos arriving at 
INO-ICAL would have > 3(4) GeV. The energy resolution of the detector is expected to be 
reasonable and we assume that the neutrino energy will be reconstructed with an uncertainty 
parameterized by the Gaussian energy resolution function with aE = 0.15E, where E is the 
energy of the neutrino and as is related with the half width at half maximum (HWHM) in 
the fashion HWHM = 1.1 7cr. We will present and compare the sensitivity of this experimental 

*In case of a beta-beam, the neutrino and antineutrino events produced in distinct bunches can be dis- 
tinguished from each other using the nano-second time resolution of the detector. The charge identification 
efficiency is incorporated since that helps in reducing the neutral current backgrounds. 
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Figure 3.8: The expected number of events in 5 years running time, as a function of sin^ 2^13. 
The value of 7 and the hierarchy chosen corresponding to each curve is shown in the figure 
legend. 



set-up with and without the full spectral analysis using both the neutrino (^B) and antineutrino 
(*Li) run. 



3.4.2 Oscillation Signal at ICAL@INO 

The z/^ induced /i~ event spectrum at INO is estimated using 

N, = T Tin fine dE ' dEAKiE)(T,,{E)RiE,EA)Pe>.iE), (3.19) 

Jo Jet'^ 

where T is the total running time (taken as five years), (l)y^{E) (given by Eq. 13. 2p is the 
unoscillated beta-beam flux at the detector, e is the detector efficiency, are the number 
of target nucleons in the detector, is the charge identification efficiency and R{E,Ea) is 
the energy resolution function of the detector, for which we assume a Gaussian function. The 
quantities E and Ea are the true and reconstructed neutrino energy respectively. E"^^^ is the 
maximum energy of the neutrinos for a given Lorentz factor 7 (given by Eq. fl3.3ap ) and Pg^ 
is the z/e — > oscillation probability. For muon events, a^^ is the neutrino interaction cross- 
section. The expression for the /i+ signal in the detector from a beta-beam flux is given 
by replacing (f)^^ by a^,^ by Up^ and Pg^i by Pep,. For the neutrino-nucleon interaction we 
consider quasi-elastic scattering, single-pion production and deep inelastic scattering and use 
the cross sections given in the Globes package [118] which are taken from [119,120]. 



1. 



Event Rates vs. sin^ 29 
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Fig. 13.81 portrays the number of events expected in INO-ICAL from 5 years exposure 
of an (for v^., shown in the left-hand panel) or ^Li (for v^., shown in the right-hand 
panel) beta-beam from CERN. The expected number of events is presented as a function 
of sin^ 2^13 for both NH and IH for three benchmark values of the boost factor 7. Large 
resonant matter effects in the neutrino channel for NH drives the number of expected 
events to very large values, compared to what would be expected for IH. Similarly, in the 
antineutrino channel we have resonant matter effects for IH and the number of predicted 
events is many times larger than for NH|^. This difference in the number of events is seen 
to increase with sin^ 2^13. 

From the left-hand panel of Fig. I3.8l it can be noted that for 7 = 500 and sin^ 2^13 = 0.05, 
the predicted number of neutrino events for NH is 137, while that for IH is only 4. This 
implies that if the NH was true, we could comprehensively rule out the wrong IH using 
the neutrino run. If IH is true then antineutrino run will be very useful (see the right 
panel of Fig. 13. 8p to rule out the wrong NH. Note that while the interaction cross section 
for the i/gS are smaller, the flux itself is larger owing to the (assumed) larger number of 
decays per year for ®Li. Thus the statistics expected in both the neutrino as well as the 
antineutrino channel is comparable. As discussed above, resonant matter effects in the 
neutrino (antineutrino) channel for the normal (inverted) hierarchy, results in substantial 
enhancement in the observed number of events. In particular, we note from Fig. 13.81 
that the number of events in the near-resonant channels depends strongly on the value of 
sin^2^i3, since the extent of matter effects is dictated directly by 6'i3. We can see from 
the figure that the dependence of the event rate on sin^ 26'i3 is much enhanced due to 
matter effects. 



Iso-event Curves in the sin^ 20i^-5cp plane 



Fig. 13.91 shows the effect of 5c p on the variation of the measured rate with the value of 
the mixing angle ^13, for four different baselines 732 km, 4000 km, 7152 km and 10000 
km. The left-hand panel shows the results for the Ve beta-beam assuming a normal mass 
hierarchy, while the right-hand panel gives the same for the z/g beta-beam flux with IH. 
For all cases in any panel, a 50 kton magnetized iron calorimeter as the far detector 
and the same flux created at the source have been considered. Each curve gives the 
sets of values of {sin^ 26'i3, ^cp} which give the same observed rate in the detector as 
the set {sin^ 2^13 = 0.04, 5(7p = 0°}. In other words, if the true value of sin^ 26*13 and 
5cp were 0.04 and 0° respectively, then every point on a given curve would also be a 
solution for that experiment. It is clear from this figure that combining results from 
experiments at different baselines helps solve/reduce the problem. However, the most 
important issue exemplified here is the fact that for the baseline 7152 km, which is the 
CERN-INO distance, the effect of the unknown value of 5cp on the measurement of ^13 
and the mass hierarchy, is very small. This happens because this distance corresponds 
to a near-magic baseline for which, as noted earlier, the 5c p dependent terms are almost 
vanishing. 



^In fact, matter effects are seen to suppress tlie oscillation probability for neutrinos (antineutrinos) compared 
to that in vacuum when the hierarchy is inverted (normal) . 
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Figure 3.9: The iso-event curves in the 6cp-siia^ 26'i3 plane for four basehnes are shown in the 
figure. The true values of 6cp and sin^ 26*13 are assumed to be 0° and 0.043 respectively. Left- 
(right-) hand panel is for the Ue {ve) beta-beam. The assumed hierarchy is mentioned in the 
figure. 



3. Event Rates vs. L 



Fig. 13.101 depicts the number of events expected in five years as a function of the baseline 
L, if we run the experiment in the neutrino mode with 7 = 500. A similar figure is 
expected for the antineutrino beam. The upper black hatched area shows the events for 
sin^ 2^13 = 0.05 and the lower red hatched area corresponds to sin^ 26*13 = 0.01. For each 
baseline L, the range covered by the hatched area shows the uncertainty in the expected 
value of the number of events due to the completely unknown 5c p-, which could take any 
value from to 27r. The baseline L where the width of this band reduces to (almost) zero 
is the magic baseline. We see from the figure that the magic baseline appears at about 
L ~ 7500 km. Note that while for sin^ 26*13 = 0.01 the magic baseline is very clearly 
defined with the CP dependence going completely to zero, for the higher value of sin^ 2^13 
of 0.05, the "magic" is not complete. The reason for this anomaly can be traced to the 
fact that Eq. (13. 4p was derived for only very small values of 6*13. For larger values of 
this angle, higher order terms become important. These terms might depend on 5c p and 
remain non-zero even at the magic baseline. 



Event Rates vs. 7 



Figs. 13.111 and 13.121 demonstrate the impact of the magic baseline on the mass hierarchy 
and 6*13 sensitivity respectively. In each of the four panels of both the figures the expected 
events in five years as a function of 7 have been shown. Each panel is for a certain fixed 
value of L, shown in the corresponding panel. In Fig. I3.11l the black hatched area shows 
the band for NH while the open band delimited by the dashed red lines are for the IH. 
As in Fig. I3.10[ the band corresponds to the uncertainty in the event rate due to the 
unknown 5c p- The effect of the uncertainty of 5c p almost vanishes for L = 7500 km 
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Figure 3.10: Total number of expected events in five years as a function of the baseline L for 
the source with 7 = 500 and for two values of sin^ 2^13 and assuming that the NH is true. 
The hatched areas show the expected uncertainty due to the CP phase. 



which is very close to the magic baseline. It can be seen that for the smaller baseline 
L = 1000 km, NH and IH predictions are largely overlapping, making it almost impossible 
for these experiments to give sensitivity to the mass hierarchy unless sin^ 29i3 turns out 
to be extremely large and Sep favorable. The hierarchy sensitivity is expected to improve 
as we go to larger baselines and this is reflected from the two bands for NH and IH 
separating out. It turns out that because the matter effects are very large for the magic 
baseline and effect of CP uncertainty is zero, this baseline gives the best sensitivity to the 
mass hierarchy. For L larger than magic, matter effects are higher but the flux is lower, 
while for L lower than magic, flux is higher but the matter effects are lower. For both 
above and below the magic baseline, the effect of 6cp is expected to further reduce the 
sensitivity. This is particularly true for the lower L baselines. Fig. 13.121 shows the bands 
for NH but with two different choices for sin^2^i3. Here the effect of the magic baseline 
is seen even more clearly. 

Fig. 13.131 shows as a function of 7, the number of events expected in five years in the 
CERN-INO beta-beam set-up. The solid (dashed) lines are for neutrino (antineutrino) 
events, with the thick line showing the event rate for NH (IH) while the thin line is for the 
IH (NH). We have assumed sin^ 26'i3 = 0.01 and 6cp = 0. One point which is transparent 
from this figure and which will be very relevant in understanding the behavior of the 
CERN-INO beta-beam set-up is the following: For a given value of ^13 and for NH (IH), 
we expect a large number of events in the neutrino (antineutrino) channel and almost 
negligible events in the antineutrino (neutrino) channel. This means that for NH (IH) it 
will be the neutrino (antineutrino) channel which will be statistically more important. 
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Figure 3.11: Total number of events as a function of 7 for the B source, for different values 
of L are shown in the four panels. The black hatched area shows the uncertainty range due 
to the CP phase when NH is true, while the area between the maroon dashed lines shows the 
corresponding uncertainty when IH is true. For all cases we assume sin^2^i3 = 0.05. 



3.4.3 Background Rejection 



Some technical issues pertaining to rejection of the beam related and atmospheric neutrino 
backgrounds will be discussed here. 



Beam Related Backgrounds 



The possible backgrounds in a z/g beta-beam experimenio using /i as an oscillation signal 
come from NC events such as 

v^ + d{u) v^ + d{u) (3.20) 

+ d{u) i'x + d{u) + qq (3-21) 

and z/g charged current events 

z/g + (i + u or c (Cabibbo suppressed) (3.22) 

z/e + (i + u + qq (3.23) 

"'^''The discussion concerning the v^. beta-beam is similar and hence is not repeated here. 
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Figure 3.12: Total number of events as a function of 7 for the source, for different values of L 
are shown in the four panels. The black hatched area shows the uncertainty range in the events 
due to CP phase when sin^ 2^13 = 0.05, while the red hatched area shows the corresponding 
uncertainty when sin^ 2^13 = 0.01. For all cases we assume NH to be true. 



The quarks in the final state could produce mesons as a part of the hadronic junk. These 
mesons might constitute backgrounds in the following way : 

1. Some energetic mesons can give long tracks inside the detector which can mimic the muon 
tracks. 

2. These mesons can decay in flight producing secondary muons, a possible source of back- 
ground. 

As discussed earlier, ICAL@INO will have 6 cm thick iron plates. Such a dense tracking detector 
will have excellent muon/pion and muon/electron separation capability in the energy range we 
are working with. The simulations carried out by the INO collaboration have shown that after 
the standard kinematical cuts are imposed, the electrons do not give any signal at all, while in 
99% of the cases, the pions and kaons for the energy range of interest to us get absorbed very 
quickly in the iron via strong interactions and therefore do not hit enough RPCs to give a signal. 
At the energies of a beta-beam considered here, production cross section of D mesons (also 
Cabibbo suppressed) is small and they do not constitute a problem for the experiment. The 
associated strange or charm production is also highly suppressed at these energies. In addition, 
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Figure 3.13: Total number of events as a function of 7 for the *B (solid lines) and the ^Li 
(dashed lines) sources. Results for both normal and inverted hierarchies are shown. 



the fact that the detector will have a charge identification capability means that secondary 
fi^ produced can be safely discarded, reducing the background even further. Therefore, in the 
analysis we estimate the NC background^ by assuming an energy independent background 
suppression factor of ~ 10~^ which agrees with [117]. We have noted that after five years of 
running of the CERN-INO beta-beam experiment with 7 = 650, we expect only about 0.1 
NC background events. Nevertheless we take this background into account in our numerical 
analysis. The background is assumed to have the same shape as the signal. But one should 
keep in mind the fact that this shape is not much of an issue since anyway the background is 
very small, we do not consider any backgrounds coming from the CC events of u^. 

Since the oscillation probability —>■ Vt (silver channel) is about the same as that for z/g — > 
(golden channel), we expect almost as many z/^ arriving at the detector as v^. The r~ produced 
through CC interaction may decay further inside the detector producing secondary /i~ with a 
branching ratio of 17.36%. But, the r threshold (3.5 GeV) is high and the production cross 
section suppressed compared to /i. So we do not expect any significant background from 
this source either. We have estimated the number of secondary muons produced from the Vr 
component of the beam. For sin^ 2^13 = 0.01, we expect 0.3 muon events per year with 7 = 500. 
In addition, these secondary muons will be severely degraded in energy and therefore can be 
eliminated through energy cuts. We therefore neglect the backgrounds from this source. 

-'^^ Mesons produced in NC processes are degraded in energy. Note that backgrounds from these mesons are 
very important in the case of the neutrino factory. However, since our relevant energy range is lower, the mesons 
produced in each event are much lower in energy and hence can be easily rejected by putting suitable cuts. 
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Atmospheric Backgrounds 



The atmospheric neutrino flux falls steeply with energy and is expected to produce much 
fewer events for the energy range that we are interested iij^. The fact that INO has charge 
identiflcation capabihty further reduces the atmospheric background. The most important 
handle on the reduction of this background comes from the timing information of the ion 
bunches inside the storage ring. For 5T magnetic fleld and 7 = 650 for ions, the total length 
of the storage ring turns out to be 19564 m. We have checked that with eight bunches inside 
this ring at any given time, a bunch size of about 40 ns would give an atmospheric background 
to signal ratio of about 10~^, even for a very low sin^ 2^13 of 10^^. For a smaller bunch span, 
this will go down even further. In addition, atmospheric neutrinos will be measured in INO 
during deadtime and this can also be used to subtract them out. Hence we do not include this 
negligible background here. 



3.5 Details of the Statistical Method 



In order to quantify the sensitivity of this CERN-INO beta-beam experimental set-up to the 

>^sis of the "data" generated in 
For our statistical analysis we 



mixing angle 613 and sgn{Am'^j), we perform a statistical analysis of the "data" generated in 
ICAL@INO, assuming certain true values of the parameters 
employ a function deflned as 



Xtotal 



+ x 



prior ' 



(3.24) 



where the first term is the contribution from the neutrino channel, the second term comes from 
the antineutrino channel, while the last term comes from imposing priors on the oscillation 
parameters which we allow to vary freely in our fit and which we expect will be determined 
better from other experiments at the time when the data from the CERN-INO beta-beam 
set-up would be finally available. The for the neutrino channel is given by 



Xv 



2^(?/j -Xi - Xjln 



1=1 



Xi 



+ es + a 



where n is the total number of bins. 



(3.25) 



(3.26) 



N^^{{u}) given by Eq. (13.191) being the predicted number of events in the energy bin i for a set 
of oscillation parameters uj and are the number of background events in bin i. The quantities 
TT'^ and n'' in Eq. f l3.26p are the systematical errors on signals and backgrounds respectively. 



We have taken vr* = 2.5% and 



5% (see Table 13. 4p . The quantities and are the 



"pulls" due to the systematical error on signal and background respectively. The data in Eq. 



^^We will show in the next section that even a threshold energy of 4 GeV is easily admissible in our set-up, 
and above 4 GeV there are much fewer atmospheric events. 
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13.251 enters through the variables Xi = Nf^ + N^, where Nf^ are the number of observed signal 
events in the detector and is the background, as mentioned earlier. We simulate the signal 
event spectrum using Eq. 13.191 for our assumed true values for the set of oscillation parameters 
which are given in the first column of Table 13.21 Different options are used for sin^ 26i^ (true) , 
5c'p(frue) and the true hierarchy and these are mentioned wherever applicable. In our fit we 
marginalize over all oscillation parameters, the Earth matter density, as well as the neutrino 
mass hierarchy, as applicable. We do this by allowing all of these to vary freely in the fit and 
picking the smallest value for the function. Of course, we expect better determination of 
some of these parameters, which are poorly constrained by this experimental set-up. Therefore, 
we impose a "prior" on these parameters through the x^rior given by 



where the lex error on these that we use are taken from [90, 112] and are given in the right 
column of Table 13. 2[ In our computation, we have used a matter profile inside the Earth with 
24 layers. In Eq. 13.271 p is a constant number by which the matter density of each layer 
has been scaled. The external information on p is assumed to come from the study of the 
tomography of the earth [113]. In Eq. 13. 27^ p varies from 0.95 to 1.05 i.e., 5% fluctuation 
around 1. 

Note that in our definition of the x^ function given by Eqs. 13.241 and 13.251 we have assumed 
that the neutrino and antineutrino channels are completely uncorrelated, all the energy bins 
for a given channel are fully correlated, and and are fully uncorrelated. We minimize the 
xlotai two stages. First it is minimized with respect to and to get Eq. I3.25[ and then 
with respect to the oscillation parameters uo to get the global best-fit. For minima with respect 
to and ^6, we require that 




(3.27) 



dx^ 



and 



dx^ 



. 



(3.28) 



From Eqs. [31251 13:261 K28\ we get 




(3.29) 



where. 
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Figure 3.14: The range of sin^ 26'i3(true) for which the wrong hierarchy can be ruled out at the 
3a C.L., as a function of 7. The left panel is for NH as true, while the right panel is when 
IH is true. The red solid curves show the sensitivity when the 7 is chosen to be the same for 
both the neutrino and the antineutrino beams. The blue dashed lines show the corresponding 
sensitivity when the 7 for the antineutrinos is scaled down by a factor of 1.67 with respect to 
the 7 of the neutrino beam. 



an 



^22 



ai2 




(3.30) 



i=l 



Using Eqs. 13.291 and 13.301 we calculate the values of and and then we use these values to 
calculate xl^-,u ■ ^ similar fashion, we estimate xl^-,u obtain the Xtotai- 



3.6 Measurement of the Neutrino Mass Ordering 



Let us first focus on the issue of neutrino mass ordering. If the true value of 613 has indeed been 
chosen to be large by Nature, sizable matter effects can be expected in this magical CERN- 
INO beta-beam experimental set-up, giving us a handle on the mass hierarchy. The mass 
hierarchy sensitivity is defined as the range of sin^ 2^^i3(true) for which the wrong hierarchy can 
be excluded at the 3a CL. The left (right) panel of Fig. 13.141 shows the hierarchy sensitivity 
when the NH (IH) is true. We show this as a function of 7. In both panels we show by the red 
solid lines the sensitivity when we add neutrino and antineutrino data with the same value of 
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Figure 3.15: Plots showing the impact of various factors on the mass hierarchy sensitivity of 
the CERN-INO beta-beam experiment. The top left panel shows the impact of changing the 
detector threshold. The lower left panel shows the effect of changing the background rejection 
factor. The top right panel shows the difference in the sensitivity between the rate and spectral 
analysis. The lower right panel shows how the density profile would impact the hierarchy 
sensitivity. 
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Figure 3.16: Effect of 5cp(true) on the hierarchy sensitivity. The black dashed hnes show the 
worst and best cases when we allow 5(7p(true) to take any value between and 27i. The red 
solid curve corresponds to the reference case where (5cp (true) = 0. The left panel shows the 
case for true NH while the right panel is for true IH. 

the Lorentz boost shown in the x-axis. The blue dashed lines on the other hand correspond to 
the sensitivity expected when the neutrino beam runs with the 7 shown in the x-axis while the 7 
for the antineutrino beam is scaled down by a factor of 1.67. We assume years of running 
of the beta-beam in both polarities and a full spectral analysis has been performed. We note 
that using the combined neutrino and antineutrino beam running at the same value of 7 the 
wrong hierarchy could be ruled out at 3a for sin^ 2^13 > 6.8 x 10~^ (sin^ 26'i3 > 4.0 x 10~"^) for 
7 = 500 if the NH (IH) is true. Presence of both neutrino and antineutrino data simultaneously 
in the analysis restricts the fitted value of ^13 to be in a range very close to the assumed true 
value. For instance, for NH true, data corresponds to a large number of events for neutrinos and 
a small number of events for antineutrinos. When this is fitted with IH, we have a small number 
of events predicted for the neutrinos. In order to minimize the disparity between the data and 
prediction for neutrinos, the fit tends to drive ^13 to its largest allowed value. However, larger 
values of ^^13 would give very large number of antineutrino events for IH and this would be in 
clear conflict with the data. Therefore, the net advantage of adding data from both neutrino 
and antineutrino channels is that one cannot artificially reduce the ^iny longer by tinkering 
with ^13 in the fit. As a result, the sensitivity of the experiment to mass hierarchy witnesses a 
substantial improvement. 

It can be noted from the plots that the hierarchy sensitivity falls when the scaled 7 option for 
the antineutrino beam has been used. This is particularly relevant when the true hierarchy is 
inverted and/or when 7 is low. Since scaling the 7 reduces it by a factor of 1.67, the statistics 
for the antineutrinos fall by nearly a factor of 1.67 for this case and this reflects in the reduced 
hierarchy sensitivity of the experiment. Its impact when true hierarchy is inverted is more 

^^Here all sensitivity figures correspond to a five years run. 
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Figure 3.17: The variation of the experimental sensitivity on the number of useful ion decays 
in the straight sections of the storage ring. Left panel shows sensitivity to the mass hierarchy 
assuming NH to be true. Right panel shows the sin^ 2^13 sensitivity reach. In both panels, the 
majenta solid vertical line corresponds to the reference value used in the rest of the analysis. 



because in that case, the data corresponds to larger events for the antineutrinos and very small 
events for the neutrinos. The antineutrino events are therefore the driving force and an increase 
in their statistical uncertainty due to the scaled down 7 accentuates the adverse effect on the 
hierarchy sensitivity. In the case of NH, the events in the neutrino channel are the dominant 
factor and the role of the antineutrinos is only to prevent the 6*13 values in the fit to run to very 
large values, as discussed before. As long as the antineutrino events have enough statistical 
power to restrict 6*13 to values close to the true value at which the data was generated, the 
hierarchy sensitivity remains reasonably good. Therefore, for the NH only for very low values 
of 7 the hierarchy sensitivity gets seriously affected by the Lorentz boost scaling. 

Fig. 13.151 shows how the hierarchy sensitivity depends on diverse input factors. As in Fig. 13.141 
we show the 3a limit for sin^ 29i3 as a function of 7 in all the four panels and we assume that 
NH is true. The reference curve (red solid line) in all panels corresponds to the result obtained 
with a Uf. and beam with a spectral analysis. The upper left hand panel shows the effect 
of changing the threshold energy of the detector. The sensitivity of the experiment is seen 
to remain almost stable against the variation of the threshold energy upto 4 GeV. Only for 
a threshold of 6 GeV and above the sensitivity falls, the lower 7 values getting more affected 
since they correspond to lower neutrino energies. In the lower left hand panel of the figure 
the effect of the chosen background fraction on the hierarchy sensitivity has been shown. The 
red solid line shows the sensitivity for our assumed background factor of 10~^ while the black 
dashed line shows the corresponding sensitivity when the background rejection is poorer and we 
have a higher residual background fraction of 10~^. The upper right hand panel shows how our 
sensitivity increases by taking into account the spectral information of the events. It also shows 
how much improvement we get by combining the antineutrino data with the neutrino data. The 
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black dashed line shows how the sensitivity falls when we use the total event rates instead of 
the events spectrum. The blue dashed-dotted and green dashed-triple-dotted lines show the 
sensitivity expected from the neutrino data alone. The blue dashed-dotted line is for binned 
neutrino data while the green dashed-triple-dotted lines shows the sensitivity for the total 
event rate for neutrinos alone. It can be seen that the effect of using the spectral information is 
only marginal when both neutrino and antineutrino are used, while the effect of combining the 
antineutrino data with the neutrino data on the sensitivity is huge. For the neutrino data alone, 
the sensitivity improves significantly when one uses the spectral information. In the lower right 
hand panel we show how the sensitivity of the experiment to hierarchy would get affected if we 
use a different profile for the Earth matter density instead of PREM. The red solid line is for 
earth density according to the PREM profile while the blue dotted and black dashed lines are 
when the matter density is 5% lower and 5% higher respectively than the density predicted by 
the PREM profile. When the density is higher (lower) the matter effects are higher (lower) and 
therefore the sensitivity improves (deteriorates). 

One crucial point that we have not stressed so far concerns the dependence of the detector 
performance on the true value of 6c p- AH the earlier plots were presented assuming that 
5(7p(true) = 0. At exactly the magic baseline, it is expected that the sensitivity of the exper- 
iment to be completely independent of 6cp- The CERN-INO distance of 7152 km is almost 
magical, but it is not the exact magic baseline. Therefore, we do expect some remnant impact 
of Scpit^vie) on our result^. To show how our results get affected by Scp{tTue), we show in 
Fig. 13.161 the hierarchy sensitivity just as in Fig. 13.141 but here we show the full band corre- 
sponding to all values of (5cp(true) from to 27r. As before, the left panel is for NH true while 
the the right panel is for IH true, and we have taken in the analysis the full spectral data for the 
neutrinos as well as the antineutrinos, with the same 7. The lower edge of this band shows the 
best possible scenario where the experiment is most sensitive, while the upper edge shows the 
worst possible sensitivity. The red solid lines in both panels show for comparison the hierarchy 
sensitivity corresponding to 5cp(true) = 0, which we had presented in Fig. I3.14[ We note from 
the figure that the hierarchy sensitivity is nearly the best for 6cp{tTne) = when IH is true 
while if NH is true then it would give us almost the worst sensitivity. For NH (IH) as true the 
best possible sensitivity would be sin^ 2^13 > 3.96 x 10"^ (sin^ 2^13 > 2.96 x 10^^) for 7 = 650 
to be compared with sin^ 2^13 > 5.51 x 10"'^ (sin^ 26113 > 3.05 x lO"'*) when 5cp(true) = 0. 
Therefore, we conclude that if NH is true then it would not be unfair to expect an even better 
hierarchy sensitivity than what was reported in Fig. 13.141 while if IH is true then the best 
sensitivity will be returned for 6cp{tTue) ~ 0. 

It has been noted from Figs. 13.111 and 13.121 that the total number of events in the detector 
increases roughly linearly with 7, except for extremely long baselines. Increasing the number 
of ion decays per year will also bring about a simple linear increase in the statistics. It is 
therefore pertinent to make a fair comparison between the dependence of the mass hierarchy 
sensitivity to the Lorentz boost 7 and the number of useful ion decays in the rin^. In the left 

"'^'^Note that in all our results presented in this work, we have fully marginalized over all the oscillation 
parameters in the fit, including Sep- 

-'^^Note that this is also equivalent to increasing the total exposure time of the experiment. Both number of 
ion decays per year and exposure appear as a normalization factor for the event rate and hence increasing the 
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7\N 


Mass Hierarchy {3a) 


sin^ 26'i3 sensitivity {3a) 


1.1 X W 


2.043 X 10^« 


1.1 X W 


2.043 X 10^« 


350 
650 


1.3 X 10-3 
5.6 X 10"^ 


9.3 X 10"^ 
4.1 X 10"^ 


3.8 X 10-3 
1.1 X 10-3 


2.3 X 10-3 
7.3 X 10-^ 



Table 3.5: Comparison of the variation of the detector sensitivity to mass hierarchy (columns 
2 and 3) and sin^ 2^13 sensitivity (columns 4 and 5) with 7 and N, the number of useful ion 
decays per year. 



panel of Fig. 13.171 we show the effect of increasing the number of ion decays on the hierarchy 
sensitivitjcj. The plots exhibit the dependence of the sensitivity on the number of useful ion 
decays per year for an exposure of five years, for three different values of 7. The same Lorentz 
boost for the neutrino and antineutrino beams have been assumed. We present in Table 13.51 
the relative increase in the hierarchy sensitivity when we increase the 7 by a factor of 1.86 and 
compare it against the increase in the sensitivity when the number of ion decays are increased 
by the same factor. It can be noted that while the hierarchy sensitivity improves by a factor of 
2.54 in going from 7 = 350 to 650 keeping the number of ion decays per year as 1.1 x 10^^, it 
increases 1.5-fold when we raise the number of ion decays per year from 1.1 x 10^*^ to 2.04 x 10^^ 
keeping 7 = 350. However, we would like to stress that the improvement of the hierarchy 
sensitivity is not linear with either 7 or number of ion decays per year. The crucial thing is 
that the behavior of the sensitivity dependence on both 7 and number of ion decays per year 
is very similar. It increases very fast initially and then comparatively flattens out. 



3.7 



Measurement of sin^ 29 



13 



The CERN-INO beta-beam set-up is also expected to give very good sensitivity to the 6*13 
measurement. In what follows, we will quantify our results in terms of three "performance 
indicators" , 



sin 2^13 sensitivity reach, 



2. 



sin^ 2^^i3 discovery reach, 



sin 26'i3 precision. 



number of ion decays by a factor n keeping the exposure same is equivalent to increasing the exposure by a 
factor n keeping the number of ion decays per year fixed. 

^^We assume that the number of useful ion decays for both and ^Li have been scaled by the same factor. 
In the figure along the a;-axis only the numbers are shown. 
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Figure 3.18: Left panel shows the 3a sensitivity hmit for sin^26'i3. Right panel shows the 3a 
discovery reach for sin^ 2^13 (true). The red solid lines in the left and right panels show the 
sensitivity reach and discovery potential respectively, when the 7 is assumed to be the same for 
both the neutrino and the antineutrino beams. The blue dashed lines show the corresponding 
limits when the 7 for the ^Li is scaled down by a factor of 1.67 with respect to the 7 of the 
neutrino beam, which is plotted in the x-axis. 

Below a detailed description of our definitions of these performance indicators has been given. 
All results in this section have been obtained by taking into account the full event spectrum 
and combining five years data from both the neutrino and antineutrino channels. 

3.7.1 sin^2^i3 Sensitivity Reach 

The 9i3 sensitivity reach is defined as the range of sin^ 26*13 which is incompatible with the 
data generated for sin^ 26*13 (true) = at the 3a CL. This performance indicator corresponds 
to the new sin^ 2^13 limit if the experiment does not see a signal for 6'13-driven oscillation^. 
In that case, we can exclude some allowed values of sin^ 26*13, which we call our "^13 sensitivity 
reach". We simulate this situation in our analysis by generating the data at sin^ 2^13 (true) = 
and fitting it with some non-zero value of sin^ 2^13 by means of the technique. In our fit 
we marginalize over all the oscillation parameters including Sep and the mass hierarchy!^ and 
choose the value of sin^ 2^13 for which the fit yields = 9. The result is shown in the left 
panel of Fig. 13.181 as a function of 7. The red solid line shows the sin^ 26*13 sensitivity when 7 

^^Note from Eq. 13.41 while the first three terms go to zero when ^13 0, the last term, which depends 
only on the solar parameters and 023, remains non- vanishing. Therefore, when the flux is high, i.e., for large 7 
and/or enhanced luminosity, we expect a sizable number of events even when sin^ 2^13 (true) — 0. 

^^Note that since sin^ 20i3(true) = 0, the data is independent of the Scp{true) and the true mass hierarchy. 
However, since we allow for non-zero sin^ 2f?i3 in the fit, the predicted event rates in our "theory" depend on 
Sep and the mass hierarchy. 
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is assumed to be the same for both the neutrino and the antineutrino beams. The blue dashed 
hues shows the corresponding Sa upper hmit when 7 for the ^Li is scaled down by a factor 
of 1.67 with respect to that for the neutrino beam. In generating the data we have assumed 
that sin^ 2^13 (true) = 0, which means that we have negligible events in both the neutrino as 
well as the antineutrino channels, irrespective of the mass hierarchy. When this data is fitted 
allowing for non-zero sin^2^i3, the neutrino (antineutrino) channel plays a dominating role 
when NH (IH) is assumed in the fit. We reiterate that Fig. 13.181 shows the sin^26'i3 sensitivity 
after marginalizing over hierarchy as well. In other words, the sensitivity shown in this figure 
corresponds to the statistically weaker channel. For the case where the same 7 for and 
^Li has been used, the neutrino channel is weaker since the event rate is about 1.5 times less 
than antineutrino events with the same 7. On the other hand when we scale down the Lorentz 
boost for ^Li, the flux in the antineutrino channel goes down significantly and hence it becomes 
the statistically weaker channel as can be seen from Fig. 13.131 and therefore the marginalized 
corresponds mainly to that from antineutrinos. Indeed one can check that the sin^ 26'i3 
sensitivity that we exhibit by the blue dashed line for the scaled 7 case is comparable to what 
one can obtain for the antineutrino channel with IH and the corresponding lower 7. Similar 
feature is observed in addressing the issue of neutrino mass ordering taking IH as true hierarchy 
(see right panel of Fig. 13.141) and the 7 for the antineutrinos is scaled down by a factor of 1.67 
with respect to the 7 of the neutrino beam. 

The dependence of the sin^ 2^13 sensitivity on the number of useful radioactive ion decays per 
year in the straight section of the storage ring is shown in the right panel of Fig. 13.171 Here 
we have taken the same Lorentz boost for and ^Li and we have shown the results for three 
fixed values of 7. The relative increase in the sensitivity by increasing 7 and/or the number of 
useful ion decays per year by the same factor is quantified in the last two columns of Table 13. 5[ 

3.7.2 sin^2^i3 Discovery Reach 

How good are our chances of observing a positive signal for oscillations and hence ^13 in the 
CERN-INO beta-beam set-up? We answer this question in terms of the parameter indicator 
which we call the "discovery reach" of the experiment for sin^ 26'i3. This performance indicator 
is defined as the range of sin^ 26'i3(true) values which allow us to rule out sin^26'i3 = at the 
3(7 C.L. To find this we simulate the data at some non-zero value of sin^ 2^13 (true) and fit it 
by assuming that sin^ 2^13 = 0, allowing all other oscillations parameters to take any possible 
value in order to return back the smallest value for the . Note that since the fitted value 
of ^13 in this case always corresponds to 0, there is no need of any marginalizing over the 
hierarchy when fitting the data. However, since the data here is generated at a non-zero value 
of sin^ 2^13 (true), it depends on the true mass hierarchy. The discovery reach of the experiment 
is therefore expected to be dependent on the true mass hierarchy. Likewise, while the value 
of bcp in the fit is inconsequential as sin^25i3 = in the fit, the data itself and hence the 
discovery reach, would depend on ^c-p (true). For each sin^ 26^13 (true), we generate the data for 
all possible values of (5(7p(true) and for both the true mass hierarchies. For each case, the data 
is then fitted assuming sin^ 26'i3 = and marginalizing over the other oscillation parameters. 
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Figure 3.19: The precision with which sin^ 26*13 will be measured by the CERN-INO beta-beam 
experiment as a function of sin^ 26*13 (true). Left panel shows the 3a allowed range of sin^ 26*13 
while the right panel shows the precision defined in the text. 



returning a value of Xmin each data set. We choose the minimum amongst these xLin find 
the value of sin^ 2^13 (true) for which we could claim a signal in the detector at the 3a CL. In 
the right panel of Fig. 13.181 we show this "most conservative"!^ sin^ 26*13 discovery reach of our 
experiment as a function of 7. We assume equal 7 for both the ions for the red solid curve. One 
can see that for 7 = 650, the most conservative discovery reach is sin^ 26*13 (true) = 5.11 x 10^^ 
while if 5cp(true) = then we have checked that the reach is slightly better and this will be 
sin^ 2^13 (true) = 5.05 x 10~^. For the blue dashed line we assume that the 7 for ®Li is scaled 
down by a factor of 1.67 compared to that for ^B, plotted on the x-axis. Since for same 7, 
neutrino is the statistically weaker channel, the red line mainly corresponds to what we expect 
for the true NH. For the scaled 7 case since the antineutrino channel becomes statistically 
weaker, the lower comes from this channel and the blue dashed line corresponds to what we 
expect for the true IH. 



3.7.3 sin^ 29is Precision 



In Fig. 13.191 we show how precisely the mixing angle sin^ 2^13 will be measured, if we observe 
a 6*13 driven signal at the detector. The left panel depicts as a function of sin^ 26*13 (true) the 
corresponding range of allowed values of sin^ 2^13 at the 3a CL. We have assumed 7 = 500 and 
5c'p(true) = 0. The solid line is assuming NH to be true, while the dashed line is for IH true. 
Note that in the fit we always marginalize over the hierarchy and 6c p- The right panel shows 

^^This is "most conservative" in the sense that no matter what the choices of (5cp(true) and the true neutrino 
mass ordering, the ^13 discovery hmit cannot be worse than the value presented. 
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the variable "precision" which we define as 

precision = j^^^' ^^^3^,. - (sin^ '29,,) ^ 

where (sin^ 2^13)^^^ and (sin^ 2^13)^^^ are the maximum and minimum allowed values of 
sin^ 2^13 respectively at Scr. 



3.8 Summary and Conclusions 



Long baseline experiments which will use the golden channel for determining the neutrino 
oscillation parameters face a serious threat from the menace of clone solutions due to the so- 
called parameter degeneracies. These degeneracies come in three forms: the 5c p — ^13 intrinsic 
degeneracy, the Sep — sgn{A'mli) degeneracy and the 623 octant degeneracy, and necessarily 
result in degrading the sensitivity of the experiment. The CERN-INO near-magic distance of 
7152 km offers the possibility of setting up an experiment at a baseline where the 6c p dependent 
terms almost drop out from the expression of the golden channel probability. Thus two out 
of the three degeneracies are evaded, providing a platform for clean measurement of ^13 and 
sgn{Am^i), two major players in our understanding of the origin of neutrino masses and mixing. 
A large magnetized iron calorimeter with a total mass of at least 50 kton is expected to be 
built soon at INO. It will be ideal for detecting multi-GeV z/^ and hence can be used as the far 
detector for a high energy beta-beam. 

In this chapter we studied in detail the physics reach of the CERN-INO magical beta-beam 
set-up. For 7 = 650, 6cp{tTue) = and true NH, the sensitivity to hierarchy determination 
at 3cr is sin^ 2^13 (true) = 5.51 x 10~^ when full spectral data from neutrino and antineutrino 
channels are combined. Even though the effect of 5cp(true) on the event rate of such an 
experiment is expected to be small, there is some residual dependence on it because the CERN- 
INO distance does not conform to the exact magic baseline. We studied the change in the 
hierarchy sensitivity due to the uncertainty in 6cp- It turns out that for 7 = 650 and with 
NH (IH) true, the best sensitivity to hierarchy determination corresponds to sin'^ 2^^i3(true) = 
3.96 X 10-'' (sin^ 26'i3(true) = 2.96 x 10"^), while the worst case is sin^ 26'i3(true) = 5.58 x 10"^ 
(sin^ 2^13 (true) = 3.59 x 10"^). 

A detailed analysis of the potential of probing ^13 at this experiment has been presented here. 
We defined and studied the ^13 reach in terms of three performance indicators: the sensitivity 
reach, the discovery reach and the precision of sin^2^i3 measurement. The ^13 sensitivity 
reach is defined as the range of sin^ 2^13 which is incompatible with the data generated for 
sin^ 26'i3(true) = at the 3a C.L. The sensitivity reach corresponds to sin^ 26'i3 = 1-14 x 10"^ 
at 3(7 C.L. for 7 = 650 and this is independent of the true hierarchy and 6cp{tTue). The 
discovery reach is defined as the range of true values of the mixing angle for which we have an 
unambiguous oscillation signal in the detector. At 3<t C.L. the discovery reach corresponds to 
sin^ 2^13 (true) = 5.05 x lO^^ (sin^ 2^13 (true) = 2.96 x 10"^) for 7 = 650, 6cp{truc) = and 
NH (IH) true while the most conservative limit irrespective of 5cp(true) and the true neutrino 
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mass ordering is sin^ 2^^1,3 (true) = 5.11 x 10 We also presented the expected precision with 
which sin^26'i3 would be determined in this experiment for sin^ 26'i3(true) > 10~^. 

Neutrino physics is in wait for the next great leap forward in the decade ahead. The study of 
the golden channel probability Pe^ using a beta-beam neutrino source and an iron calorimeter 
detector at a very long magical baseline may well turn out, as we have demonstrated in this 
chapter, to be very crucial in this endeavour. In the next chapter, we will discuss the impact 
of matter effects in the i/g survival probability at long baselines in extracting the information 
about the neutrino mass ordering and the 1-3 mixing angle. 
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Chapter 4 



Neutrino parameters from matter 
effects in Pge sit long baselines 



4.1 Introduction 



Determination of the 1-3 neutrino mixing angle ^13, sgR{Amli), the three CP phases and the 
absolute neutrino mass scale are necessary for reconstruction of the neutrino mass matrix, which 
will have important consequences for nuclear and particle physics, astrophysics and cosmology. 
The golden channel Fg^t has the capability of measuring the Dirac phase Scp, sgnlAm^-i^) and 
^13 in long baseline accelerator based experiments. However, this strength of the golden channel 
also brings in the well-known problem of parameter "degeneracies" , where one gets multiple 
fake solutions in addition to the true one [93-96]. Various ways to combat this vexing issue 
have been suggested in the literature, including combining the golden channel with the "silver" 
(Per) [106] and "platinum" (P^e) channels. While each of them would have fake solutions, their 
combination helps in beating the degeneracies since each channel depends differently on 6c p, 
sgn{Am'^^) and ^13. In this chapter, we focus on the Ve survival channel. Pee, which is 

independent oi Sep and the mixing angle 623- It is therefore completely absolved of degeneracies 
and hence provides a clean laboratory for the measurement of sgn(Am'li) and ^13. This gives 
it an edge over the conversion channels, which are infested with degenerate solutions. 

The Pee survival channel has been extensively considered for measuring ^13 with Pg produced 
in nuclear reactors [121] and with detectors placed at a distance ~ 1 km. Reducing systematic 
uncertainties to the sub-percent level is a prerequisite for this program and enormous R&D is 
underway for this extremely challenging job. For accelerator based experiments, the survival 
channel. Pee, has been discussed with sub-GeV neutrinos from a beta-beam source at CERN 
and a megaton water detector in Prejus at a basehne of 130 km [62,63,65, 107]. However, no 
significant improvement on the ^13 limit was found in [107] for a systematic error of ^ 5%. 
This stems mainly from the fact that in these experiments one is trying to differentiate between 
two scenarios, both of which predict a large number of events, differing from each other by a 
small number due to the small value of ^13. Also, since sgn{Am'^i), is ascertained using earth 
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matter effects, there is no hierarchy sensitivity in these survival channel experiments due to the 
short baselines involved. 

In this chapter, we emphasize on the existence of large matter effects in the survival channel, 
Peei for an experiment with a very long baseline. Recalling that Pee = 1 — Peii ~ Per and since 
for a given sgnlArn^^) both Pefj, and Per will either increase or decrease in matter, the change 
in Pee is almost twice that in either of these channels. Using the multi-GeV z/g flux from a 
beta-beam source, we show that this large matter effect allows for significant, even maximal, 
deviation of Pee from unity. This, can thus be a convenient tool to explore ^13. This is in 
contrast to the reactor option or the beta-beam experimental set-up in [107], where increasing 
the neutrino flux and reducing the systematic uncertainties are the only ways of getting any 
improvement on the current ^13 limit. We further show, for the first time, that very good 
sensitivity to the neutrino mass ordering can also be achieved in the Pee survival channel owing 
to the large matter effects. We discuss plausible experimental set-ups with the survival channel 
and show how the large matter effect propels this channel, transforming it into a very useful 
tool to probe sgn{Am1i) and di^ even with relatively large room for systematic uncertainties. 



4.2 The z/g i^e Survival Probability 



Under one mass scale dominance and the constant density approximation, the matter conversion 
probabihtes are, 

Pe^ = 1^23 sin^ 2^™ sin^ [l.27(Am2^)"L/E] , (4.1) 
where 1^23 = sin^ ^23 for x = /i and 1^23 = cos^ ^23 for x = r. Then we have 

Pee = 1 - Pe^, - Per = I - siu^ 2^™ slu^ [l.27( Am^^ . (4.2) 

The largest deviation of Pee from unity is obtained when the conditions 13.141 and 13.151 are 
satisfied simultaneously. This ensures maximum matter effect in the conversion channels (see 
Eqs. 13.171 and 13.181) . However, there are suppression factors, sin^ 6*23 for Pe^ and cos^ 6*23 for 
Per, not present in Pge. Moreover, since Pee does not contain 6*23, the octant ambiguity as well 
as parameter correlations due to uncertainty in ^^23 [95] are absent. In addition, as mentioned 
earlier, the Pee channel does not contain the CP phase, 5cp and therefore this channel does not 
suffer from the effect of the so-called (6*13, 5cp) intrinsic degeneracy [93] and the (sgn(Am|]^), 
5cp) degeneracy [94]. Both of these remain true in the presence of non-zero Am^i [122]. 

Using Eq. I3.18[ one can estimate the distance where Pee — for a given value of ^13. For instance 
for p = and sin^ 26'i3 = 0.2 and 0.1, these distances are 7600 km and 10200 km respectively 
which can be readily seen from Fig. 13.61 For higher values of p the distance exceeds the earth's 
diameter for ^13 in the current allowed range. Using [pL)^"-^ corresponding to the PREM 
profile, from Eq. 13.181 one can estimate that the condition of maximal matter effects inside the 
earth's mantle is satisfied only for sin^ 26'i3 ^ 0.09. In our numerical work, we solve the full 
three flavour neutrino propagation equation assuming the PREM [110] profile and keep Am^^ 
and sin^ 612 fixed at their present best-fit values of 8.0 x 10~^ eV^ and 0.31 respectively [78] as 
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Figure 4.1: The survival probability Pee in matter as a function of E for four different values 
of the baseline L. For each L, the plots are given for three different values of sin^2^i3 (0.17, 
0.1 and 0.05). Thick (thin) lines are for NH (IH). 



mentioned in the first column of Table [3^ For the other oscillation parameters, the benchmark 
values are also taken from the same table. 

In Fig. I4.1l we plot Pee as a function of energy, at four different L and for three values of sin^ 26*13. 
The plots confirm that maximal matter effects come at L ^ 10000 km and L ~ 7500 km for 
sin^ 26*13 = 0.1 and 0.17 respectively for the normal hierarchy (NH). For the inverted hierarchy 
(IH) there is no significant matter effect for u^. This large difference in the probabilities for NH 
and IH can be exploited for the determination of sgn{Am1i). Further, since the matter effect 
is a sensitive function of ^13 it may also be possible to obtain information on this angle. We can 
also see that for a given value of sin 2^13 ( > 0.09) and E, the matter effect increases (almost 
linearly) with L, until the L for maximal matter effect is reached, beyond which matter effect 
falls. For values of sin^ 2^13 < 0.09 the condition for maximum matter effect is not met inside 
the earth's mantle and hence the matter effect and sensitivity to both hierarchy as well as ^13 
increase with L. 

In what follows, we will show how, in a plausible experiment, one can use this near- resonant 
matter effect in the survival channel. Pee, to constrain 6*13 and sgn{Amli). Fig. 14. II shows that 
the requirements for such a program include 



1. 



a i/g beam, 
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2. 



a baseline of at least a few thousand km, 



3. average energies around 6 GeV, 



a detector capable of observing e unambiguously at these energies. 



4.3 The Experimental Set-up 



4.3.1 Pure z/g {i^e) Source 

The various issues related to beta-beams have already been discussed in great detail in section 
3.2 of the previous chapter. The ions considered as possible sources for beta-beams are ^^Ne 
and for and ^He and ^Li for z/g. The end point energies of ^He and ^^Ne are ~ 3.5 MeV 
while for and ^Li this can be larger ~ 13-14 MeV (see Table [LT!) . For the Lorentz boost 
factor 7 = 250 (500) the ^B and ^Li sources have peak energy around ~ 4 (7) GeV. Since this 
is in the ball-park of the energy necessary for near-resonant matter effects as discussed above, 
we will work with ®B (^Li) as the source ion for the {i^e) beta-beam and 7 = 250 and 500. 



4.3.2 Water Cerenkov Detector 

Water Cerenkov detectors have excellent capability of separating electron from muon events. 
Since this technology is very well known, megaton water detectors are considered to be ideal 
for observing beta-beams. Proposals for megaton water detectors include UNO [123] in USA, 
HyperKamiokande [124] in Japan and MEMPHYS [125] in Europe. If the beta-beam is pro- 
duced at CERN, then baselines in the range 7000-8600 km would be possible at any of the 
proposed locations for the UNO detector. Likewise, if the beta-beam source be at FNAL, then 
the far detector MEMPHYS would allow for L = 7313 km. HyperKamiokande could also be 
considered as the far detector and in that case L = 10184 (9647) km if the source be at FNAL 
(CERN). Such detectors do not have any charge identification capacity. But in a beta-beam, 
the f3~ and emitters can be stacked in different bunches and the timing information at the 
detector can help to identify the e~ and e+ events [60]. 



4.3.3 Backgrounds 

It is well known that there are no beam induced backgrounds for beta-beams. In the previous 
chapter (subsection 3.4.3), these issues have been addressed in detail. In this experimental 
set-up, the process z/g — i^,- — ^ t~ — > e~ could mimic the signal. We have checked that the 
background to signal ratio for these events in the relevant energy range is ~ 10^^ and can be 
neglected for the disappearance mode. The electron events from K and tt~ decays are also 
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Figure 4.2: Events in 5 years vs. sin^ 2^13 for NH (dashed line) and IH (solid line) for L = 7500 
km and 7 = 500. The inset shows the same but for 7 = 250. 



negligible. The atmospheric background can be estimated in the beam off mode and reduced 
through directional, timing, and energy cuts. 



4.4 Simulation Details and Event Rates 



For our numerical analysis we use the standard technique with 



wher^ 



X.total X.Ue—>Ue ' ^priory 

(4.3) 



2 _ /|Am|i| - |Am|i(true)|\ 
^f"''" ~ V a(|Amii|) J ' ^ • ^ 

Xue^ue calculated using the same expression that we have used for xl^^u^ (^^e Eqs. 13.251 & 
I3.26P in the previous chapter. Here all the results have been presented using the data from only 
neutrino (^B) run without any spectral information. We have not considered any backgrounco 
in our analysis. We have taken 3% systematical error on signals. The prospective "data" is 
generated at the "true" values of oscillation parameters as given in the Table 13. 2[ assuming 440 
kton of fiducial volume for the detector with 90% detector efficiency, threshold of 4 GeV and 



^No prior is taken on 623 because Pee is independent of ^23- Here we have kept Atoji a-nd sin^ 6'i2 fixed 
at their present best-fit values in the fit. No uncertainty has been taken in the PREM profile when we fit the 
simulated data. 

^Backgrounds are not important here because we are dealing with the survival channel {Pee ) and the number 
of events is also quite large. 



72 



7 
6 

4 
3 
2 
1 



5 years neutrino run 
y =250 



7500 km> 




10,000 km/ 



0.003 



0.01 



0.03 



0.1 0.2 



sin29i3 (true) 
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7 = 250, as a function of sin^ 2^13 (true). 



energy smearing of width 15%. For the beta-beam we have assumed 1.1 x 10^^ useful 
decays per year and show results for 5 years of running of this beam. The number of events 
as a function of sin^ 2^13 at L = 7500 km with a 7 = 500 Ue beta-beam is shown in Fig. 14.21 
for NH and IH. The inset in Fig. 14.21 shows the number of events in 5 years expected from a 
lower 7 = 250. We have used the neutrino- nucleon interaction cross-sections from the Globes 
package [118] which are taken from [119,120]. 



4.5 



Sensitivity to sgn(Am3i) and ^13 



In Fig. 14. 3[ we show the sensitivity {na, n = yx^) of the survival channel to the neutrino 
mass ordering for L = 7500 and 10000 km and 7 = 250. If the true value of sin^ 2^13 = 0.05, 
then one can rule out the IH at the A.8a {5.0a) C.L. with L = 7500 (10000) km. For L = 7500 
(10000) km, the wrong IH can be disfavored at the 90% C.L. if the true value of sin^ 2^^i3 > 
0.03 (0.025). The sensitivity improves significantly if we use 7 = 500 instead of 250, since (i) 
the flux at the detector increases, and (ii) the flux peaks at E closer to 6 GeV, where we expect 
largest matter effects. For 7 = 500 and L = 7500 (10000) km, the IH can be disfavored at 2. 60" 
(3.8(j) C.L. for a lower value of sin^ 26'i3(true) = 0.015. The range of sin^ 2^13 (true) for which 
the IH can be ruled out at 90% and 3a C.L. for different values of L are shown in the upper 
panel of Fig. 14.41 for 7 = 250 and 500. From the figure one can see that for 7 = 500 and L = 
7500 (10000) km the wrong IH can be disfavored at the 90% C.L. if sin^ 2^i3(true) > 1.0 x 10"^ 
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Figure 4.4: The upper panel shows the range of sin^ 26'i3(true) for which the wrong IH can be 
excluded at 90% and 3a CL. while the lower panel gives the sensitivity to sin^2^i3 at various 
baselines at 90% and 3a C.L., for two values of 7. 



(8.0 X 10-3). If instead we use a total systematic error of 5% then we get the above sensitivity 
limits as sin^ 2^13 (true) > 1.6 x IQ-^ (1.2 x IQ-^) at 90% C.L. 

If the sin^ 2^13 (true) turns out to be smaller than the sensitivity reach shown in the upper 
panel of Fig. 14.41 for a given L, then it would not be possible to determine the hierarchy at 
the given C.L. However, we would still be able to put better constraints on sin^ 2^13 itself. The 
lower panel of Fig. 14.41 demonstrates as a function of L, the sensitivity to 613, i.e., the range of 
sin^ 26'i3 which is incompatible with the data generated for sin^ 26'i3(true) = at 90% and 3a 
C.L. Both Figs. 14.31 and 14.41 show that the sensitivity improves with L, even though the flux 
falls as This results from matter effects increasing with L, as noted before. For L = 7500 

(10000) km, we can constrain sin^ 2^13 < 6.3 x lO'^ (4.3 x lO'^) at the 90% C.L. for 7 =500. 
For a 5% systematic error the above numbers are changed to sin^ 2^13 < 1.0 x 10~^ (7.3 x lO"'^). 

How does this compare with alternate possibilities? If the energy can be reconstructed accu- 
rately, then the result can be improved further. For instance, for L=7500 km, if one could 
preferentially select the energy in the range 5 to 7.5 GeV, then the wrong IH can be excluded 
from true NH for sin^ 2^i3(true) = 7.47 x IQ-^ at 90% C.L. for 7 = 500. 
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We have presented our results using a z/g beam and assuming NH to be the true hierarchy. 
Similar results can also be obtained with a Vf. beam for IH. It is also possible to run both beams 
simultaneously. 



4.6 Discussions and Conclusions 



In conclusion, we propound the possibility of using large matter effects in the survival channel, 
Pee; cit loug bascliues for determination of the neutrino mass ordering (s(7n(Am3^)) and the 
yet unknown leptonic mixing angle ^13. Matter effects in the transition probabilites Pg/* and 
Per act in consonance to give an almost two-fold effect in the survival channel. In addition, 
the problem of spurious solutions due to the leptonic CP phase and the atmospheric mixing 
angle ^23 does not crop up. The development of beta-beams as sources of pure v^j beams 
enables one to exclusively study the Pge survival probability and adds a new direction to the 
prospects of a future beta-beam. In the next chapter, we will study the performance of the 
magical CERN-INO beta-beam set-up in exploring the signals of new physics which are present 
in the R-parity violating supersymmetric models. 



75 



Chapter 5 

Can R-parity violating supersymmetry 
be seen in long baseline Beta-beam 
experiments? 



5.1 Introduction 



Long baseline neutrino oscillation experiments using beta-beam hold promise of refining our 
knowledge on the third mixing angle ^^13, sgn(Am3]^) and the CP phase 5cp^ vital missing 
parameters of the neutrino mass matrix. The phenomenon of neutrino oscillation firmly estab- 
lishes the evidence for physics beyond the Standard Model and the sub dominant effects due to 
possible new physics can leave their imprints in the future long baseline beta-beam experiments 
which are supposed to give us precision below 1%. In this chapter, we will try to address the 
following questions : 

1. Can the proposed magical CERN-INO long baseline beta-beam neutrino ex- 
periment probe non-standard interactions (NSI) that are present in RPVSM? 

2. Can these NSI become fatal in attempts to further sheirpen our understanding 
of the neutrino properties? 

3. Can new physics leave its imprints at the propagation stage of neutrinos? 

Interaction of neutrinos with matter affect long baseline experiments and this becomes more 
prominent at higher values of ^^13. Various authors [108, 116, 126] have considered this effect 
for atmospheric neutrinos. Apart from the electroweak effects, there may well be NSI leading 
to flavour diagonal and flavour changing neutral currents. Here we have in mind interactions 
with quarks and leptons involving an initial and a final neutrino. If there is no change in the 
neutrino fiavour - as, for example, in exchange - this is classified as an FDNC process, while 
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it would be FCNC otherwise. RPVSM [80], which have such interactions aheady built-iij^, will 
be the main focus of our work. Recently a model in which couplings associated with FCNC and 
FDNC can be quite a bit higher than permitted in RPVSM has also been considered [127,128]. 
Naturally, here the matter effect will be further enhanced. However, as RPVSM is a well- 
studied, renormalizable model which can satisfy all phenomenological constraints currently 
available, we shall restrict our main analysis only to it and shall make qualitative remarks 
about the other model, for which our results can be easily extended. 

Consequences of FCNC and FDNC for solar and atmospheric neutrinos [10, 129, 130], and 
neutrino factory experiments [131] have been looked into. Our focus is on long baseline beta- 
beam experiments. Our analysis encompasses both NH and IH and we also incorporate all 
relevant trilinear R-parity violating couplings leading to FCNC and FDNC. Huber et al [132] 
have a somewhat similar analysis using neutrino beams obtained from muon decays. 

The very long baseline from CERN to INO will capture a significant matter effect and offers 
a scope to signal NSI. We examine whether the presence of R-parity violating (^) interactions 
will come in the way of constraining the mixing angle 6*13 or unraveling the neutrino mass 
ordering. The possibility to obtain bounds on some ^ couplings is also probed. 



5.2 ^ Supersymmetry 



In supersymmetric theories [80], gauge invariance does not imply baryon number [B) and 
lepton number (L) conservation. In the minimal supersymmetric Standard Model, L and B 
conservation is ensured by invoking 'i?-parity' (defined as = (—1) 3-8+^+25 -^j-^g^-g g jg ^j-^g 
spin). It is a discrete Z2 symmetry under which the Standard Model particles are even and 
their superpartners are odd. The imposition of such a symmetry, while it serves a purpose, is 
rather ad hoc. In general, from the naive theoretical point of view it is expected that L and B 
conservation does not hold in supersymmetric theories. However, as this leads to a very fast 
proton decay, we follow a common practice and assume that B is conserved. This can be ensured 
by replacing the Z2 symmetry of -R-parity by a symmetry, the so-called 'baryon triality' [133]. 
In such a scenario, in addition to the usual Yukawa interactions, the superpotential contains 
renormalizable L-violating trilinear A- and A'-type couplings and bilinear /ij couplings : 

^f' = Y^ l^ijkLiLjEl + X'.^^LiQ^Dl + ^iL,H^, (5.1) 

where j, /c = 1, 2, 3 are generation indices and colour and SU{2) indices are suppressed. Here 
Li and Qi are S'f/( 2) -doublet lepton and quark superfields respectively; Ei, Di denote the 
right-handed S'[/(2)-singlet charged lepton and down-type quark superfields respectively; Hu is 
the Higgs superfield which gives masses to up-type quarks. Particularly, Ajj/. is antisymmetric 
under the interchange of the first two generation indices. The bilinear couplings, /ij, are severely 



^ e.g. through squark (A'-type couphngs) or slepton (A-type couphngs) exchange. 
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constrained by the small neutrino masses. So we will discuss the phenomenology of A- and A'- 
type couplings only. Then, expanding the above superpotential in standard four- component 
Dirac notation, we have for A'-type couplings 

- eldVL - «i 4^1 - X] + h.c, (5.2) 

where only the first two terms and their hermitian conjugates are relevant for the quark-neutrino 
interactions inside the matter via squark exchange. Above, the sfermion fields are characterized 
by the tilde sign. For A-type couplings, one can write 

= ^A,,, [ ei e\vl + (e~^)*(z^l)'=ei + z>l4ei - (z - j)] + h.c, (5.3) 

where only the first two terms with i ^ j and their hermitian conjugates are responsible for the 
interactions of neutrinos with the charged leptons inside the matter via slepton exchange. The 
interaction terms given by Eq. 15.21 and [ 5. 31 violate L, as well as lepton flavour number. Suitable 
combinations of two such terms can lead to processes which are lepton flavour violating but 
L-conserving. In what follows, all these new physics couplings are assumed to be real but will 
entertain both positive and negative values. The interactions of neutrinos with electrons and 
d-quarks in matter induce transitions (i) Ui + d uj + d and (ii) Ui + e ^ uj + e. (i) is possible 
through A' couplings via squark exchange for all i,j and through Z exchange for i = j while 
(ii) can proceed via W and Z exchange for i = j, as well as through A couplings via slepton 
exchange for all 



5.3 Golden Channel Oscillations including NSI 



The golden channel probability Pg^ in the presence of standard Earth matter effects has been 
discussed in great detail in section 3.3. Now let us see what are the changes that we would 
have in the presence of NSI inside the Earth matter. In the mass basis of neutrinos 

= diag(m2, mj, ml) = U^M+M^U, (5.4) 

where M^, is the neutrino mass matrix in the flavour basis and mi, m2, and correspond to 
masses of three neutrinos. U is the mixing matrix defined by Eq. 12.181 Eq. 12.281 depicts the 
time evolution of flavour eigenstates (in the presence of three neutrinos) in matter which can 
be rewritten in the form 

(5.5) 



where 







( 















H = U[%\v^^rJJL. (5.6) 
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^ is the effective mass squared matrix, i? is a 3 x 3 matrix reffecting the matter effect in the 
form 

Rij = i?ij(Standard Model) + RijiX') + RijiX). (5.7) 

Specifically, 

(Standard Model) = y^GpNeSijii, j = 1), (5.8) 

^ \4m\dm) 4m\d.m) J 
(A) = y bhlh^N, + y ^!l!^iVe, (5.10) 

where Nd is the down-quark density in Earth matter. Note that i? is a symmetric matrix and 
also that antineutrinos will have an overall opposite sign for Riy Assuming Earth matter to be 
neutral and isoscalar, = Np = Nn and A^^^ = 2iNf,. 

The current bounds on the ^ couplings [80] imply that the A' induced contributions to 
i?i2 and -Ri3 are several orders less than \/2GFNf,. We neglect those terms in our analysis. The 
upper bounds on all couplings in Rij{X) are also very tight [80] in comparison to v^GiriVe and 
their effect will be discussed later. So, first we consider, in addition to the Standard Model 
contribution, only 

R23 = R32 = ~ (A2ml-^3ml + '^21m^31m) ' 

Am^{dm) 

R22 = ~ {Kill + -^21™) 5 -R33 = „ , 7 , {Kill + Kim) J (5-11) 

Am\dm) Am\dm) 

which are comparable to \/2GpNe. One can see from Eq. 15.111 that R22, 7^ implies both R22 
and i?33 are non-zercl§. 

The current bounds on the relevant couplings follows [80] : 

|A22i,A'23il < 0.18; < 0.06; \K^^\ < 0.58; \K^^\ < 0.52; \K^J < 0.12, (5.12) 

for down squark mass = 100 GeV. The chosen limits on X^im and Kim do not conflict with 
the ratio Rr-w = r(r — >• 7rz/T-)/r(7r /iz/^) [80]. However, the recently published BELLE bound 
on the mode r /ivr^ [134] tightly constrains precisely those products of the A' couplings which 
enter in i?23 = -R32 in Eq. (15.111) . It has been shown that |A2i„A'3\^| and |A2^^A3„^| both must 
be < 1.8 X 10^^( ioo"^ev )^ [135]. This effectively makes -R23 negligible for our purposes. 

In general, it is cumbersome to write an analytical form of the oscillation probability in the 
three-flavour scenario with matter effects. However, under certain reasonable approximations 
it is somewhat tractable. For the energies and baselines under consideration, Am\^L/ E « 1 
and we can use the one mass scale dominance approximation. Under OMSD, the V12 part of 
U (see Eq. I2.17P drops out from Eq. 15.61 In the special case where i?22 = -R33, if one uses the 



^However, in other models [127] this may not be the case. 
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best-fit value of tlie vacuum mixing angle 6'23 = 7r/4 then the neutrino mass squared eigenvalues 
are: 

K22 — -K23, 



2E 



''31 



2E 



Rn + R22 + -R23 =F B 



(5.13) 



where 
B 



Am 



2 \ 2 

31 



2E 



2 /^ 7T7,^ 

+ (— -Rii + -R22 + -R23) — 2— — — 



cos 2^13 {R11-R22-R: 



23 J 



The matter induced neutrino mixing matrix is given by 



1/2 



(5.14) 




(5.15) 



Here 



iVl,: 



2 -i?n + B 



22 



^23 + 



COS 2^13 ±S 



2E 



sm' 



2^ 



13 



-1/2 



(5.16) 



where A^^i {N-i) corresponds to +(-) sign in the above expression. Neglecting the CP phase 5c p 
in the standard parameterization of f/™, one may write = sin 6'^ and = sin 6'^ cos ^]|^. 
From eq. fl5.15p it follows that = 623, the vacuum mixing angle. 6^, on the other hand, 
changes from its vacuum value and it is it /A for 



-Ri 1 — R22 — R 



L23 



^cos2^^i3. 



(5.17) 



In the absence of non-standard interactions, R22 = R23 = R33 = and -Rn = y^GpNg, this is 
the well-known condition for matter induced maximal mixing (see Eq. I2.26p . Since in Eq. 15.151 
^12 ~ 0) to Ufj, oscillation probability the terms involving (M| — Mf ) and (M| 

will not survive and we get : 



M|) 



4 (f/^ 



(5.18) 



where ^'^ii and L are expressed in GeV, eV^, and km, respectively. This expression is 
also valid for antineutrinos. Using Eqs. 15.131 and 15.151 one can easily obtain the oscillation 
probabilities for other channels. 

The above analytical formulation has been used as a cross-check on our numerical results. For 
example. Fig. 15.11 which shows the variation of P^^^^^ as a function of the energy, is obtained 
by solving the full three flavour neutrino propagation equation assuming the PREM [110] 
profile and including NSI. The range of energy is chosen in line with the discussions in the rest 
of this work. The probability falls with decreasing and, for illustration, we have chosen 
sin^ 26'i3 = 0.043 and for all other oscillation parameters we assume the benchmark values 
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Figure 5.1: Py^^y^ for the NH and IH. SM corresponds to only standard electroweak interac- 
tions. The values of A' are given in parentheses, m can take any value, n = 2 or 3. 



given in Table 13.21 Apart from some qualitative remarks, all the results have been presented 
considering 5cp = ^ corresponding to the CP conserving case. The purpose of Fig. 15.11 is 
twofold: (a) to show how the distinguishability between the NH and IH may get blurred by the 
RPVSM interactions, and (b) how irrespective of the hierarchy chosen by Nature the results 
may be completely altered by the presence of these interactions. Each panel of Fig. 15.11 has 
three curves: the solid line (only electroweak interactions), dot-dashed line (in addition, R^'j, 
gets a non-zero RPVSM contribution), and dashed line {R22 = R33 are nonzero, in addition to 
the electroweak contribution). Only in the last case is the analytical formula we have presented 
above applicable. We find excellent agreement. Two aspects of the results are worth pointing 
out. 

First, in the absence of NSI, for the IH with neutrinos, the resonance condition Eq. 15.171 is not 
satisfied and the oscillation probability is negligible (right panel solid line). This could be al- 
tered prominently by the RPVSM interactions (dot-dashed curve) so that the distinguishability 
between the two hierarchy scenarios may well get marred by ^ supersymmetric interactions. 

Secondly, for the NH, it is seen that the peak in the probability may shift to a different energy in 
the presence of the RPVSM interactions. This is because the condition for maximal mixing in 
Eq. 15. 171 is affected by the ^ interactions. For the IH, the oscillation probability is considerably 
enhanced for some energies. Thus, physics expectations for both hierarchies will get affected 
by RPVSM. 

In the following section we dwell on the full impact of this physics on a long baseline CERN-INO 
beta-beam experiment. 
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5.4 Results 



In this work the same magical CERN-INO long baseline beta-beam experimental set-up has 
been considered with which we dealt earlier in chapter 3. Here all the results will be presented 
taking as the candidate ion for the Ue beta-beam|§ source and with 7 = 350. All the 
presented results are based on a five years of ICAL@INC10 data sample. The event rates have 
been estimated using Eq. 13.191 Details about the backgrounds and neutrino-nucleon cross 
sections may be found in section 3.4. 

At the production and detection levels, FCNC and FDNC effects can change the spectrum and 
detection cross sections by a small 0.1%) amount but this would not alter the conclusions. 
At the source and detector, they may also mimic the oscillation signal itself, but these effects 
are tinj@ (~ C»(10-^^)). Here we discuss how FCNC and FDNC may significantly modify the 
propagation of neutrinos through matter over large distances. 



5.4.1 Extraction of ^13 and determination of hierarchy 

If neutrinos have only Standard Model interactions then the expected number of muon events is 
fixeci^l for a particular value of ^13 with either NH or IH as may be seen from the solid lines in Fig. 
15. 2[ The vast difference for the alternate hierarchies picks out such long baseline experiments 
as good laboratories for addressing this open question of the neutrino mass spectrum. 

If non-standard interactions are present then, depending on their coupling strength, the picture 
can change dramatically. In Fig. 15.21 the shaded region corresponds to the allowed values when 
supersymmetric FCNC and FDNC interactions are at play. It is obtained by letting the A' 
coupling^ vary over their entire allowed range - both positive and negative - given in Eq. 15.121 
subject to the further constraints on particular products. 

It is seen that to a significant extent the distinguishability of the two hierarchies is obstructed 
by the ^ interactions unless the number of events is more than about 60. Also, the one-to-one 
correspondence is lost between 6*13 and the number of events and, at best, a lower bound can 
now be placed on ^^13 from the observed number. Of course, if the neutrino mass ordering is 
known from other experiments, then this lower bound can be strengthened, especially for the 
inverted hierarchy. 

^Broadly speaking, the results obtained with ®B neutrino beta-beam taking NH (IH) are similar to that with 
®Li antineutrino beta-beam source considering IH (NH) but details do differ. 

detailed description of the ICAL@INO detector has already been given in section 3.4 and we use the 
same detector characteristics as depicted in Table 13.41 

^This is due to the very tight constraints from — s- e transition limits in atoms [136]. 

^Recall we assume that, but for and the mass hierarchy, the other neutrino mass and mixing parameters 
are known. 

''in fact, we have chosen the subscript m in the A' couplings in Eq. IS.lll to be any one of 1, 2 or 3. 
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Figure 5.2: Number of muon events for the NH (left panel) and IH (right panel) as a function 
of sin^ 2^13 for a five years of ICAL@INO run. The solid lines correspond to the absence of any 
NSI. The hatched area is covered if the A' couplings are varied over their entire allowed range. 



It is also noteworthy that for some values of A'-couplings there may be more events than can 
be expected from the Standard Model interactions, no matter what the value of ^13. Thus, 
observation of more than 112 (5) events for the NH (IH) would be a clear signal of new physics. 



5.4.2 Constraining A' 

If ^13 is determined from other experiments then it will be easier to look for non-standard signals 
from this beta-beam experiment. However, even if the precise value remains unknown at the 
time, considering the upper bound on 613 one may tighten the constraints on the A' couplings. 
Fig. I5.2l reflects the overall sensitivity of the event rate to the ^ interactions obtained by letting 
all RPVSM couplings vary over their entire allowed ranges. In this subsection, we want to be 
more specific and ask how the event rate depends on any chosen A' coupling. 

At the outset, it may be worth recalling that the BELLE bound on r ^ fiir^ [134] severely 
limits the products Ag^^Ag^^ and Ajim-^sim- Thus, R23 can be dropped in the effective neutrino 
mass matrix (see Eq. I5.6p . If only X2mi and/or Ag^^ (-^smi and/or X'^^i^) is non-zero, then 
R22 (-R33) alone receives an RPVSM contribution. Both R22 and i?33 can be simultaneously 
non-zero if \2im and X'^^i (o^ '^2mi and X'^im) are non-zero at the same time. 

In the light of this, we consider the situation where only one of the above ^ coupling is non-zero. 
In such an event, only one of R22, R33 is non-zero. The dependence of the number of events 
on a non-zero A33^ or Ag^^, for a chosen sin^2^i3, can be seen from Fig. 15. 3[ In this figure, we 
use the fact that if only one of these ^ couplings is non-zero, it appears in the results through 
|A'|. For the NH, the curves for Ag^i) for = 2,3, are terminated at the maximum allowed 
value of 0.18. Fig. 15.31 can also be used for A321, X'sim and X2imy bearing in mind their different 
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Figure 5.3: The number of events as a function of a coupling |A'|, present singly, for the NH 
(left panel) and IH (right panel). The thick (thin) lines are for |A33^| (lA'a^j^l, m = 2, 3). The 
chosen sin^ 2^13 are indicated next to the curves. 



upper bounds. For the IH, the number of events is small for y^2m\ ^21^1 insensitive 
to the magnitude of the coupling. These are not shown. It is seen that for the NH there is a 
good chance to determine the ^ couplings from the number of events. In fact, if the number of 
events is less than about 38 there is a disallowed region for |A'|, while for larger numbers there 
is only an upper bound. For the IH, more than about five events will set a lower bound on the 
coupling. 



5.4.3 Effect of A 

The A couplings which can contribute in Eq. 15.101 have strong existing bounds [80] and their 
contribution to R is rather small in comparison to \p2GpN(.. Among them, the bounds A121 < 
0.05 and A321 < 0.07 for m^- = 100 GeV are relatively less stringent [80]. We show their very 
modest impact in Fig. 15.41 It is clear from this figure that (a) the A-type couplings cannot 
seriously deter the extraction of 6'i3 or the determination of the neutrino mass ordering, and 
(b) when ^13 is known in future it will still not be possible to constrain these couplings through 
long baseline experiments. 



5.5 Conclusions 



R-parity violating supersymmetry is among several extensions of the Standard Model crying 
out for experimental verification. The model has flavour diagonal and flavour changing neutral 
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Figure 5.4: Number of muon events for the NH (left panel) and IH (right panel) as a function 
of sin^ 2^13 for a five years of ICAL(QlINO run. The solid lines correspond to the absence of any 
NSI. The hatched area is covered if the A couplings are varied over their entire allowed range. 

currents which can affect neutrino masses and mixing and can leave their imprints in long 
baseline experiments. This is the focus of this work. 

We consider a beta-beam experiment with the source at CERN and the detector at INO. 
We find that the ^ interactions may obstruct a clean extraction of the mixing angle ^13 or 
determination of the mass ordering unless the bounds on the A' couplings are tightened. On 
the other hand, one might be able to see a clean signal of new physics. Here, the long baseline 
comes as a boon over experiments like MINOS which cover shorter distances. Two experiments 
of these contrasting types, taken together, can expose the presence of NSI like RPVSM. 

There are other non-standard models [127] where four-fermion neutrino couplings with greater 
strength have been invoked. The signals we consider will be much enhanced in such cases. 

Our results are presented for the CP conserving case. As 61^ is small, the CP violating effect 
is expected to be suppressed. We have checked this for the Standard Model, where the 'magic' 
nature of the baseline [97] also plays a role. 

Finally, in this work we have restricted ourselves to a beta-beam neutrino source. Much the 
same could be done for antineutrinos as well; then the signs of all terms in R (see Eq. 15. 6p 
will be reversed. It follows from Eq. 15.171 that 6'^ can then be maximal only for the IH and 
as such more events are expected here than in the NH. Broadly, results similar to the ones 
presented here with neutrinos can be obtained with antineutrinos if NH is replaced by IH and 
vice-versa. In the next chapter, we will show that a small detector placed near a beta-beam 
storage ring can probe L violating interactions, as predicted by supersymmetric theories with 
R-parity non-conservation at the production and detection stages of neutrinos. 
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Chapter 6 

Probing Lepton Number Violating 
Interactions with Beta-beam using a 
Near-Detector 



6.1 Introduction 



In the Standard Model, lepton number conservation is only accidental; the particle content 

and the requirement of renormalizability ensure that each lepton flavour number is conserved 
separately. However, non-zero neutrino masses, as indicated by recent neutrino oscillation 
experiments, have proved that the success of the Standard Model should be viewed as that 
of a low energy effective theory. It is not unreasonable to expect that in some extensions of 
the Standard Model, L conservation may not hold. Indeed, a Majorana mass term for the 
neutrinos violates total lepton number. The non-observation of direct L violation in the past 
experiments have put stringent constraints on some of these interactions. In this chapter we 
show that beta-beams and a nearby detector can be a good further probe of such interactions. 
For this work, the advantage of a 'near' detector is twofold. 

1. Firstly, due to the short base-length, neutrinos do not get much scope to 
oscillate before being detected, which could otherwise mimic signals of the 
L- violating {p} interactions. 

2. The other obvious advantage is that a Icirger pcirt of the beam can be picked 
up with a smaller detector. 

We consider placing a 5 kton cylindrical detector, aligned with the beam axis, within 1 km 
from the beta-beam storage ring. The // interactions can lead to tau leptons in near-detector 
experiments in two ways. A Vr can be produced due to such interactions during beta-decay, 
yielding a r through weak CC interactions in the detector. Alternatively, the electron neutrinos 
in the beam, produced through usual beta-decay, can undergo ^ interactions with the detector, 
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Figure 6.1: A schematic diagram of the proposed detector (a part only). The incoming z/g beam 
may have a very small contamination of neutrinos of other flavours in the presence of lepton 
flavour violating interactions. 



leading to tau leptons. The taus promptly decay, part of the time in a muonic channel. Iron 
calorimeters with active detector elements serve well for identifying these muons, which leave 
long tracks in the detector, and for filtering out backgrounds. We will also briefly comment on 
water Cerenkov and other detectors. 

In the following section a brief account of the experimental set-up is presented. In section ESI 
interactions are discussed in the context of R-parity violating supersymmetric models [80,137]. 
We stress how beta-decay can be affected in the presence of such interactions, yielding Vr in 
a few cases in place of the standard z/g. The processes via which produce tau leptons in 
the detector are also described. The expected number of muon events from tau decay and the 
constraints ensuing in the event of their non-observation will be presented in section 16. 4[ 



6.2 Beta-beam Flux at a Near-Detector 



The choice of ions for a beta-beam is predicated by the intended physics. The low end-point 
energies (see Table II. ip of the ®He and ^^Ne ions restrict the energy reach of the beam; a 
threshold energy of 3.5 GeV is necessary to produce a r-lepton from an incoming neutrino. 
Therefore, we will work with and ^Li ions, offering higher end-point energies. We pick only 
the neutrino beam for our discussion. 

The geometry of the beta-beam storage ring determines the neutrino flux at a near-detector. 
For a I0W-7 design, a 6880 m decay ring with straight sections of length (= S) 2500 m each 
(36% useful length for ion decays) has been proposed. In such a configuration, Np = 1.1 x 10^^ 
useful decays (decays in one of the straight sections) per year can be obtained with ^^Ne ions. 
We have used this same luminosity for *B and higher 7. To settle these issues a dedicated study 
is on at CERN. 
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6.2.1 Detector Simulation Study 



We consider a cylindrical 5 kton detector (as in Fig. 16.11) aligned with one of the straight sections 
of the storage ring. The detector is made of iron slabs (thickness 8 cm) with interleaved active 
detector elements (thickness 2.5 cm) such as resistive plate chambers (RPCs). The readouts 
from these RPCs will be concentric annular strips of small width with further segmentation 
to improve the position resolution. In this proposal, iron is the main content of the detectoi0. 
The thickness of the slabs ensures that electrons do not propagate in the detector. The signal 
muons are of sufficient energy to give rise to long tracks. To eliminate possible beam-induced 
backgrounds (see below) from pions produced in CC and NC processes, typically 6 to 13 hits 
(depending on the boost 7 ranging from 250 to 450) are required of a putative muon track. 

As noted earlier, the signature of new physics we consider is the appearance of prompt tau 
leptons which decay into muons with a branching fraction of 17.36% [138]. The tau production 
threshold is around 3.5 GeV. This is what necessitates the higher boost 7. 

Backgrounds, other than those of the beam-induced variety discussed below, are controllable, 
as we now point out. A beam-off run will help make a first estimate of these backgrounds. 
Further, an important aspect of the beta-beam source is its capability of eliminating back- 
grounds through timing information. The beam itself will consist of bunches of typically 10ns 
size and the number of bunches will be chosen so as to ensure that the ratio of the active- 
to the dead-time is (9(10^^). Backgrounds from other sources, namely, atmospheric neutrinos, 
spallation neutrons, cosmic rays, etc. can thus be largely rejected from the time-stamp of a 
recorded event. Even further reductions of the backgrounds of external origin can be envisioned 
through fiducial and directionality cuts. 

Now let us turn our attention to the issue of beam-induced backgrounds caused by NC and CC 
interactions of unoscillated v^. Electrons produced through weak interactions by the incoming 
Vf. are quickly absorbed and do not leave any track. Formation of prompt muons through ^ 
supersymmetric interactions is suppressed by strong bounds on the relevant couplings arising 
from limits on /i — e transitions in atoms [136]. However, the beta-beam neutrinos can produce 
pions along with other hadrons at the detector via CC and NC processes. They undergo strong 
interactions with the detector material and are quickly absorbed before they can decay. But as 
numerous pions are produced, it needs to be checked whether some of them can fake the signal. 

We have checked our naive expectations with a detector simulation study using GEANT [139] 
aided by NUANCE [140]. It is observed that for neutrino- nucleon interactions at energies 
interesting for our study, the produced lepton preferentially carries most of the energy of the 
incident neutrino. Moreover, pions are usually produced with multiplicity more than unity. 
Hence it is not unreasonable to expect that the pions will be less energetic compared to the 
taus produced via fj interactions and hence in detectors of this genre, it is possible to distinguish 
hadronic showers from a muon track. 

"'^Lead may be an interesting alternative material to enhance the event rate. 



88 



95 




90 - N 



y = 250 (6 hits) 

Y =350 (10 hits) 

y = 450 ( 13 hits) 



65 



20 40 60 80 100 120 140 160 180 200 220 



Detector length (m) 



Figure 6.2: Detector efficiency for 7 = 250, 350, 450. The corresponding cuts on muon hits used 
are 6, 10 and 13 respectively. 



However, a conservative approach is followed in pion background estimation. Although pions do 
not leave behind a straight track like a muon, we still count the number of hits as a measure of 
the distance traversed by a pion. We impose a criterion of minimum number of hits to identify 
a track to be a muon one. We find from a simulation that, for 7 = 250/350/450, imposing a 
cut of 6/10/13 hits will reduce the pion background at least to the 10"'^ level. 

The detector geometry plays a role in determining the signal efficiency after imposition of these 
cuts. Since the muons produced from boosted tau lepton decay carry transverse momentum, 
some of them may exit the detector through its sides, failing to satisfy the cuts. For a fixed 
detector mass (5 kton), a longer detector has a smaller cross-sectional area, resulting in a drop 
in the detector efficiency for the above reason. As the detector length increases from 20 m 
to 200 m, with our set of cuts, the efficiency factor reduces from 85% to 70% approximately, 
showing little dependence on 7 (see Fig. 16. 2p . 

6.2.2 Neutrino Fluxes 

Neglecting small Coulomb interactions, the lab frame neutrino beta-beam flux (per unit solid 
angle per unit energy bin per unit time per unit length of the straight section) emitted at an 
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angle 9 with the beam axis is described by (see also Eq. I3.2p 



^Near^j^^Q^ = - ^ (^E, - E*)E*^ ^ {E, - E^f - ml (6.1) 

47r J 7(1 — p cost^) 

where g = Np/S is the number of useful decays per unit time per unit length of the straight 
section. 

To calculate the resulting number of events at a cylindrical near-detector of radius R and length 
D aligned with the beam axis it is necessary to integrate over the length S of the straight section 
of the storage ring and the volume of the detector. The event rate at a detector placed at a 
distance L from the storage ring is given by [141] 



pS pD i-e' pE- 
ne dx d£ rf^27rsin^ / dE (f)^^''''{E,e) a{E), (6.2) 

Jo Jo Jo JE^i„ 

tanO' (x,£) = - — — and E' = .^^ (6-3) 

^ ' L^x^i 7(1 -/5 cos 6*) ^ ' 

Here n represents the number of target nucleons per unit detector volume, e is the detector 
efficiency as presented in Fig. 16.21 Emin denotes the tau production threshold, and (y{E') stands 
for the neutrino-nucleon cross section. Note that the source of L-violation may lie either in 
(pNear^^j^^ 6*) (in case of ^ beta-decay) or in <7{E) (in case of ^ tau production). 

To help subsequent discussion, following [141], we rewrite the above formula isolating the 
geometry integrated total flux $(£'; S,D,R, L) (per unit time per unit energy bin) falling on 
the detector and emitted from the whole length of the straight section as follows : 



dt 

where 



dN^ 
dt 



ne I """^ dE^{E;S,D,R,L)a{E), (6.4) 

Em. in 



where 



pS pD pa- 

<^>{E-S,D,R,L)= dx d£ d9 2n sin 9 ^^^""'{E, 9) (6.5) 
Jo Jo Jo 

Eq — rrie 

= ^(13^5) • ^'-'^ 

The beta-beam also involves a few small uncertainties which we neglect in our analysis. However 
for completeness, these are listed here : 



and 



• There exist different excited states of the daughter nuclei of the decaying ion, which 
additionally lead to small contributions to the spectra with different end-point energies. 

• The ion beam has a finite transverse size. However, as this size varies [142] between only 
3.0 cm to 5.1 cm, with an average of 4 cm {3a), in both transverse directions inside the 
ring, the variation in flux at the detector due to this is negligible. 



90 




Figure 6.3: Feynman diagrams for ^ interactions during beta-decay through (a) A'A' and (b) 
AA' type trihnear product couphngs. Substantial event rates are obtained in (a) when k = 2,3. 



• The decaying ions may have small transverse momentum due to thermal fluctuations 
[ksT ~ 2.6 X 10"'^ eV), but this can be safely ignored in comparison with the end-point 
energy of the beta-decay. 



6.3 ft Processes 



A detailed description of the RPVSM [80, 137] containing renormalizable trilinear A- and 
A'-type couplings is given in section 5.2. To impose conservative upper bounds, we work in 
a minimal ^ framework where only a pair of such couplings are assumed to be non-zero at a 
time. For a near-detector, ^ interactions can come into effect in two ways as described in the 
following subsections. 



6.3.1 Ur Production in Beta-decay via ^ Interactions 

I? interactions can drive beta-decay producing Ur instead of z/g. so produced give rise to r 
leptons in the detector which may decay in the leptonic channel producing muons. 

Simultaneous presence of X'^^j^ and A'^^^^ couplings (see Fig. 16.3( a)) can be responsible for pro- 
ducing a in beta-decay. Of these, X^^^ is tightly constrained from neutrinoless double 
beta-decay [143]. But the upper bound on the combination \X'^ij^X[ii.\, k = 2,3 is rather re- 
laxed; a limit of 2.4 x 10^^(m/100 GeV)^, m being a common sfermion mass, follows from 
r~ e~p° [136,138]. The corresponding decay amplitude can be written as. 

Alternatively, z/^ can be produced in beta-decay if another combination of ^ couplings X*^iX'm 
{i = 1, 2) is non-zero (see Fig. 16.3( b)). As mentioned earlier, X'm is severely constrained. The 
combination |A23iA2]^^| is bounded from above by 1.6 x 10^^(m/100 GeV)^ arising from the 
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Figure 6.4: Feynman diagrams for tau production from an incoming Ve beta-beam through (a) 
A' A' and (b) AA' type trihnear product couphngs. Substantial event rates are obtained in (a) 
when = 2, 3. 



decay channel r —>■ e rj^ [136, 138], which is not too small to produce an observable effect. 
The corresponding decay amplitude is given by, 

M^{u^de+u,) = ^hl^[u^^Pj,u,][udPLUu]. (6.8) 

[t — m^j 

6.3.2 r production from via ^ Interactions 

Ue produced through ordinary beta-decay driven by weak interactions can undergo ^ interac- 
tions with the detector producing r which subsequently decay into muons. 

Simultaneous presence of Ag^^ and A'^^^j^. couplings can give rise to r~ in the final state from 
an incoming Ue beta-beam (see Fig. 16.4( a)). The amplitude for the corresponding s-channel 
diagram can be written, after a Fierz transformation, as 

M^{u,d-^r-u) = [Urlf^PLUu^ [uurPLUd]. (6.9) 

An alternative channel of tau production from an incoming Ue beam exists (see Fig. 16.4( b)) if 
a particular combination of the A and A' couplings AjisA'^^^ {i = 2, 3) is non-zero. Here again, 
A313 is severely constrained from neutrinoless double beta-decay experiments [144]. An upper 
bound of 1.6 x 10^'^(m/100 GeV)^ applies to the combination |A2i3A2\^|, from the decay channel 
T~ e^rf [136,138]. The amplitude for this t-channel process is 

Mf^{u, d^T- u) = [urPLU,:^ [UuPnUd]. (6.10) 

(t — m^j 

In what follows, we categorize the above two kinds of diagrams (a) and (b) in Figs. 6.3 and 6.4 
as A'A' and AA' processes, respectively. 
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Note that, if l^'iik^'ukl^ /c = 2, 3 is non-zero, tau leptons can be produced at the detector either 
due to ^ interactions affecting beta-decay or due to ^ interactions of a u^, with the detector 
material. These two equal contributions add in the total rate of tau production. 

However, for the A A' process, we see that the ^ combinations IA231A211I (which drive the ^ 
beta-decay) and |A2i3A2\^| (which is responsible for producing a tau from an incoming z/g in 
the detector) are different. As we are following the strategy of taking only two ^ couplings 
non-zero at a time, these contributions, which are of the same magnitude, cannot be present 
at the same time. 

In passing, a few comments are in order : 

• In both the diagrams of Fig. 16.41 the incoming can interact with a u quark from the sea 
to produce a tau. Due to the smallness of the corresponding parton distribution function, 
this contribution is suppressed but we do include it in the numerical evaluations. 

• Here we should mention that the FCNC process, —>■ n^uu [145] puts stringent bounds 
on all the A' couplings. However, these are basis dependent and hence can be evaded. 

• As already noted, the non-observation of the process /i — > e (Ti) severely restricts [80,136] 
the possibility of emitting a in beta-decay and direct production of muons from an 
incoming z/g beam. 

• Since the beta-beam energy is ~ a few GeV, the expected event rate will be essentially 
independent of the sfermion mass as the bounds on A, A' scale with (m/100 GeV)^. 

At this energy range it is important to consider contributions from deep-inelastic, quasi-elastic, 
and single-pion production channels. To estimate the ^ deep-inelastic scattering cross section, 
we have used CTEQ4LQ parton distributions [146]. ^ quasi-elastic scattering and single-pion 
production cross sections have been evaluated from the corresponding Standard Model cross 
section^ [120,147] by a rescaling of the couplings. It is noticed that, as Eq. I6.10l is not reducible 
to a Standard Model-like {V — A)^{V — A) Lorentz structure, in calculating deep-inelastic cross 
section a factor ~ 1/3 appears from polar integration compared to that for Eq. 16.91 For the 
AA' process we have adopted the same suppression factor for the Standard Model quasi-elastic 
and single-pion production cross sections as well. Conservatively, we assume that a similar 
suppression also applies to the case of ^ beta-decay. It bears stressing that the effect of the 
tau mass is felt on the neutrino-nucleus cross section throughout the energy range beyond the 
T-threshold and this is included in the analysis. 



6.4 Results 



A near-detector set-up is qualitatively different from a far-detector as in the former case the 
storage ring and the detector really 'sees' each other and relative geometric considerations are 



^ These cross sections include all nuclear effects for an iron target. 



93 



of much relevance. The observed number of events in a given period of time depends on the 
choice of the radioactive ion, the boost factor 7 and the details of the set-up (which include 
storage ring parameters, detector configuration and the short base- length between them). As 
alluded to earlier, the maximum 7 available is limited by the storage ring configuration. With 
a view to optimizing the set-up, the essential inputs are summarized as follows : 



Storage Ring Parameters : 



Total length 6880 m, length of a straight section, S = 2500 m, number of beta-decays in 
the straight section, A^^ = 1.1 x 10^*^ per year. 



Detector Configuration : 



The detector material iqj iron (p = 7.87 gm/ cm^). A detector of mass 5 kton is considered. 
For a given material, this fixes the length of the detector as the radius is changed. It 
varies from 202.13 m to 12.63 m as the radius ranges over 1 m to 4 m. 



Base-length : 



Results are presented for three representative values of the distance of the detector from 
the storage ring, L = 200 m, 500 m, 1 km. 



Boost factor 7 : 



The tau production threshold (3.5 GeV) calls for a high 7. We consider 7 = 250, 350, 450 
for and as large as 800 for ^^Ne with fixed Njj. 

The high collimation achievable in the beta-beams encourages the choice of a detector of cylin- 
drical shape coaxial with the storage ring straight section. As 7 increases, the ^ event rates 
increase for the following reasons : 

1. an increasingly larger part of the beam falls onto the detector, 

2. more neutrinos have enough energy to produce a tau lepton, 

3. with the more energetic neutrinos the cross section is larger. 



The first two effects are demonstrated by Fig. 16.51 The geometry integrated flux, $, as defined 
in Eq. 16.51 represents the beta-beam neutrino flux spectrum falling onto the detector per unit 
time. It is seen that as 7 increases, the total area under the curve also increases, illustrating 
the first effect. The area under the curves on the right side of the vertical line (the threshold) 
also increases with 7, in conformity with the second expectation. For a high 7 the beam should 
saturate. However, with the 7 used in Fig. 16.51 this is not evident due to the enormous length 
of the straight section of the storage ring. To coUimate the flux emanating from the rear part 
of the ring a very high 7 will be needed. 

•^Brief comments are made about a water Cerenkov detector in section 6.4.2. 
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E [GeV] 

Figure 6.5: Geometry integrated flux 5*, D, R,L)/g taking as the decaying ion is plotted 
against neutrino energy E for different 7 for S = 2500 m, D = 202.13 m, i? = 1 m, and L = 200 
m. The vertical line at 3.5 GeV indicates the tau production threshold energy. 



As geometry plays a crucial role in optimizing the near-detector set-up, we study the detector 
length dependence of the expected number of ^ events for different base-lengths and different 7. 
We consider the contribution coming from the two options - the A'A' and AA' processes - in dif- 
ferent panels for every figure, assuming the ^ coupling constants saturate present experimental 
upper limits. 

6.4.1 Choice of Ion Source and Detector 

The choice of as the ion source provides the most attractive option due to its high end-point 
energy. Iron calorimeters are preferred for the smaller size and significant background removal. 

To get a glimpse of the number of events one might expect in such a set-up, let us present the 
following estimate. A 5 kton Fe detector of radius 1 m (length 202.13 m) placed at a distance 
200 m from the decay ring can give rise to 92 (24) muon events via the A'A' (AA') process in 5 
years foi0 7 = 250. 

■^The corresponding numbers for 7 — 350 are 421 (103). 
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Detector length (m) Detector length (m) 



Figure 6.6: Expected number of ^ muon events in five years for a 5 kton iron detector vs. the 
detector length for 7 = 250, 350, and 450 for beta-beam flux. The left (right) panel is for 
the A' A' (AA') driven process, /c = 2, 3. 

In Fig. I6.6l we exhibit the 7 dependence of the expected number of muon events over a five-year 
period for a fixed base-length of 200 m. Collimation plays a role as is demonstrated by the 
increase in the number of events for higher 7. As expected, a long detector serves better as it 
provides more opportunity for a neutrino interaction to occur. However, this increase with the 
length is not linear; a part of the beam is lost due to the concomitant decrease in the radius 
(to keep the total mass fixed at 5 kton). In addition, with the increase in detector length as 
the detector efficiency decreases, the increase in the rates is also somewhat restricted. 

It is also of interest to study the base-length dependence of the number of events. The beam 
spreads with an increase in the base-length, reducing the effective flux hitting the detector. 
This causes a fall in the number of events (other parameters remaining the same) as shown in 
Fig. 16.71 It is interesting to note that the increase in the number of events with increase in the 
length of the detector gets severely diluted at larger base-lengths. 

While presenting the expected number of events we assumed the ^ couplings saturate the 
present experimental upper bounds. In case less or even no events are seen, the existing limits 
on the combinations jAg'^^A']^^^!, A; = 2,3, IA231A211I and |A2i3A2*i]^| will be improved. Choosing 
the minimum number of non-zero ^ couplings, one can put conservative upper bounds. In Fig. 
16.81 we show the bounds - the region above the curves are disallowed - achievable in the case 
of 'no-showj^l. It is seen that to put stringent bounds it is necessary to go for a higher 7 and a 
longer detector. 

^At 95% C.L. this corresponds to not more than 3 events. 



96 



u 



a 

0) 

> 



450 
400 
350 
300 
250 
200 
150 
100 
50 



— \ \ 1 \ 1 1 \ — 

l4lk^'llk'= •'■•'"24 200 m, 



1=350 



500 m 



1km 



20 40 60 80 100 120 140 160 180 200 220 

Detector length (m) 



0) 

•H 



no 

100 
90 
80 
70 
60 
50 
40 
30 
20 
10 



1 — ^ 

.1 A,213 A,211 1 = 0.0016 



IA,23i A,2iil = 0.0016 



1=350 




20 40 60 80 100 120 140 160 180 200 220 

Detector length (m) 



Figure 6.7: Muon signal event rate in 5 years as a function of the detector (Fe) length for three 
different choices of base-length have been shown for beta-beam flux. The left (right) panel 
corresponds to the A'A' (AA') driven process, /c = 2, 3. 



6.4.2 Alternative Set-ups 

Although so far we have presented results with as the beta-beam source, ^^Ne is the most 
discussed decaying ion in the literature. As mentioned earlier, due to the smaller end-point 
energy of ^^Ne, a high 7 is required to cross the r threshold. Fig. I6.9l depicts the variation in the 
expected event rate with detector length for ^^Ne with 7 = 800 using a 5 kton iron calorimeter. 
We see that due to high 7 for ^^Ne, the beam is so collimated that the event rates increase 
almost linearly with increasing detector length in contrast to the ^B case we have presented. 
However even in such an extreme scenario, where we use the same storage ring configuration 
to reach such a high 7, the expected event rates are comparable to that in the ^B case. Hence 
we conclude that ®B is preferred to ^^Ne in exploring lepton number violating signatures with 
beta-beams. 

The use of water Cerenkov detectors with good capability of muon-electron separation and 
moderate efficiency of neutral current rejection may be an interesting option to see the signals 
of new physics and to normalize the incoming flux. The disadvantage of this set-up turns out 
to be the huge background. Consider a 5 kton water Cerenkov detector with radius 2.5 m at a 
distance 200 m from the decay ring. In five years, this will lead to 45 (12) muon events from 
r-lepton decay for A'A' (AA') driven processes from an incoming ^B z/g beam accelerated with 
a 7 of 250 and with a muon detection threshold of 200 MeV. For the same configuration and 
duration, one expects roughly 10^ pions produced from charged and neutral current interactions 
of the z/e beam. Muons produced from vr decay will thus completely swamp the signal. 

The number of signal events may be increased by designing a very long water detector with 
small radius though this could be technologically challenging. In any case, the background 
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events will continue to be very high. So, this option also does not hold much promise. The 
basic problem of high backgrounds, avoided in the Fe detector, will also plague totally active 
scintillator based detectors. 



6.5 Discussion and Conclusion 



Beta-beam experiments may be sensitive to the ^ interactions. In the previous chapter it was 
shown that ^ interactions can interfere with pure oscillation signals in long-baseline beta-beam 
experiments. In this chapter we explore a complementary scenario. It has been proposed that 
to probe such interactions, an iron calorimeter detector placed close to the storage ring holds 
promise as it provides essentially a neutrino oscillation free environment. In particular, the 
combination of a 5 kton cylindrical iron detector placed within a distance of 200 m to 1 km 
from the decay ring and a neutrino beam from an ^B ion source with 7 in the range 250 to 450, 
running for 5 years is well-suited in this regard. We have examined the impact of non-trivial 
design details of such a near-detector setup. 

At production, low energy beta-decay experiments may get contaminated by tau neutrinos 
through ^ interactions. We show that, this contamination, though small, can be probed us- 
ing the above set-up. 1^ interactions can also play a role in such an experiment during the 
interactions of the beta-beam electron neutrinos with the detector. 

It is interesting to explore if 1^ interactions can affect beta-beam experiments in other ways. 
For example, we have checked that the impact of these interactions on the 11 detection cross 
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Figure 6.9: Comparison of the muon signal event rates as a function of the detector length for a 5 
kton iron calorimeter placed at a distance of 200 m from the storage ring for 7 = 800 (250) with 
i8]\jg (^B). The left and right panels correspond to A' A' and AA' driven processes, respectively. 



section is insignificant. As mentioned earlier, u^j, may be produced in beta-decay through ^ 
interactions but this also is severely suppressed as the corresponding couplings have stringent 
upper limits. 

Results are presented for a neutrino beam. Antineutrino beams can also be produced using 
^Li or ^He as sources. In fact, a storage ring design may allow both beams to be present 
simultaneously. The expected event rates for antineutrinos are of similar order as for the 
neutrinos. 

In conclusion, we find a near-detector setup can be useful for exploring lepton number violating 
interactions with beta-beams. It may allow us to put stringent bounds on some of these 
couplings. In the next chapter, we will present a summary of this doctoral work. 
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Chapter 7 

Summary and Conclusions 



Neutrino physics is traversing through an exciting phase with lots of new data making this field 
more interesting. There is an excellent progress in our understanding of the neutrinos over the 
last ten years or so, thanks to the experiments on neutrino oscillations which confirm the fact 
that neutrinos have a tiny, but non-zero, mass quite against the expectations of our best theory 
- the Standard Model. We have now a rough picture of the parameters governing three-fiavour 
oscillations and we are all set to move into the precision regime. There is no doubt that the 
use of artificial neutrino sources is mandatory in the era of high precision experiments. In 
this direction the beta-beam is a recently proposed technique of producing a pure, intense and 
collimated beam of or Vf, through the beta-decay of completely ionized radioactive ions. My 
thesis sheds light on the importance of beta-beams in probing some unknown territories in the 
field of neutrino physics. 

We started with a brief description of neutrino properties in the first half of chapter 1 and then 
we took a glance at the fascinating journey of discovery into one of Nature's most impalpable 
particles - the neutrino which is the main theme of this thesis. In the second half of chapter 1 
we focused on the main sources of neutrinos. 

The first part of chapter 2 dealt with a basic introduction to the quantum mechanics of neutrino 
oscillation in vacuum under both two and three fiavour frameworks. Then the importance of 
matter effects in neutrino oscillations was discussed. In the second half of chapter 2 we took a 
look at our present understanding of neutrino parameters and we identified the major unknowns 
in the neutrino sector. Finally we closed chapter 2 by giving a brief note on the future neutrino 
road-map based on long baseline experiments. 

In chapter 3 we expounded in detail the physics reach of an experimental set-up in which the 
proposed large magnetized iron detector at the INO would serve as the far detector for a beta- 
beam. If this pure and/or Pg beam is shot from some source location like CERN such that 
the source-detector distance L ~ 7500 km, the impact of the CP phase 5c p on the oscillation 
probability and associated parameter correlation and degeneracies are almost negligible. This 
"magical" long baseline beta-beam experiment would have unprecedented sensitivity to the 
neutrino mass hierarchy and ^13, two of the missing ingredients needed for our understanding 
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of the neutrino sector. With Lorentz boost 7 = 650 and irrespective of the true value of 6cp, 
the neutrino mass hierarchy could be determined at 3a C.L. if sin^ 2^13 (true) > 5.6 x 10~^ 
and we can expect an unambiguous signal for ^13 at 3a C.L. if sin^ 2^13 (true) > 5.1 x 10~^ 
independent of the true neutrino mass ordering. 

In chapter 4 we showed that the earth matter effects in the i/g ~^ survival probability can be 
used to cleanly determine the third leptonic mixing angle ^13 and the sign of the atmospheric 
neutrino mass squared difference, Amg^, using a beta-beam as a source. 

In chapter 5 we discussed possible new physics scenarios with a beta-beam neutrino source 
and made a strong physics case for these high-precision experiments. Long baseline neutrino 
oscillation experiments may well emerge as test beds for neutrino interactions as are present 
in R-parity violating supersymmetry. We showed that flavour diagonal and flavour changing 
neutral currents arising therefrom prominently impact a neutrino beta-beam experiment with 
the source at CERN and the detector at the proposed India-based Neutrino Observatory. These 
interactions may preclude any improvement of the present limit on ^13 and cloud the hierarchy 
determination unless the upper bounds on ^ couplings, particularly A', become signiflcantly 
tighter. If ^ interactions are independently established then from the event rate a lower bound 
on 6*13 may be set. It has been shown that there is scope to see a clear signal of non-standard 
FCNC and FDNC interactions, particularly in the IH scenario and also sometimes for the NH. 
In favourable cases, it may be possible to set lower and upper bounds on A' couphngs. FCNC 
and FDNC interactions due to A type ^ couplings are unimportant. 

A detector placed near a beta-beam storage ring can probe ^ interactions, as predicted by 
supersymmetric theories with R-parity non-conservation. This issue has been discussed in 
great detail in chapter 6. In the presence of ^ interactions, z/^ can be produced during beta- 
decay leading to tau leptons through weak interactions. Alternatively, electron neutrinos from 
beta-decay of radioactive ions can produce tau leptons in a nearby detector through these 
interactions. The muons from the decay of these tau leptons can be readily identified in a small 
iron calorimeter detector and will signal violation of R-parity. 

The study of neutrinos has always landed up with lots of surprises. So we can expect further 
surprises in store and, obviously, beta-beams can play a leading role in this direction. This 
doctoral work is an effort to judge the expected performance of a beta-beam neutrino source 
in the future progress of neutrino physics and also for the hunt for signals of non-standard new 
physics. We hope that this work will provide a boost to the detailed and thorough R&D of the 
novel beta-beam neutrino source in future. 
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